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STRUCTURAL  DYNAMIC  THEORIES  OF  FRACTURE' DAMAGE  OIAGNOSTS 

ABSTRACT 

In  this  report,  two  modal  theories  of  damage  diagnosis  in  simple  and 
complex  structures  are  presented.  The  modal  frequency  theory,  based  on  the 
changes  in  the  modal  frequencies  of  a  structure  due  to  damage,  is  shown  to 
be  a  suitable  method  for  simple  fracture  configurations,  but  not  for  complex 
ones.  The  transmi ssi bi  1  i ty  theory,  which  utilizes  the  changes  in  trans- 
missibi  lities  across  a  structure,  is,  on  the  other  hand,  proposed  as  a 
feasible  method  for  damage  diagnosis  in  large  structures. 

The  modal  frequency  theory  is  applied  to  multiple  cracks  in  simple 
beams  to  demonstrate  its  intrinsic  uncertainties.  It  is  shown  that  closely 
spaced  multiple  cracks  behave  as  one  single  crack  as  far  as  the  modal 
frequencies  are  concerned.  Uncertainties  may  arise  in  the  diagnosis  of 
multiple  cracks  when  an  inadequate  number  of  frequency  measurements  are 
available.  When  the  structure  has  one  major  crack  in  addition  to  minor 
ones,  hcwever,  the  major  crack  is  diagnosable. 

In  the  transmi ssi bi 1 i ty  theory,  the  effects  of  (1)  the  excitation 
location  and  frequency,  (2)  the  locations  of  the  response  stations  where 
transmi  ssi  bi  1  i  ti  es  are  computed,  (3)  the  crack  location  and  severity,  and 
(4)  damping  on  the  transmissibility  changes  are  studied  for  a  three-story 
four-span  frame  structure.  It  is  found  that  best  results  are  obtained  when 
the  frame  is  excited  near  the  joints  at  its  modal  frequencies  and  when  the 
response  stations  are  located  in  the  vicinity  of  the  minimum  deflection 
points.  The  largest  transmissibility  changes  are  shown  to  occur  at  the 
response  stations  nearest  the  crack.  This  result  allows  the  crack  to  be 


located  within  one  or  two  cells  of  a  frame  structure.  The  regions  of  the 
frame  which,  when  cracked,  are  undiagnosable  at  the  fundamental  frequency 
are  established  in  the  study.  The  advantage  of  exciting  the  structure  at  a 
higher  modal  frequency,  in  such  cases,  is  demonstrated.  It  is  shown  that 
undiagnosable  regions  become  larger  for  less  severe  cracks. 

The  structures  are  modeled  without  loss  of  generality  as  composed  of 
generalized  Bernoulli-Euler  beams,  undamped  in  the  application  of  the  modal 
frequency  theory  and  damped  in  that  of  the  transmissibility  theory.  Damping 
is  incorporated  by  means  of  the  complex  Young's  modulus.  In  the  analysis, 
the  general  theory  of  circuit  is  utilized  to  represent  simple  beams  with 
basic  electric  circuits,  and  the  fracture  hinge  by  an  electrical  resistor. 
The  damage  analysis  of  a  complex  frame  structure  is  thus  formalized  and 
adaptive  for  computer  programming.  A  formalized  program  for  the  computation 
of  frame  deformation  is  enclosed. 
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CHAPTER  1 

INTRODUCTION 

Fracture  damage  diagnosis  stems  from  the  need  to  ensure  safety  of 
structural  components  or  structures,  or  from  the  requirements  of  quality 
control  in  production.  Structures  may  develop  cracks  from  aging  but  most 
likely  as  a  result  of  some  strong  excitations  such  as  in  earthquake,  blast 
loading,  wind  loading  and  the  like.  Diagnosis  of  flaws  is  thus  essential  in 
assessing  the  reliability  and  integrity  of  the  structure.  The  diagnostic 
method  must  be  of  such  a  nature  that  it  is  nondestructive.  The  present 
report  summarizes  two  diagnostic  theories  for  fracture  damage  in  structures 
based  on  the  modal  theory  of  structure  dynamics.  The  first  one  is  the  modal 
frequency  theory.  The  second  one  is  the  transmissibil ity  theory. 

Nondestructive  testing  began  gaining  importance  in  industry  after  World 
War  I  even  though  there  were  earlier  studies  to  detect  cracks  in  metals 
nondestructi vely  [I],  Nondestructive  testing  had  no  great  impetus,  however, 
until  World  War  II.  Since  then,  it  has  become  a  major  area  of  research  and 
development.  Nondestructive  test  techniques  include  the  categories  of 
radiography,  sonic-ultrasonic  methods,  visual  methods-penetrants,  electro¬ 
magnetic  methods  and  thermal  methods.  The  current  ASME  code  for  the  inspec¬ 
tion  of  boilers  and  pressure  vessels  recognizes  some  of  the  existing  methods 
which  fall  into  the  first  four  of  these  categories  [2].  None  of  these 
methods,  however,  yield  good  results  when  the  surfaces  of  the  crack  are 
nearly  parallel  to  the  direction  of  propagation  of  waves  or  magnetic  flux 
lines.  On  the  other  hand,  all  of  these  methods  are  local  methods,  that  is, 
the  excitation  (X-rays,  ultrasonic  transmitter,  magnetic  coil,  or  heat  flux) 
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and  the  sensor  must  be  physically  near  the  flaw.  As  a  result,  a  scanning  of 
the  structure  or  the  component  is  necessary  to  cover  all  areas  which  are 
likely  to  contain  flaws.  A  global  full-scale  testing  method  is,  therefore, 
needed  for  diagnosis  of  damage  in  large  structures  such  as  bridges,  offshore 
platforms,  buildings  and  the  likes.  By  a  global  method  is  meant  a  method 
whereby  the  location  and  possibly  the  severity  of  flaws  in  a  structure  can 
be  determined  by  exciting  the  structure  at  an  arbitrary  point  and  monitoring 
the  responses  at  various  stations  on  the  structure,  without  the  need  for 
scanning.  An  ultrasonic  technique  described  as  global  was  proposed  [3,4] 
for  the  inspection  of  tubular  K-joints,  where  globalism  is  defined  as  the 
capability  to  monitor  an  entire  joint  rather  than  an  entire  structure.  The 
proposed  method,  however,  is  not  capable  of  detecting  cracks  which  are 
located  around  the  periphery  of  the  major  tubular  column  near  the  joint  [4]. 

A  global  method  of  damage  detection  studied  recently  Is  the  random 
decrement  technique  [5,6].  The  random  decrement  process  is  a  signal  pro¬ 
cessing  technique  which  extracts  the  free  decay  responses  from  the  random 
dynamic  responses.  The  changes  in  the  free  decay  responses  at  several 
locations  on  an  experimental  scale  model  of  an  offshore  platform  were  used 
to  detect  the  damage  in  the  structure  [5].  The  scale  model  consisted  of  six 
levels  with  K-joints  at  the  mid-span  of  each  girder.  The  predictions  of  the 
damage  locations  in  various  damage  scenarios  were  not  very  specific  with  the 
random  decrement  method.  The  cross  random  decrement  method  was  later 
employed  as  an  improvement  [6].  This  technique  utilizes  two  simultaneous 
response  measurements  from  two  different  positions  on  the  structure. 
Experiments  were  carried  out  on  a  similar  platform  model .  Relative  phase 
shifts  between  every  two  neighboring  response  locations  were  computed  using 
the  cross  random  decrement  technique  before  and  after  damage  was  introduced. 
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The  largest  changes  in  the  relative  phase  shifts,  in  general,  occurred  at 
locations  near  the  damage.  Both  studies  by  Yang  et  al.  [5]  and  Tsai  et  al . 
[6]  were  purely  experimental.  Experiments  of  this  nature  are  costly  and 
cannot  be  carried  out  for  a  significant  number  of  damage  configurations. 
Consequently,  neither  of  the  studies  was  able  to  reach  a  general  conclusion. 
It  was  concluded,  "due  to  the  complexity  of  the  large  structure  configur¬ 
ation,  the  relationship  between  phase  changes  and  the  damage  location  still 
needs  more  research  effort"  [6]. 

The  report  will  cover  two  modal  theories  of  damage  diagnosis  developed 
by  the  author  and  his  co-workers,  using  the  change  of  dynamic  character¬ 
istics  of  the  damaged  structure.  The  modal  frequency  theory  uses  the  change 
of  modal  frequencies  to  establish  damage  functions,  from  which  the  intensity 
and  location  of  individual  crack  damage  are  diagnosed  [7,8].  The  appli¬ 
cation  of  the  theory  to  large  structures  and  to  multiple  crack  problems  -was 
made  possible  with  the  development  of  a  general  theory  of  circuit  analogy 
[9,10,11].  For  complex  structures  or  crack  damage  configurations,  there  are 
intrinsic  uncertainties  of  the  modal  frequency  theory  to  be  included  in  the 
present  report,  portion  of  which  was  presented  by  Ju  [12].  The  transmissi- 
bility  theory  uses  the  change  of  the  transmissibility  between  a  known  exci¬ 
tation  station  and  a  number  of  response  stations  to  diagnose  the  fracture 
damage.  Details  of  the  theory  and  application  are  included  in  this  report. 
The  feasibility  and  application  of  the  theory  were  earlier  presented  by 
Akgun  and  Ju  [13,14]. 

It  is  well-known  that  a  structural  member  deminlshes  its  load-carrying 
capacity  when  a  crack  is  developed  in  the  member  [15,16,17].  It  was  shown 
in  the  case  history  of  a  large  structure  that  the  frequencies  of  the  struc¬ 
ture  decreased  after  an  earthquake  [18].  Numerous  studies  have  been 
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reported  in  the  literature  on  the  effect  of  notches  and  cracks  on  beam 
behavior.  Static  deflection  of  beams  with  abruptly  changing  cross  sections 
was  studied  by  analyzing  an  equivalent  uniform  beam  with  modified  loading 

[19] .  This  approach  was  extended  to  vibration  of  bars  with  narrow  slots 

[20] .  An  experimental  and  analytical  study  was  conducted  earlier  on  narrow 

cracks,  in  which  bending  frequencies  of  a  beam  with  a  1.6  rrm  wide  crack  at 
the  center  were  measured  [21].  It  was  found,  for  this  particular  configur¬ 
ation,  that  analytical  and  experimental  natural  frequencies  agreed  when  an 
equivalent  slot  width  of  five  times  the  actual  width  was  used.  Such  a 
relationship  cannot  be  generalized,  however,  since  it  omits  the  effect  of 
slot  depth  [20].  A  combination  of  the  finite  element  and  transfer  matrix 
methods  was  used  recently  to  determine  the  effect  of  grooves  on  the  natural 
frequencies  of  beams  and,  as  an  example,  changes  in  the  fundamental  fre¬ 
quency  of  a  free-free  beam  with  a  slot  at  mid-span  were  determined  as  a  . 
function  of  the  slot  depth  and  width  [22].  The  method  of  representing  a 
slot  with  modified  loading  [19,20]  was  extended  to  a  study  of  the  effect  of 

a  crack  on  the  stresses  and  deflections  in  cylindrical  shells  and  beams 
under  dynamic  loading  [23,24].  In  the  model,  the  cracked  structure  was 
represented  by  an  uncracked  structure  with  the  local  reduction  in  stiffness 
due  to  crack  being  accounted  for  by  a  pair  of  concentrated  couples  M  located 
a  small  distance  c/2  on  either  side  of  the  physical  crack  location.  The 
method  requires  the  knowledge  of  static  radial  deflections  at  two  points  in 
the  case  of  a  circular  cylindrical  shell  from  which  the  product  Me  is  deter¬ 
mined  and  the  model  is  calibrated  [23].  The  static  modes  of  deformation  are 
then  used  to  compute  the  dynamic  response  of  the  structure.  The  required 
static  deflections  of  the  cracked  structure  need  to  be  measured  experi¬ 
mentally  or  computed  numerically  to  establish  the  dependence  of  Me  on  the 
crack  depth  for  the  particular  structure. 
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Chondros  et  al.  [25]  investigated  the  effect  of  a  crack  in  a  welded 
joint  on  the  dynamic  flexural  behavior  of  beams.  Their  model  was  a  massless 
cantilever  beam  with  a  tip  mass.  The  crack  at  the  welded  joint  was  modeled 
as  a  torsional  spring.  The  spring  constant  was  experimentally  determined 
from  the  changes  in  the  natural  frequency  for  different  crack  depths. 
Gudmundson  [26]  using  a  first  order  perturbation  method  for  small  cracks, 
showed  that  the  frequency  changes  are  functions  of  the  static  strain  energy 
due  to  the  crack.  Gudnundson  like  Ju  et  al .  [7,8]  related  the  strain  energy 
to  the  crack  depth  via  the  stress  intensity  factor.  In  the  case  of  a 
cantilever  beam,  his  analytical  results  showed  excellent  agreement  with  the 
results  of  an  experimental  study  except  when  the  crack  was  close  to  the 
built-in  end.  He  explained  the  disagreement  near  the  built-in  end  by  the 
influence  of  the  built-in  boundary  on  the  static  stress  intensity  factor. 
Gudmundson  [27]  recently  modeled  a  crack  by  a  static  flexibility  matrix 
which  accounts  for  the  dicontinuity  in  the  generalized  displacements  at  the 
crack  location.  He  obtained  the  flexibility  matrix  for  a  cantilever  beam  by 
using  static  stress  intensity  factors.  For  the  case  of  torsional  vibra¬ 
tions,  he  used  the  finite  element  mehtod  to  obtain  the  flexibility  matrix. 
Gudmundson's  experiments  with  a  cantilever  beam  confirmed  his  analytical 
results  very  well . 

Ju  [7,8]  proposed  the  theory  that  the  stiffness-softening  effect  of 
crack  damage  in  structures  can  be  effectively  represented  quantitatively  by 
a  spring-loaded  hinge  at  the  cracked  section  of  the  structure.  He  further 
postulated  in  his  theory  that  the  spring  constant  is  completely  defined  by 
the  configuration  of  the  crack  and  is  Independent  of  its  location  or  of  the 
frequency  of  osicllation.  The  determination  of  the  effective  spring  con¬ 
stant,  hence,  can  be  made  experimentally  or  analytically.  One  of  the  ana¬ 
lytical  techniques  for  the  determination  of  the  spring  constant  is  by 
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equating  the  strain  energy  stored  in  the  spring  with  the  surface  energy 
gained  through  the  creation  of  the  crack  surface.  The  fracture-hinge  theory 
has  recently  been  verified  ’by  experiment  to  be  reported  later.  A  portion  of 
the  work  will  be  presented  by  Ju  and  Mimovich  [28]  at  the  International 
Modal  Analysis  Conference  and  published  in  its  proceedings. 

The  numerical  analysis,  used  in  the  present  report,  is  based  on  the 
general  theory  of  circuit  analogy,  which  was  initially  develcpoed  in  Russia 
[29]  and  was  later  generalized  for  structural  dynamics  by  Ju  and  AkgUn 
[9,10,11].  The  fundamental  principle  lies  in  that  the  beam  equations  with 
the  assumption  of  modal  shapes  become  algebraic  and  are  mathematically 
equivalent  to  the  Kirchhoff's  circuit  theory  equations.  Oetails  were 
reported  in  [11].  The  work  has  been  generalized  to  include  damping  and 
extended  to  beams  of  non-uniform  properties  (geometrical  or  material).  In 
this  revised  form,  the  general  theory  of  circuit  analogy  is  included. 
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CHAPTER  2 
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GENERAL  THEORY  OF  CIRCUIT  ANALOGY 


The  structures  considered  in  this  study  are  beams  and  plane  frames  whose 
motions  are  described  by  the  Bemoulli-Euier  equation  with  damping.  The  method 
used  to  analyze  the  structures  is  the  electrical  analogy  method.  The  present  chapter 
is  devoted  to  the  development  of  the  analog  circuits  and  circuit  equations  for  a  beam 
and  for  a  general  frame  with  and  without  cracks.  First,  the  complex  Young's  modulus 
and  the  types  of  damping  which  can  be  incorporated  in  the  structure  model  are 
studied  in  the  next  section. 

2. 1.  Types  of  Damping  Encompassed  by  the  Model 

2.1.1.  Internal  Damping  -  Complex  Young’s  Modulus 

Metals  at  low  stress  display  linear  viscoelastic  behavior  [  30]  The  general  stress- 
strain  relation  for  a  linear  viscoelastic  material  is  given  by  [  31  ]  . 

(c0 +aIi?l+  •  •  •  va* Dn*...) a  =  (&„+•& (2.1) 

where  a*  and  bk  are  constants  and  denotes  the  L'th  partial  derivative  with  respect 
to  time.  When  c  and  c  vary  sinusoidally  with  time,  D*  -  (it;)*,  and 

~=  E(d)  -  Z*(y)  [l  +  i  d(y)]  (2.2) 

where  £{a)  is  the  complex  Young’s  modulus,  is  the  dynamic  Yeung’s  modulus 

and  6(a)  is  the  Loss  or  damping  factor.  The  real  part  of  £(u)  is  termed  the  storage 
modulus  and  the  imaginary  part  is  the  loss  modulus.  For  low  damping  materials  (5  of 
the  order  0.1  or  less),  variation  of  £  with  frequency  is  slow.  For  many  materials  of 
engineering  interest,  E*  and  o  may  be  treated  as  constants.  Common  metals  such  as 
aluminum  and  magnesium,  for  which  5  is  very  small,  fall  into  this  category  [  31  , 
pp.27,13?].  A  type  of  damping  which  is  sometimes  used  in  Sexural  and  longitudinal 
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vibrations  is  the  internal  strain  velocity  damping  (viscous)  for  which  the  constitutive 
equation  is  ■  . 


c  =  E*  t  +bDi  s  (2.3) 

where  the  second  term  on  the  right  is  the  damping  stress  [32]  For  harmonic 
vibrations, 


iboi 


(2-4) 


which  yields  a  damping  factor  proportional  to  frequency.  Experimental  evidence 
contradicts  such  dependence  of  damping  on  frequency.  A  more  realistic  damping 
factor  may  result  if  higher  order  time  derivatives  are  used  in  Eq.  (2.3).  However,  the 
damping  force  in  real  structures  appears  to  be  in  phase  with  the  velocity  but 
proportional  to  the  displacement  [  32p.7S].  Such  damping  is  termed  structural  or 
hysteretic  damping.  The  constitutive  equation  can  then  be  written  as 

a  -  Ei  £  +  ids.  (2.5) 


The  complex  modulus  for  structural  damping  then  follows  from  Eq.  (2.5);  that  is. 

E  -  Ed  (1  +  Jr-)  (2.5) 

where  -or— -is  the  structural  damping  factor  which  is  independent  of  frecuency.  In  this 

hi 

study,  E  will  be  used  In  the  form  given  by  eq.  (2.2)  without  referring  to  the  type  of 
damping. 

The  defection  of  a  beam  under  transverse  vibration,  excited  by  a  sinusoidal 
force,  is  governed  by  [  31  p.  196] 


,d? 


eiut- 


0  ' 


(2.7) 


where  £  =  —is  the  normalised  axial  coordinate,  y  =  y(£)  is  the  complex  modal  shape 
and 


2  f  ^ 

where  °  and  5*.  at  resonance,  is  the  und&moed  characteristic  value.  /?, 


near  resonance,  will  be  called  the  "complex  characteristic  value".  "When  5  is  small, 

(2.9) 


jff  «  £d(l  -i 


In  the  study,  however,  exact  values  of  £  are  used. 

At  this  point,  a  clarification  on  terminology  will  be  made.  It  is  appropriate  to  call 
y  the  -complex  modal  shape,  and  this  term  has  been  used  in  the  literature.  The 
physical  shape  of  a  damped  beam  under  forced  vibration,  however,  cannot,  in  general, 
be  referred  to  as  the  modal  shape.  This  is  a  consequence  of  the  fact  that  the  physical 
beam  deflection  cannot,  in  general,  be  expressed  as  the  product  of  a  time  function 
and  a  space  function.  The  argument  can  be  illustrated  by  taking  the  real  part  of  the 
complex  deflection.  Namely, 

Wr  =  Re  [y  siui  ]  =  ym(£)  cos  [o  t  -  «($)]  (2.10) 

where  uv  is  the  physical  deflection,  and  ym  and  p  are  the  magnitude  and  phase  angle 
of  the  complex  modal  shape  y,  both  of  which  are  functions  of  £  and  u.  Only  when 
p  -  0  (Le.,  no  damping),  can  ym  be  referred  to  as  the  modal  shape.  In  this  study,  ym 
is  termed  the  amplitude  of  deflection. 


2.1.2.  External  Viscous  Damping 

■When  there  is  external  viscous  resistance  to  the  transverse  harmonic  vibration  of 
a  uniform  beam,  the  governing  equation  is  again  of  the  form  of  Equation  (2.7)  where  £ 
is  now  deflned  as 

plil  (2.U) 


where  c  is  the  coefficient  of  viscous  damping.  Therefore,  the  methods  cf  this  study 
are  also  valid  when  the  damping  is  of  the  external  viscous  type  or  a  combination  cf 
external  and  internal  types. 
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In  the  remaining  part  of  the  chapter,  the  basic  circuit  analogs  of  a  beam  and  a 
general  frame  with  cracks  are  derived  based  on  Eq.  (2.7).  In  the  frequency  theory 
(Chapter  3),  the  structures  are  assumed  undamped  and  therefore  p  is  real  In  the 
transmissibiiity  theory  (Chapter  4),  damped  structures  are  considered,  for  which  p 
is  complex  The  same  analog  circuits,  however, '  are  used  in  both  methods,  with  p 
chosen  properly. 

2.2.  Basic  Analog  Circuits 

Keropyan  and  Chegoiin  gave  the  T-  and  Il-circuit  analogs  cf  vibrating  beams  [29] . 
The  analog  circuits  will  be  rederived  here.  The  complex  mode  shapes  of  a  Bemoulli- 
Euler  beam  element  under  transverse  vibration  are  given  by  the  solution  to  Equation 
(2.7).  That  is, 

y(q)  =  AcashPf  4  5  sinh  p£  4  Ceos  p$- 4  D  sin  p$  (2.12) 

Four  variables  are  associated  with  each  end  of  the  beam  element,  namely,  deflection, 
y,  slope  (or  angle  of  rotation),  y\  resisting  moment,  M,  and  shear  force,  V  (Figure  1). 
From  Equation  (2.12)  and  its  appropriate  derivatives  at  £  =  0,  the  variables  at  the  left 
end  are  obtained: 

y  x  ~  A  -r  C  yl'  =  (5  4i?)£- 

Jml 

■  =  (A  -  C) Vl*(E-3)&£-  (2.13) 


From  the  equations  in  (2.13),  the  four  coeScients  can  be  solved  fer  in  terms  cf  the 
four  variables.  Namely, 


-  _  (y,  4  h  Vi )L 

2  .2 p 


*  _  (yi-kJii)  tt  _  (yi’  4  h  vql 

c  ~  2  ■  ~  2  P 


where  h  — 


L a 


EIp* 


(2.14) 

The  variables  at  the  right-end  of  the  beam  element  are  obtained 


from  Equations  (2. 12)  and  (2. 14)  with  £  =  1.  The  result  is 


n 


yt 

yi 

Mi 
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Figure  1.  Beam  element  under  free  vibration. 


Vz  =  Ay,  +  Gy ,'|-+  CM,k  -  DV,kj- 

3  BM:L 

yi  -  Py  i2"+  -*/if  + 

a  =  +  Dyi'Slj-  +  AJTX  -  ■ 

v*s-&ifc-cQm~DSllj-+Ayl  (2.15) 


-where 


.  _  (casks  +■  cos  3) 

„  _  (sinks  +  sin  3) 

2 

2 

~  _  ( cash  S  -  cas  S) 

n  (sinks  —  sin  3 ) 

2 

2 

(2.16) 


Any  four  of  the  eight  variables  In  Equation  (2.15)  can  be  soived  for  in  terms  of  the 
remaining  four  variables.  The  mathematical  analogy  between  electrical  circuits  and 
beams  is  based  on  the  liner  transform  (2.15).  The  type  and  properties  of  the  resulting 
circuit  depend  on  the  choice  of  the  independent  variables.  In  this  study,  angles  of 
rotation  are  analogous  to  voltages  and  resisting  moments  are  analogous  to  electrical 
currents. 


2.2.1.  T  Circuit 

When  the  condition  sin^^O  is  imposed,  (,5  is  real  for  an  undamped  beam  )  slopes 
and  shears  can  be  ezpressed  in  terms  of  resisting  moments  and  defections.  Thus, 
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froca  Equation  (2.15) 
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where  h  — 


S  =  cath  §  —  cat  § 
S’  =  cath  p  +  cat  (3 
T  —  esc  (3  -  csch  (3 
T  —  esc  jS  +  csch  (3 


The  first  two  equations  in  (2.17)  can  be  rewritten  as 


y  {  —  (Z-rZo)IrI\  —  Z^Uz  *  £*;  .  Vz  =  ZC.V 1  --  (Z  +•  Zz)J/z  +  Z2 

where 


Z  o  = 


7X 

2Z/,3 


Ei*(T'yi-S’y  :)£- 


(2.17) 


(2.13) 


(2.  IS) 


(2.20) 


(2.21) 


The  two  equations  in  (2.17'  are  the  XirchhcfTs  equations  for  the  active  three 
terminal  network  shewn  in  Figure  2  (F  circuit)  with  siope  and  moment  (y\  .«/)  being 
analogous  to  eiectrioai  voltage  and  current,  respectively.  Tne  quantities  Z  and  Z. 
denote  electrical  impedance  and  voltage  source,  respectively.  Negative  resistance 
(real  part  of  Z)  poses  no  difficulty  in  analytical  and  numerical  analyses.  The  ends  of  a 


beam  segment  are  simulated  by  the  ports  of  the  circuit.  It  must  be  observed  that  the 
electromechanical  analogy  described  above  does  not  simulate  the  differential 
equation  of  motion,  but  instead  the  solution  based  on  the  assumption  of  harmonic 
motion  [29]. 
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Figure  2.  T-circuit  analog  of  a  beam  under  free  vibration 

Boundary  conditions  at  the  ends  of  a  beam  element  can  be  simulated  as  follows: 

a.  Free  end  or  simply-supported  end,  .rJ=Q:  the  corresponding  pert  of  the  circuit 
is  left  open,  which  results  in  zero  electrical  current. 

b.  Fixed  end,  y'~ 0:  The  corresponding  port  is  short-circuited,  which  results  in 
zero  voltage.  The  other  boundary  conditions  at  these  ends  and  the  last  two  Equations 
in  (2.17)  are  utilized  to  express  the  voltage  sources,  Equation  (2.2l),  in  terms  of  the 
moments,  as  will  be  illustrated  with  an  undamped  cantilever  beam  vibrating  freely 
(Figure  3),  Ail  of  the  circuit  variables  are  real  in  this  case.  The  KirchofTs  voltage  law 
applied  to  the  Iccp  yields 

(Z  +  Zc)SIl  +  Ex  =  0  (2.22) 

where  Z  is  the  analog  resistance.  From  the  last  equation  in  (2.17)  with  y  i  =  Jfj  =  0, 

y2  =  0  =  (rtf:-S^^-  (2.23) 

from  which 

ar.d  '  (Z2i) 

The  second  equation  in  (2.24)  Is  obtained  from  Equation  (2.21). 
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Figure  3.  Cantilever  beam  and  Its  T-circuit  analog 


The  substitution  of  Equations  (2.20)  and  (2.24)  into  (2.22)  and  the  fact  that 
result  in 

5  =  0  (2.25) 

which  yields  the  characteristic  equation  of  a  cantilever  beam.  1  +•  coshjS  csc/3  =  0. 

f 

In  the  derivation  of  Equation  (2.17),  it  is  theoretically  sound  to  impose  sin/2  ^  0. 
However,  during  the  numerical  search  of  the  undamped  natural  frequencies,  sin/3 
may  get  very  small  at  some  point,  or  the  structure  may  actually  have  frequencies  at 
or  near  sin£  =  0.  The  soiution  to  the  problem  lies  in  realizing  that  /3  is  directly 
proportional  to  the  length  L  of  the  beam  element.  The  problem  can  thus  be  solved  by 
dividing  the  beam  Into  two  sections  and  representing  each  section  by  an  individual 
circuit.  To  illustrate  the  approach,  an  undamped  simply-supported  beam  is 

arbitrarily  divided  into  two  elements  with  a  length  ratio  of  --  (Figure  4).  With 

6 

reference  to  the  figure,  the  continuity  conditions  in  slope  and  moment  (voltage  and 
current)  at  the  common  boundary  is  preser/ed  under  the  cascade  connection  of  the 
two  basic  circuits.  The  loop  equation  then  yields 
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2.2.2.  n  Circuit 


Equation  (2.32)  Is  the  Kirchhoffs  current  equation  for  the  active  three  terminal 
network  shown  in  Figure  5  (II  circuit)  with  G  and  /  denoting  an  admittance  and  a 
current  source,  respectively.  Moment,  M,  and  slope,  y‘,  are  again  analogous  to 
current  and  voltage.  As  an  illustration,  the  cantilever  beam  of  Figure  3  will  now  'je 
simulated  with  a  II  circuit.  The  left  port  of  the  circuit  in  Figure  5  is  shorted  yielding 
the  circuit  of  Figure  S.  The  node  voltage  equation  is 
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2.3.  Simulation  of  Crack  with  Circuit  Analogy 

A  cracked  section  in  a  beam  is  modeled  following  the  method  of  fracture 
binge  [7,8] .  The  effect  of  a  crack  is  to  introduce  local  flexibility  in  the  beam.  The. 
slope  of  the  mode  shape  Is  discontinuous  at  the  cracked  section.  The  crack  can  thus 
be  mechanically  represented  by  a  torsional  spring  of  spring  constant  «.  The 

M 

discontinuity  in  the  slope  is  given  by  Ay'  =  — — where  Mc  is  the  resisting  moment  at 

the  cracked  section.  Such  a  discontinuity  is  analogous  to  a  voltage  drop  in  the  circuit 
theory,  when  slope  and  moment  are  analogous  to  voltage  and  current,  respectively. 
The  crack  can  then  be  simulated  with  a  resistor  of  resistance  —1/  k.  Hence,  a  beam 
with  a  single  crack  is  represented  by  two  circuits  (T  or  II)  joined  by -a  "crack 
resistance,"  —1/  K. 

The  derivation  cf.«  is  based  on  equating  the  energy  stored  in  the  torsional  spring 
to  the  increase  in  strain  energy  of  the  beam  due  to  the  crack.  For  a  plane  strain 
crack,  this  approach  yields  £7,81. 

(2.38) 

where  y  -  is  the  relative  crack  depth  (Figure  7)  and  /  is  the  dimensionless  stress 

intensity  factor  for  symmetric  cracks  given  by 

f(y)  -  *[1.122  -  2.3837  +  4.3S7-/2  -  4.5S-/3  4  2.3457*  -  0.36375]  (2.33) 

The  analysis  of  the  next  chapter  will  indicate  that  the  crack  intensity  is  quantified  by 
the  following  nondimensional  number; 


Equation  (2.38)  reveals  that  9  is  a  function  of  the  beam  slenderness  ratio,  2b/ L,  and 
the  relative  crack  depth,  y.  The  analysis  indicates  that  as  the  value  of  9  increases, 
the  reduction  in  the  natural  frequencies  of  a  beam  also  increases.  Hence,  of  the  two 
beams  with  a  crack  of  the  same  relative  depth,  at  the  same  relative  location,  the 


frequencies  of  the  more  slender  beam  Trill  experience  smaller  changes.  Therefore.  0 
is  called  the  "sensitivity  number". 

'  2.4.  Analog  Circuit  for  A  Plane  Frame 

Xeropyan  and-  Chegciin  presented  some  examples  of  the  application  of  eiectrical 
analogy  circuits  to  the  solution  of  static  and  dynamic  frame  problems  [29].  The 
approaches  presented  by  them  to  frames  vrith  freedom  of  lateral  motion  Involve  an 
iterative  solution  or  a  several  step  procedure.  In  the  former,  the  amount  of  sidesrray 
is  guessed  and  iterated  on.  In  the  latter,  the  frame  is  first  restricted  laterally,  the 
reactions  at  the  restrictions  are  solved  for,  end  the  sidssvrav  is  then  computed.  Both 
of  these  approaches  become  complicated  and  uneconomical  for  multi-story  frames. 


Keropyan  and  Chegolin  illustrated  the  latter  approach  with  a  static  one-story  frame 
problem-  Their  examples  of  dynamic  frames  involved  laterally  constrained  frames. 

This  section  expounds  a  generic  model  of  a  two-dimensional  multi-story  frame 
(Figure  S)  with  n  stories  and  m  spans  (thus  m+1  anchors).  Damage  diagnosis  on  such 
a  frame  using  both  the  frequency  and  the  transmissibility  methods  will  be  developed 
in  later  chapters.  In  the  study  of  the  frequency  method,  free  vibration  of  the  frame 
with  and  without  a  crack  is  analyzed  from  which  the  characteristic  equation  is 
derived.  In  the  development  of  the  transmissibility  method,  transverse  defection  at  a 
certain  location  is  computed  under  a  sinusoidal  excitation  force.  The  analog  circuits 
are  the  same  in  both  methods  with  the  exception  that,  in  the  second  method,  0  is 
complex  and  some  additions  are  made  to  the  circuit.  In  this  section  the  circuit 
equations  are  derived  in  their  most  general  form  to  be  applicable  to  both  methods.  A 
formalized  scheme  Is  developed  for  obtaining  the  circuit  equations  of  a  general  frame 
with  or  without  cracks,  from  which  the  modal  frequencies  or  desired  kinetic  and 
-kinematical  quantities  can  be  computed.  In  the  7*-circuit  analogy,  the  circuit 
equations  for  the  structure  are  developed  through  the  mesh  current  (moment) 
equations.  Hence,  the  unknowns  are  the  moments  at  the  frame  joints  and  at  the 
cracks.  There  are  p  -1  unknown  moments  at  a  frame  node  where  p  beams  join 
together.  Because  the  moment  is  continuous  across  a  crack,  there  is  one  unknown 
moment  for  each  crack.  It  follows  that  the  order  of  the  resulting  linear  system,  which 
is  equal  to  the  number  of  unknowns,  is  [n(3m.  +  l)-r«  +2]  for  a  frame  with  cracks  and 
one  excitation,  one  of  the  unknowns  being  the  value  of  the  kinetic  or  kinematical 
quantity  at  the  desired  location  on  the  structure.  Cn  ths  other  hand,  in  the  H-circuit 
analogy,  node  voltage  (siope)  equations  lead  to  the  circuit  equations.  Slope  is 
continuous  at  a  frame  joint;  therefore,  there  is  one  unknown  slope  for  each  frame 
joint  regardless  of  the  number  of  beams  connected  there.  However,  there  are  two 
unknown  slopes  at  each  crack  location,  one  on  each  side  of  the  crack.  In  other  words, 
slope  is  discontinuous  across  a  crack.  The  structure  is  assumed  rigidly  fixed  at  the 
ground  level;  hence,  ths  slopes  are  zero  at  the  anchors.  The  order  of  the  system  is 
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t 


thus  [n(m  +  l)+2.fe+2]  •with  II  circuits.  The  H-circuit  analogy  is  therefore  preferred 
over  the  T -circuit  analogy,  since  the  order  of  the  system  for  k  <  2 rrsi  is  smaller  in 
the  former. 

The  procedure  to  obtain  the  circuit  equations  of  a  frame  structure  are  outlined 
next. 

(l)  II  circuits  simulating  individual  columns  and  girders  joining  at  right  angles  are 
interconnected  such  that  the  boundary  conditions  at  the  analog-frame  joints  are 
satisfied.  (It  will  be  shown  that  a  network  diagram  need  not  actually  be  drawn.) 


(2)  For  each  crack  present  on  any  element,  one  more  II  circuit  and  one  crack 
resistance  (hence,  two  mere  unknown  slopes-voitages)  are  added  to  the  network. 
A  girder  or  a  column  with  p  cracks  is  thus  simulated  with  p  -1  circuits 
connected  in  cascade  via  the  crack  resistance.  The  first  and  the  (p  +  l)t/i 
circuits  are  then  connected  with  the  rest  of  the  network  according  to  Step  1. 

(3)  The  beam  element  with  the  excitation  is  divided  into  two  elements  at  the  point  of 
application  of  the  excitation  and  is  thus  represented  by  two  IT  circuits  instead  of 
one. 

(4)  When  a  kinetic  or  kinematical  quantity  is.  desired  at  some  location  on  the  frame, 
the  corresponding  beam  element  is  again  divided  into  two  beam  elements  at  the 
desired  location  and  represented  by  two  II  circuits.  Steps  3  and  4  are  skipped  in 
the  frequency  method. 

(5)  XLrchhofiTs  current  Law  is  written  in  terms  of  the  node  voltages  which  correspond 
to  the  angles  of  rotation  at  the  frame-joints  and  at  the  cracks. 

(6)  Kinetic  conditions  at  each  floor  level  and  shear  conditions  at  the  cracks,  at  the 
excitation  location,  and  at  the  location  where  deflection  or  any  other  quantity  is 
desired  are  applied.  This  allows  the  lateral  motion  of  the  floors  to  be  computed 
directly. 

(7)  This  step  is  performed  when  the  frequency  method  of  diagnosis  is  used.  Under 
free  vibration,  the  coefficient  matrix  obtained  in  Step  3  must  be  singular.  If  the 
damage  parameters  (Le.,  location  and  severity)  for  each  crack  are  known,  the 
characteristic  ’values  of  the  structure,  j,5j,  can  be  determined  from  the  zero 
determinant  of  the  coefficient  matrix  (the  forward  problem).  If  the 
characteristic  values  are  known,  the  damage  parameters  can  be  solved  by  using 
the  characteristic  equation  (the  inverse  problem). 

(9)  When  the  transmissibility  method  of  diagnosis  is  used,  the  linear  system  of 
equations  is  solved  for  the  deflections  in  terms  of  the  excitation  force. 
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The  approach  taken  here  results  in  a  formal  scheme  ■which  can  be  applied  to 
frame  structures  without  referring  to  an  actual  network  diagram.  Reference  is  made 
to  Figure  S.  The  word  "walT  will  denote  the  union  of  all  the  columns  on  the  same 
vertical  Line.  Each  node  and  each  beam  element  is  identified  by  a  number.  Numbers 
corresponding  to  the  nodes  are  encircled  in  the  figure.  The  numbering  order  for  the 
beam  elements  begins  at  the  left  lowest  column,  proceeds  up  through  the  columns  on 
the  first  wall,  returns  back  to  the  second  anchor,  proceeds  up  vertically,  and 
continues  in  that  order.  Once  the  columns  are  finished,  girders  on  each  floor  are 
numbered  progressively  from  left  to  right  starting  with  the  first  floor  and  continuing 
on  with  the  upper  floors.  Quantities  related  to  a  beam  element  such  as  G,  P,  Q,  etc. 
are  subscripted  by  the  number  of  the  element.  Each  node  is  numbered  by  the 
column  under  the  node.  If  there  are  v  cracks  on  a  column  (girder),  the  original 
number  of  the  beam  element  refers  to  the  uppermost  (rightmost)  segment  cf  the 
column  (girder)  which  is  now  represented  by  v  +  1  beam  segments.  The  other 
segments  are  numbered  following  the  largest*  number  in  the  scheme.  Due  to 
discontinuity  cf  slope  at  a  crack,  two  new  nodes  are  created  on  the  two  sides  cf  each 
crack.  The  location  t  on  the  structure,  henceforth  referred  to  as  the  response 
station,  denotes  the  location  at  which  deflection  or  any  other  quantity  is  to  be 
computed.  One  more  node  is  thus  created  and  one  more  II  circuit  is  added.  Another 
node  is  added  on  the  column  2n+2  where  the  force  Is  applied.  These  last  two  nodas 
and  the  beam  segments  below  (to  the  left  of)  them  are  numbered  the  last  in  the 
numbering  scheme.  • 

2.4. 1.  Analog  Circuit  Equations  and  Boundary  Conditions 
a)  Nodal  Ecuatior.s: 

Figure  9a  shows  a  typical  node  r  at  the  intersection'  of  column  elements  r,  s  end 
girder  elements  p,  q.  In  reference  to  Figures  9  and  9a,  the  sum  of  the  branch 
currents  entering  the  node  r  is  set  equal  to  zero.  Namely, 
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(yy.  -Vt')^  (y/  ~yr')G7c  +  (yjt’  -yr')$-o +  (vt‘  -yrO^o 

~yr'{Gp  £7  +  £■  +  Cra)  +  TpZ  -  I<!\  *  IrZ  -  h  \  ~  Q  (2.41) 

« 

■where  (yi’  —  yr’)  is  the  analog  voltage  difference  between  the  nodes  i  and  r  and  £,0  is 
the  conductance  which  connects  these  two  nodes.  The  current  sources  appearing  in 
equation  (2.41)  are  dependent  on  the  transverse  deflections  of  the  two  ends  of  the 
beam  elements  to  which  they  correspond.  Under  small  deformation  theory,  vertical 
deflection  of  the  nodes  shown  in  Figure  S  is,  by  second  order  approximation,  equal  to 
zero.  Thus,  in  Figure  9a.  y?2  -  y7i  =  0.  where  yp2  and  y?1  denote  the  deflections  of 
the  right  end  of  the  pth  girder  and  Left  end  of  the  gth  girder,  respectively.  If,  in 
addition,  the  nodes  t  and  j  do  not  bound  cracks,  then  y?i  yq2  -  0.  Consequently, 
I-oi  -  \  ~  0  from  Equation  (2.34).  On  the  other  and,  /r2,  for  instance,  is 

not  zero  since,  in  general.  yri  *  0  and  yr2  =  y,i  -A  0  (i.e.,  sidesway  is  allowed). 

Upon  substitution  of  equation  (2.33)  into  (2.41), 

-Vr'ilplQ?  +VqlQl  +~rlGr  +  -Vi'Vpl^p  -y/Vql^ 

-yk'iri^r  -yi'VsiWi  +  {Ipz  -  Art  Itz  “At)  =  o  (2.42) 

jo 


0 


(*)  ®  Cb) 


Figure  3.  Typical  Nodes  on  the  Frame 
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where  rjt  t  = 


(■ EIfi/L)i 

(Zip/ Do 


and 


El  6 


L 


is  the  characteristic  cr  the  reference  value 


chosen  for  the  frame. 

In  the  formation  of  crack  nodes.  Figure  3b  snows  two  nodes,  r  and  r',  across  a 
crack  on  a  beam  element.  Xirchhoffs  current  law  is  written  for  the  node  r: 

(yi  -yr')OrQ  +  {y/  -Vr ')(-«)  -Vr'Gr  +  Ui  =  0  (2.43) 


where  yr‘  is  the  rotation  of  the  node  r1  and  the  conductance  -«  represents  the  crack. 
The  nodal  equation  for  the  node  r*  is  of  a  similar  form.  Equation  (2.43)  is  valid 
whether  the  crack  is  on  a  girder  or  a  column.  The  same  steps  which  led  to  equation 
(2.42)  yield 


“Vr' 


FI  r\Qr 


3fr(l^o)J 


Vr 


3f;(  1  -rl  <5  ) 


/r2  =  0  (2.44) 


The  node  f  (Le.,  the  response  station)  or  the  node  under  the  force  in  Figure  3  is 
the  same  as  that  in  Figure  3a  when  the  branches  p  and  c  are  removed.  The  ncdai 
equation  for  such  a  node  is  then  obtained  from  Equation  (2.41)  cr  (2.42)  by  setting  all 
of  the  variables  with  a  subscript  p  and  q  equal  to  zero. 


b)  gour.d3rv_.and  Continuity  Conditions: 

From  the  small  deformation  theory,  points  on  girders  on  the  same  doer  level  will 
have  the  same  horizontal  displacement.  Namely, 


Vij  —  Vn  *-iJ  =  l/ti-KJ:  =  '  =  l/mn+ij  i  =  1,....7»  j  —  1.2  (2.45) 

Also,  y»2  =  Vu-IA  .  i  ~  “1.  On  the  other  hand,  defection  continuity  across  a 

crack  implies  that  defections  of  the  nodes  r  and  r1  in  Figure  3b  are  equal  Defection 
Is  continuous  at  the  excitation  node  ai  d  at  the  node  f  also  (Figure  9).  Thus,  there  are 
n+k+2  unknown  defections  associated  with  the  frame,  k  denoting  the  total  number 
of  cracks  on  the  frame.  If  the  horizontal  displacement  of  the  ith  floor  is  designated 

by  Vt  (Vi  s  Viz .  i  =  l n),  an  [n+i+2]-vector  of  displacements  ^and  an 

[n(m.+i)-i-2.fc+2]  -  vector  of  rotations  are  obtained  as 
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7  =  Vi  yz  •  •  •  Vn  Vei  •  -  •  Vc*  Vf  yt 


T  s  ■  vi  •  •  ■  y‘n(m*D  Ve,'  &x*  •  •  •  y/  yt' 


(2.45) 


•where  yCl  designates  the  transverse  deflection  at  the  ith  crack,  and  yCi'  and  yei‘ 
designate  the  rotations  immediately  to  the  left  of  (belcrw)  and  to  the  right  of  (above) 
the  ith  crack  "which  is  on  a  girder  (column).  The  subscripts  /  and  t  denote  the  node 
under  the  force  and  the  node  i,  respectively.  The  unknown  deflections  can  be  related 
to  the  rotations  via  the  kinetic  equations,  defining  the  sideway  motion  of  the  floors, 
and-  the  shear  continuity  conditions  at  the  cracks.  Cross  shears  at  the  column  ends 
are  the  axial  driving  forces  for  the  floors.  The  axial  acceleration  of  a  floor!  which  is 
composed  of  all  the  girders  on  the  same  level  is  the  same  as  the  transverse 
acceleration  of  the  column  ends  at  the  nodes  where  the  columns  join  the  floor.  Thus, 
for  the  ith  floor  (Figure  10) 


TO  P+-TO-I  ■ 

E['(T>n.J  -  =  y,  Pi  Li*  E 

}-C  ;=?  J=? 


(2.47) 


where  p=(m-l)r.4-(i— l)rn.41  in  accordance  with  the  numbering  order  described 
earlier  and  k.  is  the  total  number  of  cracks  on  the  ith  floor.  There  is  one  such 
equation  for  each  floor.  The  term  in  the  braces  on  the  right-hand  side  of  Equation 
(2.47)  is  the  total  mass  of  the  ith  floor.  When  there  is  no  crack  on  the  ith  floor,  Le., 
2.  =  0,  the  second  summation  term  on  the  right-hand  side  is  absent.  When  there  is  a 
crack  in,  for  example,  girder  p-i-1,  £p  +  j  denotes  the  length  of  the  right  segment  of 
this  girder.  Mass  P74j  of  the  left  segment  is  then  Included  in  the  second  summation 
term  in  Equation  (2.47).  Equation  (2.47),  upon  substitution  of  the  iast  two  equations 
in  (2.30)  for  shears,  will  relate  deflections  to  rotations.  Figure  11b  shews  the  diagram 
of  a  frame  with  rn  =  1,  n  =  2,  fc=  2.  A  harmonic  load  is  applied  to  the  frame  and  a 
response  station,  t,  is  shewn.  There  are  tea  nodes,  ten  beam  elements  (hence  ten 
interconnected  II  circuits)  and  six  unknown  deflections.  The  first  floor  has  one  crack 


? 


on  it.  The  kinetic  equation  for  the  first  fiocr  with  pe  -  P5.  has  the  form 


V12+  V«2-  721  —  7?1  =  *  la) 

where,  for  instance, 


^32  s  ~  Pg/lPz/  ^2) 


Q?'yz'  +  ^'y?' 


(■P7I/1  ~  PiVi)fa 

L 7 


(2.48) 


(2.49) 


from  Equation  (2.30)  with  ySJ  »  ysi'  =  0,  y>  =1/32=  y?:,  Vs'  =■  722'  =  y?:\  y7  =  y?* 
y?'  sty-2’.  Upon  substitution  of  shear  expressions,  Equation  (2.49), 'and  the  similar 
equations  for  7l2  and  721  into  Equation  (2.48)  the  following  equation  is  obtained 
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|  foia^i  +  V22P2  +  77.3^3  +  ijnP7  +  Vsz(?a  +  £s)]Vi 

‘"’VzaPz'yz  ~  772-^?'y?  j  ~  (VizQi  ~~  VzzQz)y  i’ 
—  7722^2*1/2’  +  (vzzQz  ~“r)7zQi)y i  ““Vn^iyi 


(2.50) 


•where 


Vi]  = 


(EI)i 

ft  V 
A 

(El)  0 

>0  )J 

ia  j 

(2.51) 


In  Equation  (2.51).  i  is  the  number  of  the  beam  element  or  the  beam  segment  and  j  is 
an  exponent  The  kinetic  equation  for  the  second  floor  can  be  developed  similarly. 
Application  of  shear  continuity  at  a  crack  will  now  be  demonstrated  for  the  crack  on 
the  flrst  floor  (Figure  lib).  Hence  V&  =  7S;,  cr  from  the  last  two  of  Equation  (2.30) 

^  (Wa'yi  +  Ge'Va'  -  ^— )  =  ($s’SV  +  ^’Vs'  +  )  (2.52) 

since  t/31  =  7/32  =  0.  The  following  expression  is  obtained  by  rearranging  Equation 

^3  ft a 

(2.52)  and  noting  that  h.  =  A,a  and  7—  =  7 — : 

I13  L  a 

f.^3  •  n  —  or  .  1  si/  ...  1  •  >«.  1  /1  t  1  /<-> 


7—  (Pa  +  P a)Va  =  J*Vyi’  -  '’Vya'  *  <?a'Va’  -  Qa'V s' 
pa 


(2.53) 


Shear  continuity  at  the  response  station  yields  an  equation  similar  to  (2.53).  The 
shear  discontinuity  at  the  excitation  location  on  the  other  hand,  can  be  expressed  as 

=  Vn  +  (2.54) 


which  yields 


T~  C( Pi  +  ^a)y/  -  =  ($V  -  <?i')V/'  -  AW  +  (2.55) 
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For  the  general  frame  structure,  there  are  n  kinetic  equations  of  the  form  cf 
equation  (2,50),  /fc  +  1  shear-continuity  equations  (k  equations  for  the  cracks  and  one 
for  the  response  station)  of  the  form  of  equation  .(2-53),  and  one  shear  discontinuity, 
equation  given  by  (2.55).  These  n±k  +2  equations  can  be  arranged  in  a  matrix  form  to 
solve  for  the  deflections  in  terms  of  the  rotations,  namely. 


(2.58) 


where  e  is  a  column  vector  whose  entries  are  zero  except  for  the  last  entry  which  is 
unity  and  hj  is  the  value  pertaining  to  the  beam  element  with  the  excitation  (for  the 
example  of  Figure  11.  h ^  -  Ax).  On  the  other  hand,  the  nodal  equations  (2.42)  and 
(2.44)  can  be  arranged  in  the  form 


Hji-Uy' 


*0 


(2.57) 


where  1  is  the  vector  of  current  sources  and  is  related  to  the  deflections  via  Equation 
(2.34);  that  is, 

(2.58) 

After  Equation  (2.5S)  is  substituted  into  Equation  (2.59),  then  Equation  (2.52)  into 
Equation  (2.57),  and  it  is  noted  that  HiHg  =  X* ,  the  slope  equation  results: 

(Xf  2TlX-  U)y '  =  h.fFmXr7Tle  (2.59) 

This  is  the  general  form  of  the  result  for  any  planar  .frame  with  any  number  of  cracks 
on  it.  The  coefficient  matrix  in  (2.59)  is  symmetric.  The  individual  matrices,  in 
general,  have  the  forms 
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n(m+l)  2Jt,  Zfcj  2 


(2.80) 


n(m+l)  2Jb«  2^  2 


n  kc  k,  2 
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Ccr  |  Dc  |  0 
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Nf 
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Ns 


(2.81) 


■with  y  and  y '  arranged  as  in  equation  (2.4S) 


* 

r  i 

7c 

7r 

7« 

7s ' 

7, 

7'  = 

7/ 

yj 

V/’ 

Vt 

i  . 

< 

(2.52) 


where  the  subscripts  c  and  g  refer  to  cracks  on  the  columns  and  on  the  girders, 
respectively,  with  ks  and  k3  being  the  total  number  of  cracks  on  the  columns  and  on 
the  girders  ( k  =  kc  4-  k3y  Hence,  ye  and  y3,  for  example,  are  the  vectors  of 
deflections  at  the  column  and  girder  cracks,  respectively,  y 5  is  the  vector  of  the 
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horizontal  floor  displacements.  y.  and.  y}  are  the  transverse  deflections  at  the 
response  station  and  at  the  excitation  Location,  respectively.  It  should  be  noted  that 
xf",  u„.  and  2c  for  the  case  of  no  crack  are  modified  when  cracks  are  introduced  on 
the  frame.  The  matrices  U  and  Z  are  symmetric. 

"When  the  rotations  are  to  be  computed,  j'  can  be  solved  for  from  Equation 
(2.59).  When  the  deflections  are  needed,  y’  can  be  solved  for  from  (2.57)  and 
substituted  into  (2.56)  from  which  y  cam  be  obtained.  This  yields 


y  =  hfFm±  (xtrY-zre 

P  a 


(2.63) 


With  y  and  y'  known,  resisting  moments  and  shears  can  be  computed  from  Equation 
(2.30).  When  the  free  vibration  of  the  structure  is  of  interest,  Fm  ~  o  and  the 
coefficient  matrix  on  the  left-hand  side  of  (2.59)  is  singular.  The  determinant  of  the 
coefficient  matrix  then  yields  the  characteristic  equation  of  the  structure.  A 
computer  program  is  written  to  compute  the  deflections  on  a  general  frame 
structure.  The  listing  of  the  program  is  given  in  Appendix. 

A  formal  procedure  will  now  be  described  to  establish  the  above  matrices  without 
the  need  to  draw  the  actual  analog  circuit. 

2.5.  Procedure  to  Establish  the  Matrices 

In  the  following,  for  simplicity  of  explanation,  properties  of  the  beam  elements 
are  assumed  uniform  throughout  the  frame  except  for  the  lengths.  That  is,  rt :*,■  =  1  for 
all  i. ; . 


U  is  composed  of  the  coefficients  of  the  rotations  yf'  in  equations  such  as  (2.42) 
and  (2.44).  The  ith  row  in  the  matrix  stems  from  KLrchhoffs  current  law  written  for 
the  ith  node,  the  first  Ti(m-rl)  rows  being  for  the  nodes  on  the  frame  joints,  the  next 
2k  rows  for  the  nodes  at  the  cracks,  and  the  last  two  rows  for  the  nodes  at  the 
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excitation  and  response  locations.  The  procedure  to  establish  the  entries  of  the 
matrix  is  outlined  below. 


a.  The  diagonal  entry  u*  will  be  the  sum  of  the  Q  vaiues(Ecuation  2.31)  of  the 
beam  elements  adjoining  at  the  node  for  which  Klrchhoffs  current  law  is  being 
written.  If  this  node  is  at  a  crack,  then  is  given  by  the  Q  value  of  the  element 

ending  at  this  node  minus  the  A  value  for  the  corresponding  crack,  where 


0ftj(lTio)  ‘ 

And, 


The  diagonal  entries  of  B-  and  B7  will  thus  be  of  the  form  {Qt  —  .V). 


b.  The  ofi-diagonal  entry  tq;-  will  be  %  (Equation  2.31)  if  the  nodes  corresponding 
to  the  columns  i  and  j  of  the  matrix  are  linked  directly  by  the  k  th  beam  element  If 
these  two  nodes  are  linked  by  a  crack.  Le.,  they  are  the  nodes  on  the  two  sides  of  a 
crack,  then  tq;-  is  given  by  the  M  value  of  the  crack.  Otherwise  il}-  =  0.  Due  to 
symmetry.  iiy  =  u^.  The  ith  row  of  U0  in  Equation  (2.5S)  is  unchanged  relative  to  U3 
of  the  no-crack  ( k  =0)  case  if  the  ith  node  is  not  adjacent  to  a  cracked  element. 


2.  Tc-ssSe^lUh-to  Z  isaSris. 

The  matrices  Z  and  X  stem  from  the  kinetic  equations  and  shear  conditions.  The 
first  n  rows  in  them  represent  the  kinetic  equations  governing  sidesway  of  the  n 
floors.  The  nest  fc  +  1  rows  correspond  to  the  shear  continuity  conditions  at  the 
cracks  and  at  the  response  station.  The  last  row  results  from  the  shear  discontinuity 
at  the  escitation  location.  Z  consists  of  the  coefficients  of  deflections  y*  and  its 
entries  can  be  generated  as  follows: 

a.  The  diagonal  entry  z a  for  i=l,2,...n  will  be  the  sum  of  the  ft.  values  of  all  the 
girder  elements  on  the  ith  floor  plus  the  sum  of  the  P  values  (Equation  2.31)  of  all  the 
columns  adjoining  at  the  ith  floor.  If  a  column  adjacent  to  the  ith  floor  is  cracked, 
then  only  the  P  value  of  the  column  segment  nearest  to  that  floor  will  be  included,  ft. 
can  be  expressed  in  terms  of  the  characteristic  value  ft>  as  ft*  -  etftj.  If  the  frame  is 
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homogeneous  in  properties,  then  e*  -  Lk/  Lc.  The  diagonal  entry  z*  for  the  next  k 
rows  (Le.,  the  diagonal  entries  of  the  D  matrices)  will  be  the  sum  of  the  P  values  of 
the  two  beam  segments  on  the  two  sides  of  the  corresponding  crack.  z~  for  the  last 
two  rows(Le.,  the  diagonal  entries  of  the  N2  matrix)  will  be  given  by  the  sum  of  the  P 
values  of  the  two  beam  segments  joining  at  the  excitation  node  and  the  response 
node,  respectively. 

b.  The  ofl-diagonal  entry  Zi.i+l  of  Zc  will  be  the  negative  of  the  sum  of  the  P' 
values(Equation  2.31)  of  all  uncracked  columns  joining  ith  and  (i+l)th  floors. 

=  0  if  all  the  columns  between  those  floors  are  cracked,  z^-  =  0  for  n  &  j  ^.i+2 
and  Zij  -  Zji.  Za  is,  hence,  a  tridiagonal  symmetric  matrix.  If  any  one  of  the  columns 
adjacent  to  the  ith  floor  is  cracked,  then  z^  for  1  stsn,  n  <j  sin-rfe,  (Le.,  the 
entries  of  C.  ),  will  be  zero  except  for  that  (those)  value(s)  which  correspond(s)  to 
the  deflection(s)  at  the  crack(s)  nearest  to  the  ith  floor,  in  which  case  z i;-  is  equal  to 
the  negative  of  the  ?'  value  for  the  column  segment  which  Knk-«t  the  ith  floor  to  that 
crack.  If  acne  of  the  columns  adjacent  to  the  ith  floor  is  cracked,  then  ith  row  of  C. 
is  zero,  and  the  ith  row  of  Zg  is  unchanged  relative  to  Zq  of  the  no-crack,  k- 0,  case. 

c.  Zy-  =  0  for  li  isn,  n+k  <j  &  n+fc+2,  (the  entries  of  Nj),  where  k=k.+k7 
,  unless  there  is  an  excitation  and/ or  a  response  station  on  any  column(s)  adjacent  to 
the  ith  floor>  In  the  latter  case,  z.y  is  given  by  the  negative  of  the  P  value  for  the 
column  segment  linking  the  ith  floor  to  the  excitation  or  the  response  node. 

71  <  i  ^7t.+a:c  (the  entries  of  De)  will  be  zero  if  there  is  no  other 
crack  between  the  (i  — n)th  column  crack  and  the  floor  level  above  this  crack. 
Otherwise,  Zi,t>i  will  be  equal  to  the  negative  of  the  P‘  value  for  the  column  segment 
which  links  the  (i-n  )th  and  the  (i-n  4- i)th  cracks. 

*•  <  i  i  n+kj  (the  entries  of  Dy)  will  be  zero  if  there  is  no  other 

crack  between  (i~n.— fe»)th  girder  crack  and  the  wall  to  the  night  of  this  crack. 
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Otherwise  z.#<i+1  will  be  equal  to  the  negative  of  the  P’  value  for  the  girder  segment 
which  links  the  (i-n—k.)  th  crack  and  the  crack  on  its  right. 

Hence,  D-  and  D7  are  diagonal  if  at  most  one  crack  exists  on  each  beam 
element.  If  more  than  one  crack  exists  on  any  column  or  girder,  then  De  and/or  D7 
are  accordingly  tridiagonal  symmetric  matrices. 

f.  Zy-  for  n  <i  £n-rk,  n+k  <j  £n+i-r2  (the  entries  of  N\)  will  be  nonzero 
only  when  the  corresponding  crack  is  neighboring  an  excitation  or  response  node,  in 
which  case  z*y  is  equal  to  the  negative  of  the  P'  value  for  the  beam  segment  linking 
the  crack  to  the  excitation  or  the  response  node. 

3.  To  establish  the  X matrix: 

X  consists  of  the  coefficients  of  rotations  yi  in  equations  such  as  (2.50),  (2.53) 
and  (2.55).  A  square  submatrix  Xffc)  in  Equation  (2.56)  corresponds  to  rotations  of  the 
frame  joiizis  on  the  A:th  wail  (Figure  S). 

a.  The  diagonal  entry  z^  of  X'6',  which  stems  from  the  kinetic  equation  for  the 
ith  floor,  will  be  given  by  the  Q'  value  (Equation  2.31)  of  the  column  (column  segment) 
under  the  frame-joint  node  n(k  —  l)-ri  minus  the  Q  value  of  the  column  (column 
segment)  above  the  same  nods.  If  there  is  no  cracked  column  on  the  k  th  wall,  then 
=ik)  =  $*»(*-!)♦<  -  Q'n{ Jb-X)H+1  for  1  =  1 n-1  andr^*  = 

b.  The  ofi-diagonal  entry  2$$!  =  .  that  Is,  the  negative  of  the  W 

value  for  the  column  above  node  7i(k—  l)-K  for  i=l,...n  —  1  if  the  column  of  the  fcth 
•wall  between  the  ith  and  (i  +  l)th  floors  (Le.,  the  column  above  node  n{k—  l)+i)  has 
neither  a  crack,  nor  an  excitation,  nor  a  response  station  on  it.  Otherwise  =  0. 
and  that  entry  on  the  ith  row  of  Fe  or  In  which  corresponds  to  the  rotation  of  the 
crack  node  or  the  excitation  or  the  response  node  nearest  to  and  above  the  node 
n(<-l}-)-i  will  be  equal  to  the  negative  of  the  W'  value  of  the  column  segment  linking 
the  node  n(i—  l)-ri  and  the  said  crack,  excitation  or  response  node,  z ^  -  0  for 
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and  rW  =  sjf)  ;or  (The  entries  of  each  are  numbered 

independently  of  the  other  submatrices.)  X^  is  thus  diagonal  if  each  column  of  the 
fcth  wall  has  at  least  one  crack.  Otherwise  it  is  tridiagonaL  2?^  is  unchanged  relative 
iotk)  of  the  no-crack  (k  =0)  case  if  there  is  no  crack  on  the  fcth  walk 

If  the  column  below  node  n(k-l)+i  is  cracked,  or  has  an  excitation  or  a 
response  node,  then  that  entry  on  the  ith  row  of  Fe  or  Li  which  corresponds  to  the 
rotation  of  the  crack  nods,  or  the  excitation  or  the  response  node  nearest  to  and 
below  the  node  n(k— 1)-K  will  be  equal  to  the  W  value  of  the  column  segment  linking 
the  two  nodes.  Except  for  this  and  the  above  mentioned  cases,  entries  of  Fc  are  zero. 

c.  The  E  matrices  in  Equation  (2.5S)  couple  the  crack  nodes  with  the  frame-joint 
nodes.  Rows  n-rl  through  of  X  (Le.,  E  .  H  and  L[  matrices)  are  filled  in  as 
follows:  The  entry  on  the  matrix  column  which  corresponds  to  the  rotation  of  the 
node  on  the  left  (or  lower)  side  of  the  crack,  at  which  the  shear  continuity  condition  is 
being  written,  is  equal  to  the  Q'  value  of  the  girder  (or  column)  segment  which  links 
this  node  to  the  one  on  its  left  (cr  below  it).  The  entry  corresponding  to  the  node  on 
the  right  (or  upper)  side  of  the  crack  is  equal  to  minus  the  Q'  value  of  the  segment 
linking  this  node  to  the  node  on  its  right  (or  above  it).  These  two  entries  are  within 
the  submatrix  (or  H*).  The  entry  corresponding  to  the  node  on  the  left-hand  side 
of  (or  below)  the  left  (or  lower)  crack-node  is  equal  to  the  7 V  value  of  the  segment 
linking  these  two  nodes.  Finally,  the  entry  corresponding  to  the  node  on  the  right- 
hand  side  of  (or  above)  the  right-side  (or  upper)  crack-node  is  equal  to  the  negative  of 
the  V'  value  of  the  segment  Lin  long  the  two  nodes.  These  two  entries  can  be  within 
the  submatrices  E  H  or  Lj  depending  on  whether  there  are  one  or  more  cracks  on  a 
girder  (column)  and  on  whether  the  excitation  and  response  nodes  are  neighbors  with 
any  crack.  E.^  —  0  if  there  is  no  crack  on  the  kih  walk  Ej4  =  0  if  there  Is  no  crack  on 
any  of  the  girders  adjacent  to  the  /fcth.  walk 
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d.  The  diagonal  entries  of  Lj  are  given  by  the  difference  of  the  Q  values  of  the 
two  beam  segments  on  the  left  (or  lower)  and  right  (or  upper)  side  of  the  excitation 
or  response  nodes.  L3  is  a  diagonal  submatrir  if  the  excitation  and  response  nodes 
are  not  neighbors.  Only  those  entries  of  L2  corresponding  to  the  nodes  nearest  to  the 
excitation  or  the  response  node  are  nonzero  and  are  given  by  the  W  value,  or  the 
negative  of  it,  depending  on  whether  the  node  in  question  is  on  the  left-  or  right-hand 
side,  respectively,  of  the  excitation  or  the  response  node. 

4.  Corrections  for  the,  general  case: 

If  beam  properties  are  not  uniform  throughout  the  frame  (each  column  or  girder 
still  has  uniform  properties  within  itself),  then  the  following  corrections  are 
necessary: 

a.  Multiply  each  term  in  the  U  matrix,  except  the  A  terms,  by  the  corresponding 
T7t  value,  that  is,  Q?  and  Wp  by  rjpl. 

b.  Multiply  each  term  in  the  first  n  rows  of  the  X  matrix  by  the  corresponding 
T72  value,  that  is,  Qp'  and  Wp'  by  r,p3. 

c.  Multiply  each  term  in  the  first  tl  rows  of  the  Z  matrix  by  the  corresponding  t]s 
value,  that  is,  PP  and  Pp‘  by  rlpS  and  0-  by  rjrg. 

„  e,/L 

d.  Multioly  the  ivh  row  of  Z,  nflst  ^  n+jci-2,  by  A{  =  -  ■■■  — where  L  is  the 

Pv 

number  of  the  beam  element  with  the  crack  for  which  the  tth  row  represents  the 
shear  continuity  condition,  or  of  the  beam  element  with  the  excitation  or  the 
response  station  — ^ 

The  above  procedure  will  now  be  illustrated  with  no.  example. 

Example:  The  three  matrices  will  be  written  fcr  the  frame  in  Figure  11. 
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a.  Frans  without  cracks  (Fig.  11a)  with  different  bean  properties  (each  bean  Ls 
unifom  withia  itself). . 
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where,  for  convenience,  the  rotations  and  the  deflections  at  the  nodes  are  written 


ah  eve  the  corresponding  columns  cr  the  matrices.  It  is  noted  that,  fer  example. 
rj9i~ln  and  ??S2  =  7712  (2a.  2.51),  that  is, 'the  beam  segments  9  and  1  are  the 
segments  of  the  same  column  and,  hence,  have  the  same  properties. 


b.  Frame  with  two  cracks  (Fig.  lib).  The  frame  is  now  assumed  uniform  in 

properties  except  for  the  lengths  -of  the  individual  columns  and  girders.  Hence, 

>  __ 

rju  =  1  and  a £  a  =  L- — .  Matrix  U  is  obtained  as 

Po  Li} 


Q1+Q1C+Q& 

0 

0 

0 

1 

1  0 

0 

0 

r 

i 

1 

)’ 

I 

*1 

*10  ' 

0 

Qz+Qt 

0 

*8 

1 

1  0 

i 

0 

0 

0 

0. 

*2 

0 

c 

Qz*  Qz+  Gi 

0 

I  *7 

0 

0 

*s 

1 

1 

1 

0 

0 

0 

*a 

0 

i  0 

*4 

0 

0 

1 

1 

0 

0 

0 

0 

*7 

0 

i 

‘<?7-A4 

1 

0 

0 

1 

1 

0 

0 

0 

0 

0 

*4 

A* 

$4-A4 

0 

0 

i 

1 

1 

0 

0 

*a 

0 

0 

0 

,  -  _ 

i  0 

0 

if 

O) 

,cv 

Ac 

1 

1 

0 

0 

0 

0 

*a 

0 

!  0 

0 

* 

Qz~  Ac 

l 

l 

0 

0 

*1 

0 

0 

0 

0 

0 

0 

0 

t 

0 

*10 

*2 

0 

0 

0 

0 

0 

0 

1 

t 

0 

Qz~Q\z 

4 

(2.57) 

It  may  be  noted  that  the  flrst  four  elements  Ln  the  second  row  of  U  are  the  same  as 
these  in  the  second  row  in  Equation  (2.54)  (except  for  the  factors  This  row  is 

unchanged  since  node  2  is  not  linked  to  any  crack  node.  The  last  two  rows  'and  two 
columns  in  Equation  (2.54)  are  enlarged  by  filling  with  zeros  and  again  appear  as  the 


last  rows  and  columns  in  (2.57).  If  the  beams  had  different  properties,  then  each 
term  in  U  would  be  multiplied  by  the  corresponding  7]:  value  except  the  terms 


*  s  e,/Jo(i-rt3)  ‘ 


In  the 


terms,  the  sensitivity  numbers,  @1,  are  subscripted  by  the 
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original  number  of  tbs  column  or  girder  on  which  the  corresponding  cracks  are 
located-  Hence,  &4.  for  instance,  refers  to  the  crack  located  cn  the  column  numbered 
4  in-  the  uncracked  frame  (Figure  11).  It  is  also  noted  that  Qit  fly  P4,  etc.  for  the 
cracked  frame  have  different  values  frcm  the  ones  in  the  no-crack  case.  The  Z  and  X 
matrices  of  the  cracked  frame  are  subsequently  conjugated  as  : 
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(2.69) 


(2.5S) 


where  XH  in  Equation  (2. S3)  and  Equation  (2.69)  are  basically  the  same  since  there  is 
no  crack  on  the  first  wall  The  terms  -W4  and  W4  in  X/2)  in  Equation  (2.55)  are 
repiaced  by  zeros  in  Equation  (2.69),  since  nodes  3  and  4  are  no  longer  linked 
directly.  As  the  column  under  the  node  4  is  cracked,  the  node  4  is  now  linked  to  the 


node  T  by  the  column  segment  4.  Therefore,  appears  on  the  second  row  and 

second  column  (-which  corresponds  to  y?')  of  Fe.  The  third  row  of  X  reflects  the 
continuity  of  shear  at  the  column  crack  in  Figure  11b.  The  lower  crack-node  7  is 
linked  to  the  node  3  by  the  column  segment  7  and  the  upper  crack-node  T  is  linked 
to  the  node  4  by  the  segment  4.  Q'  and  W  values  are  accordingly  placed  on  the  third 
row.  The  construction  of  the  matrices  for  the  cracked  frame  is  now  complete. 
Numerical  results  for  the  frame  will  be  given  in  Chapter  3. 


CHAPTER  3 


THE  MODAL  FREQUENCY  THEORY  OF  FRACTURE  DAMAGE  DIAGNOSIS 
The  modal  frequency  theory  was  introduced  by  Ju  et  a!  [7,8]  to  utilize 
the  measurements  of  pre-  and  post-damage  modal  frequencies  in  detecting  and 
identifying  fracture  damage.  Each  pair  of  fracture  characteri sties  (7,e),  , 
the  kth  crack  intensity  and  location,  is  defined  by  a  damage  function 


(r,e)k  =  Gk(R,)  (3.1) 

wliere  R,  -  1  -  (■.,/<.,)  Is  the  modal  frequency  variation,  and  w  ,  are 
respectively  the  post-  and  pre-damage  frequencies  of  the  ith  mode.  For 
k-number  of  cracks,  therefore,  (2k  +  I)  number  of  measurements  of  frequency 
variations  are  needed  for  diagnostic  solutions.  The  solution  of  Equation 
(3.1)  yields  the  diagnosis  of  the  fracture  damage  in  the  structure.  The 
application  of  the  modal  frequency  theory  to  the  diagnosis  of  a  single  crack 
in  simple  structures  was  treated  b.y  Ju  et  al  [7,8]  with  the  matrix  method. 
The  method  in  general,  for  a  structure  of  N-beam  elements  with  k-cracks, 
requires  the  solution  of  a  (4N  +  2k)  by  (4N  +  2k)  matrix.  The  algorithm 
becomes  excessive  for  structures  of  large  number  of  beam  elements  and/or 
cracks.  The  numerical  difficulty  was  overcome  with  the  introduction  of  a 
generalized  theory  of  circuit  analogy.  The  theory  and  its  application  to 
multiple-crack  problem  and  to  complex  planar  structure  were  presented  by 
Akglin  and  Ju  [9,10]  and  summarized  in  [11].  The  basic  modal  frequency 
theory,  used  alone,  will  encounter  uncertainties,  arising  from  (1)  frequency 
crossover,  (2)  Inadequate  measurement  of  change  in  modal  frequencies  and  (3) 
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closely  packed  cracks.  The  uncertainty  of  frequency  crossover  can  be 
illustrated  with  two  modal  frequencies,  *>.  and  (without  loss  of  gene- 

•J 

rality  let  <■».  <  u» . )  at  undamaged  state.  After  crack  occurs  the  new  modal 

3 

frequencies  are  correspondi  ng  ».  and  TJ . .  If  the  crack  occurs  near  the 

J 

inflection  point  of  the  ith  modal  shape  but  if  the  point  happens  to  be  near 
the  maximum  moment  section  of  the  jth  modal  shape,  wt  may  not  differ  too 
much  from  i»..  It  is  conceivable  that,  for  the  new  modal  frequency,  we  may 
record  ».  <  tr. .  In  that  case,  not  knowing  the  actual  damage  configura¬ 
tion,  the  diagnostic  assemblage  of  frequencies  may  well  confuse  tr.  to  be 

V* 

the  new  frequency  of  the  ith  mode  and  xr.  to  be  the  new  frequency  of  the 
jth  mode.  The  phonomenon  is  a  frequency  crossover;  the  diagnostic  result  . 
may  be  disastrous.  Such  phenomenon  is  most  likely  to  occur  in  complex 
structures,  for  which  there  are  pairs  of  modal  frequencies  with  close 
values.  The  problem  of  frequency  crossover  cannot  be  handled  determinis¬ 
tically  by  the  modal  frequency  theory.  But  the  transmi ssibi 1 i ty  theory 
presented  in  the  sequal  should  alleviate  this  very  uncertainty.  The  rest  of 
the  chapter  will  address  the  uncertainties  of  multiple  cracks,  which  can 


only  be  resolved  by  probabilistic  means.  The  uncertainties  can  be  ade¬ 
quately  illustrated  with  a  simple  beam  structure  with  k-number  of  cracks. 

A  beam  with  k  cracks  can  be  represented  by  k  +  1  circuits  joined  by 
resistors  simulating  the  cracks.  The  unknowns  in  the  T-circuit  analogy  are 
the  moments  (currents)  at  the  cracks  and  at  the  ends  of  the  beam.  Hence, 
there  are  k  +  1  unknowns.  On  the  other  hand,  the  unknowns  in  the  n-circuit 


analogy  are  the  slopes  at  the  ends  of  the  beam  and  at  the  cracks.  There  are 
two  unknown  slopes  at  each  crack.  Hence,  the- order  of  the  system  is  larger 
with  n  circuits.  The  T-circuit  analog  derived  in  Section  2.2.1  is  found  to 
be  more  suitable  for  mul ti pi  e-crack  analysis.  This  section  develops  the 
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character!' sti c  frequency  equation  for  the  general  case,  establishes  the 
conditions  under  which  multiple  cracks  become  equivalent  to  a  single  crack, 
and  illustrates  the  inverse  problem  with  k  =  1.  The  beam  i-s  undamped  and  s 
and  all  other  variables  are,  therefore,  real. 


3.1.  Cantilever  Beam  with  Multiple  Cracks 

Figure  12  depicts  the  analog  circuit  for  a  cantilever  beam  with  k 
cracks.  Continuity  of  moments  at  the  cracked  sections  is  preserved  via  the 
continuity  of  electrical  currents  through  the  "crack  resistors",  -I/*-,.  The 
order  of  the  system  of  mesh  current  equations  is  k  +  1.  Namely, 


Z1  +  Z10)M1  "  Z'10M2  +  E11  "  0 


"Zi-l,0Mi-l  +  (Zi-l  +  Zi-1,0  +  h  +  Zi,0  ‘  1/ifi )Mi  "  E1-l,  2  +  E1,'l  =  0 


i  =  2,. . .k 


“ZkOMk  +  (Zk  +  Zk0  +  Zk+L  +  Zk+'1,0  '  i/’ck)Mk+l  "  Ek2  +  E  k+1,1  =  0 

(3.2) 

From  Equation  (2.21),  the  analog  voltage  sources  {E,  A  are: 
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where  is  the  deflection  of  the  i t h  cracked  section  and  yk+1  is  the 
deflection  of  the  free  end.  These  deflections  are  solved  for  by  imposing 
the  shear  boundary  and  continuity  conditions.  '  That  is, 


V12  *  V1+l,l  ,  i  =  1,2, ... ,k 


(3.4) 


where  the  first  subscripts  refer  to  the  numbers  of  the  elements  within  the 
beam  and  the  second  subscripts  1  and  2  denote  the  left  and  right  ends  of  the 
correspondi ng  element.  Substitution  of  the  last  two  equations  in  (2.17) 
into  (3.4)  yields 

(1/h)  Zy  =  Xm  (3.5) 

where  v  and  m  are  vectors  of  deflections  and  mesh  currents,  and  Z  and  X  are 


square  matrices  of  order  k+1  given  by 
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tion  (3.5),  provided  Z  is  non-singular. 


T‘ 

1  k+1 


matrix  and  X  is  an  upper  triangular 
y  can  be  solved  for  from  Equa- 
That  is, 


» 


Substitution  of  (3.7)  into  (3.3),  and  then  the  result  into  (3.2)  along  with 
the  analog  resistances  from  Equation  (2.20)  gives 

( XTZ_1X  -  U)m  -  0  (3.8) 

where 
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with  9j  being  the  sensitivity  number  for  the  jth  crack  given  by  Equa¬ 
tion  (2.40).  9 j  Is  based  on  the  total  length  of  the  beam. 

For  non-trivial  solutions  of  (3.8),  the  determinant  of  the  coefficient 
matrix,  which  is  symmetric,  must  vanish,  yielding  the  characteristic  equa¬ 
tion 


det(XTZ‘1X  -  U)  =  0  (3.10) 

When  the  crack  spacings  {e.}  ana  sensitivities  {0.}  are  known,  the  natural 

J  J 

frequencies  can  be  computed  from  Equation  (3.10).  In  computing  the  matrices 
in  equation  (3.10),  the  variables  ,  T1 ,  S ' ^ ,  T ' ^  (Equation  2.18)  per¬ 
taining  to  the  ith  beam  segment  are  computed  using  e^B  =  (L^/LJb  where  s  is 
based  on  the  total  length  of  the  beam. 
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3.2.  Damage  Diagnosis  with  a  Single  Crack 

Damage  diagnosis  using  the  present  model  is  accomplished  with  the 
knowledge  of  frequencies  after  the  damage  has  occurred.  Since  the  charac¬ 
teristic  equations  of  structures  are  in  terms  of  dimensionless  charac¬ 
teristic  values  {Si  from  the  measured  frequencies  {8}  can  be  computed  from 


s  =  L 


(3.11) 


When  three  measurements  of  frequency  are  available  for  the  case  k  =  1 
(single  crack),  exact  location,  e,  and  severity,  9,  of  the  crack  can  be 
determined.  The  procedure  will  be  illustrated  with  the  cantilever  beam 
problem  developed  in  the  previous  section.  The  coefficient  matrix  in  equa¬ 
tion  (3.8),  with  k  =  1,  is  of  the  form 


rt-1*  -  «  ■  [£  h^f29B]  ’  (3.12) 

where  the  {8}  values  in  equation  (3.12)  are  the  post-damage  values  computed 
from  equation  (3.11).  The  damage  characteristics  e  and  9  are  to  be  deter¬ 
mined  from  the  determinant  of  H,  which  can  be  written  as 


det(H(s,e,9) ]  =  det[H(s,e,0) ]  -  208 h-^  =  0  (3.13) 

where  e  =  L^/L  is  the  normalized  crack  location.  It  is  noted  that  when  8 
assumes  the  pre-damage  values  b  ,  det[H(B,e,0)]  is  equal  to  0;  that  is,  this 
term  is  the  characteristic  equation  for  the  undamaged  (9  =  0)  beam.  There 
are  three  equations  emerging  from  Equation  (3.13)  for  the  three  known 
characteri sti c  values.-  A  numerical  code  for  damage  diagnosis  has  been 
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developed  in  which  e  is  varied  through  the  range  (0,1).  8  is  computed,  for 
a  given  e,  using  the  first  known  characteristic  value  in  Equation  (3.13), 
namely, 


.  detfH(B(1),e,0)l 

‘  ’  '  2B<% 


(3.14) 


where  the  superscript  on  8  denotes  which  characteristic  value  is  used.  Then 
with  the  second  known  characteristic  value,  8^, 


det(H(a^  ,e,9)  ]  =  0 


(3.15) 


where  the  value  of  8  is  now  substituted  from  Equation  (3.14).  A  zero¬ 
searching  routine  is  used  to  find  the  roots  of  Equation  (3.15)  which,  in 
general,  yields  multiple  solutions  for  e.  8^  can  then  be  used  to  locate 
the  crack.  With  e  known,  0  is  computed  from  (3.14). 

In  the  case  when  the  beam  deviates  somewhat  from  the  Bernoul  1  i-Euler 
theory  or  when  the  material  properties  are  uncertain,  computing  {Si  from 
Equation  (3.11)  may  not  be  suitable.  The  knowledge  of  modal  frequencies 

prior  to  the  damage  is  then  necessary.  Since  u  and  u>u  are  proportional  to 
2  2 

a  and  s|],  respectively,  where  and  bu  are  the  values  for  the  pre-damage 
structure,  and  <o  and  a  are  the  post-damage  values,  the  following  relation 
holds: 


(3.16) 


It  is  assumed  in  (3.16)  that  the  proportionality  constant  between  u»u  and 
does  not  change  after  the  damage  has  occurred.  The  characteri stlc 
equations  do  not  involve  any  material  properties  and  Bu  can  be  computed  from 
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the  characteristic  equation  for  the  undamaged  structure.  With  u>u  and  w 
measured  in  the  field,  {8}  can  then  be  computed  from  Equation  (3.16)  and  the 
above  procedure  can  again  be  employed  to  determine  e  and  9.  Only  the  beam 
length  L,  the  Poisson's  ratio  v,  and  the  slenderness  ration  (b/L)  are 
required  to  determine  the  actual  crack  location  L,  and  the  crack  depth  a 
(Equation  2.40).  Since  \>  is  nearly  the  same  for  most  metals,  the  theory 
becomes  independent  of  the  specific  material  of  the  structure  if  the 
material  is  a  metal,  if  the  damage  is  presented  in  terms  of  and  a,  and  if 
(8)  are  computed  from  Equation  (3.16).  The  theory  is  totally  independent  of 
the  material  properties  and  dimensions  if  the  damage  is  presented  in  terms 
of  e  and  9.  The  data  are  more  conveniently  presented  versus  the  relative 
frequencies  (<o/wu)  or  the  relative  frequency  changes  (1  -  w/wu)  rather  than 
versus  the  absolute  frequencies,  whether  (si  are  computed  from  Equa¬ 
tion  (3.11)  or  (3.16).  In  this  report,  relative  changes  are  chosen  as  the 
means  to  convey  the  numerical  results. 

In  practice,  the  responses  of  structures  deviate  from  the  Bernoulli- 
Euler  theory.  It  is  therefore  more  accurate  to  compute  (8)  from  Equa¬ 
tion  (3.16).  Equation  (3.16),  in  effect,  adjusts  the  parameters  of  the 
specific  structure  such  that  the  structure  frequencies  match  those  predicted 
by  the  Bernoulli-Euler  theory. 

In  practice,  especially  for  structures  whose  frequencies  are  close  to 
each  other,  the  reduction  in  a  certain  frequency  due  to  damage  may  be  small’ 
while  the  reduction  in  the  next  frequency  may  drop  below  that  of  the  former. 
This  phenomenon  is  called  crossover.  When  crossover  occurs  and  when  one  is 
not  aware  that  it  has  occurred,  the  correspondence  established  between  the 
pre-damage  and  the  post-damage  values  of  the  frequencies  will  be  in  error. 
If,  in  such  a  case,  one  used  Equation  (3.16),  which  involves  the  ratio  of 
the  pre-  and  post-damage  values  of  the  frequencies,  the  computed  charac- 
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ter 1  Stic  values  will  be  incorrect.  On  the  other  hand,  if  Equation  (3.11), 
which  involves  the  structural  properties,  is  used,  a  knowledge  of  the  pre¬ 
damage  frequencies  is  not  required.  Hence,  it  becomes  immaterial  which 
pre-damage  frequency  a  certain  measured  post-damage  frequency  corresponds 
to;  the  important  point  is  that  there  are  some  frequency  values  available 
wh '  h  satisfy  the  post-damage  characteristic  equation  (Equation  3.13).  If 
the  structural  properties  are  uncertain,  they  can  be  identified  by  measuring 
one  frequency  prior  to  the  damage.  That  is,  the  constant  K  in  the  relation 
Bu  =  can  cletern1'inecl  measuring  one  <uy  and  computing  the  corre¬ 

sponding  bu  from  the  pre-damage  characteristic  equation.  The  same  constant 
can  later  be  used  to  compute  all  the  needed  post-damage  frequencies  u  by 
means  of  b  =  Kvd.  In  this  case,  however,  the  R*  values  cannot  be  used  to 

*J 

present  the  data,  since  the  ordering  of  the  post-damage  frequencies  accord¬ 
ing  to  the  magnitudes  of  their  values  does  not  correspond  to  the  pre-damage 
ordering,  unless  the  fact  that  crossover  has  occurred  is  known. 

3.3.  Uncertainties  of  Closely  Packed  Multiple  Cracks 

Under  certain  conditions,  the  effect  of  multiple  cracks  is  not  dis¬ 
tinguishable  from  that  of  a  single  crack.  In  this  section,  these  conditions 
are  investigated.  For  this  purpose,  equivalence  of  two  cracks  on  a  simply- 
supported  beam  to  a  single  crack  is  first  established  analytically.  Numer¬ 
ical  results  for  larger  numbers  of  cracks  are  then  presented.  It  is  esta¬ 
blished  that,  when  equivalence  holds,  solution  to  the  inverse  problem  of 
damage  diagnosis  cannot  differentiate  between  single  and  multiple  cracks. 

Figure  L 3a  depicts  a  simply-supported  beam  with  a  single  crack  and  i^s 
T-circuit  analog.  The  mesh  equation  together  with  the  Equations  (2.20)  and 
(2.21)  yield 
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-(S1  +  S2  +  20S)M  -  (Sr  +  S2.)  yl h  =  0 


(3.17) 


where  M  and  y  are  the  resisting  moment  and  deflection  at  the  cracked  sec¬ 
tion,  respectively.  The  continuity  condition  =  at  crac*<  aiiows 
y  to  be  solved  for  in  terms  of  M.  The  characteristic  equation  for  a  beam 
with  a  single  crack  is  thus 


(Si  +  St) 

Wfl-%  -  <S1  +  S2  *  2a»>  *  0 


(3.18) 


For  the  case  of  two  cracks,  with  reference  to  Figure  13b,  the  characteristic 
equation  is  obtained  in  the  form  of  Equation  (3.10)  with 


-Si-S2  r2  1  .  .  |VS2  T2 

Ti  -S2-S3  ’  T2  S2+S3 


(3.19) 


S'i +$2+20-1 B 


^2+^3+^2® 


When  it  is  assumed  that  the  spacing  between  the  cracks  is  sufficiently 
small,  that  is,  e^  =  l^/l  =  e;  the  following  approximations  are  obtained. 


cosh  e,B  a  1  +  f  ■  ,  sinh  e^B  «  eB  +  ^ 


(3.20) 


cos  e^B  »  1  -  5 


f  r  P  1 

sin  e,B  *  eB  -  r 
c  o 


It  then  follows  from  Equation  (2.18)  that 

<  a  T  #  l!  S'  a  f  ss  — 
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(3.21) 


After  (3.21)  is  substituted  into  (3.19)  and  the  necessary  operations  are 


performed,  the  matrix  of  coefficients  in  (3.8)  is  obtained  as 


*7  Jit 


The  above  analysis  is  valid  at  points  where  Sp  Sj,  $p  S3  are  of  sna^e^ 
order  of  magnitude  than  i/rB.  The  characteristic  equation  for  the  two  crack 
case  is  thus 


det(XTZ_IX  -  U)  = 


4S ,  S- 


1  r  /  c  <  j.  C  '  >i  2 

2 


(3. 2d) 


-  ( S|  +  S-j)^S|  +  +  2(0|  +  0 2 ) S J  }  +  0(1/ eS )  = 


For  small  values  of  e.;,  a  comparison  of  Equations  (3.18)  and  (3. 2d)  reveals 
that  $p  S,  , ,  S.^,  $2.  in  Equation  (3.18)  are  nearly  equal  to  $p  S'p  Sp 
S'p  respectively,  in  Equation  (3. 2d).  The  difference  between  the  charac¬ 
teristic  equations  for  one-  and  two-crack  cases  is  of  the  order  of  the  pro¬ 
duct  028.  (The  difference  is  proportional  to  b  and  is,  therefore,  larger 
for  higher  modes.)  The  effective  sensitivity  number  for  the  two  cracks 
approximates,  as  shown  in  (3.18)  and  (3.24),  the  sum  of  the  individual 
sensitivity  numbers,  that  is,  9  =  8,+  9?.  Thus,  closely  spaced  multiple 
cracks,  in  general,  cannot  be  differentiated  from  a  single  equivalent  crack. 
Numerical  data  for  different  crack  conf igurations  indicate  the  same  result. 

Similar  results  are  present  for  a  cl amped-cl amped  beam.  Numerical 
results  for  a  cantilever  beam  indicate  that  the  first  few  characteristic 


values  of  a  beam  with  closely  spaced  multiple  cracks  are,  in  general,  close 
to  those  of  a  beam  with  an  equivalent  single  crack  whose  sensitivity  number, 
9  ,  aP°rox'ima1:ely  equal  to  the  sum  of  the  individual  sensitivity  numoers 
of  the  cracks  on  the  original  beam.  Furthermore,  the  location  of  the  equi¬ 


valent  crack  is  generally  within  the  region  where  the  group  of  cracks  is 


focated. 


In  conclusion,  equivalence  of  cltfsely  spaced  cracks  to  a  single 


crack  implies  that,  in  the  process  of  damage  diagnosis,  it  is  imcossible  to 


distinguish  between  closely  spaced  multiple  cracks  and  a  single  crack. 
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Figure  14  typically  illustrates  the  lower  limit  of  crack  spacing  for  a 

cantilever  beam  when  uniformly  spaced  multiple  cracks  become,  as  a  whole, 

indistinguishable  from  a  single  crack.  In  Figure  3,  as  in  subsequent 

illustrations,  each  crack  assumes  a  sensitivity  number  9^  of  value  0.174 

which  is  taken  to  be  a  major  crack  in  a  beam  of  slenderness  ratio  of  0.05 

(corresponding  to  a  relative  crack  depth  y  =  a/b  of  0.6).  The  equivalence 

is  based  on  equal  values  in  the  first  two  modal  frequency  changes  and  a 

tolerance  of  O.l'o  in  the  third  mode  characteri  sti  c  values,  i.e., 

.  i  B gg  ^  -  8^1  <  0.001,  where  the  subscripts  eq  and  ac  denote  the 

values  for  the  equivalent  single  crack  and  the  actual  damage  configuration, 

1 3 1 

respectively.  If  the  tolerance  level  in  is  increased,  the  curves 

in  Figure  14  will  shift  upward.  As  the  number  of  cracks  increases,  they 
must  be  more  closely  spaced  to  be  representable  by  a  single  crack.  The 
lower  limit  of  spacing  depends  on  where  the  group  of  cracks  is  located.  At 
the  built-in  end  and  at  =  0.6,  more  widely  spaced  cracks  can  become 
indistinguishable  than  at  =  0.1.  The  lower  limit  of  spacing  is  dependent 
also  on  9 ^  values.  For  example,  when  e^  =  0.1,  k  =  2  and  01  =»  92  =  0.01 
( correspondi ng  to  y  =  0.2  for  b/t  =  0.05),  the  smallest  crack  spacing  for 
which  the  double-crack  damage  becomes  indistinguishable  from  a  single-crack 
damage  is  e^  =  0.025. 

Figure  15  illustrates  the  relative  frequency  changes,  Rj,  for  a  canti¬ 
lever  beam  with  two  closely-spaced  cracks  of  9  =  .174  each.  A  very  similar 
set  of  curves  is  obtained  if  (R.)  is  plotted  for  a  cantilever  with  one  crack 

vJ 

of  9gq  =  .348.  For  any  given  number  of  cracks  on  a  cantilever  beam,  the 
largest  decrease  in  the  fundamental  frequency  occurs  when  all  the  cracks  are 
groupd  at  the  built-in  end.  On  the  other  hand,  for  9 ^  >  .3,  the  greatest 
change  in  the  second  frequency  is  observed  at  a  location  4).  Such  infor¬ 
mation  can  be  utilized  to  set  rough  guidelines  for  damage  diagnosis. 
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It  is  clear  that  probabilistic  methods  must  be  resorted  to  in  diagnosis 
of  multiple  cracks  when  an  insufficient  number  of  frequency  measurements  are 
available.  Nevertheless,  some  qualitative  conclusions  may  be  reached  with 
the  help  of  Figure  14  which  is  independent  of  material  oroperties.  In  par¬ 
ticular,  when  a  solution  has  been  obtained  for  a  single  crack  and  a  decision 
has  to  be  made  between  a  single  crack  and  a  group  of  closely-spaced  less 
severe  cracks,  one  may  argue  that  formation  of  a  single  crack  is  more 
likely.  In  other  words,  a  small  mode  I  type  crack  is  likely  to  propagate 
under  bending  rather  than  other  cracks  forming  nearby.  On  the  other  hand, 
one  severe  and  several  minor  cracks  may  exist  distributed  along  a  structure 
in  which  case  the  objective  is  to  be  able  to  diagnose  the  major  crack. 
Guidelines  may  be  established  qualitatively  based  on  the  frequency  change 
curves  for  a  cantilever  beam.  Figure  15,  in  which  case  we  may  propose  that 

a.  If  R^  is  considerably  smaller  than  and/or  R^,  then  the  (major 
crack(s))  is  at  a  normalized  distance  greater  than  -0.45  from  the  built-in 
end  (Fig.  15).  In  addition,  (i.)  if  R^  is  considerably  larger  than  R^.  the 
crack  is  located  around  0.45-0.65  relative  to  the  built-in  end  (see 
Section  6.0  below);  (ii.)  if  R2  is  larger  than  R£,  the  crack  is  at  a  dis¬ 
tance  greater  than  -0.7  or;  (iii.)  if  R2  and  R^  are  comparable,  then  there 
may  be  one  major  crack  at  0.65-0.70,  or  two  major  cracks  one  each  in  the 
peak  regions  of  the  R2  and  R2  curves. 

b.  If  Rj  is  significantly  smaller  than  both  R^  and  R^,  the  crack  is  in 
the  region  0.2-0. 3,  or  there  may  be  a  (major)  crack  at  0.2-0. 3  and  another 
(major)  one  at  a  distance  greater  than  0.8,  the  latter  being  more  likely  if 
Rj  is  greater  than  R^. 

c.  If  R^  is  significantly  larger  than  both  R2  and  R^,  and  R^  and  R^ 
are  rather  small,  the  crack  is  at  0.1-0. 2. 


jL 
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d.  If  the  values  of  R^,  R.,  and  R2  are  comparable,  several  possibil¬ 
ities  exist.  To  list  a  few,  there  may  be  a  number  of  cracks  distributed 
over  the  beam  (Section  3.5  below);  the  crack  may  be  at  0.35-0.40;  or  there 
may  be  a  crack  at  0.0-0. 3  and  a  few  others  at  locations  greater  than  -0.5. 

In  the  present  study,  only  the  first  three  natural  frequencies  are 
assumed  measurable.  The  peak  locations  of  R.^  and  curves  in  Figure  15  are 
weak  functions  of  0  increases.  For  example,  as  9  is  increased  from  0.01  to 
0.50  (effective  value  for  a  few  closely  spaced  cracks),  the  peak  of  the  R^ 
curve  shifts  from  the  location  0.530  to  0.563  while  the  second  peak  of  the 
curve  shifts  from  0.710  to  0.760.  Similar  guidelines  for  different 
structures  can  be  established  based  on  minimum  and  maximum  frequency  changes 
although  the  procedure  may  become  tedious  for  more  complicated  structures. 


3.4.  Uncertainty  From  Inadequate  Measurements 

To  illustrate  the  uncertainties  involved  in  damage  diagnosis  when  the 
three  R.  values  are  comparable  in  value,  the  following  relative  frequency 
changes  are  assumed  to  have  been  computed  from  the  measured  frequencies  of  a 
cantilever  beam:  R^  =  .0579,  R-^  =  .0586,  R^  =  .0591.  It  is  to  be  determined 
whether  the  damage  mainly  consists  of  one  major  crack. 

The  (actual)  characteristic  values  are  computed  from 


,m  *  „»>  (1 .  r.,1/2 


ac 


j 


where  {8^}  are  the  undamaged  characteristic  values,  the  first  three  of 
which  for  a  cantilever  are  (1.8751,  4.6941,  and  7.8543).  From  (22),  in 
particular,  8 =  7.6192.  The  first  two  actual  characteristic  values 
yield  one  solution  for  a  single  crack  located  at  egC]  =  .372  with  9gq  =  .128 
(corresponding  to  a  relative  crack  depth  of  7  =  .55  for  a  slenderness  ratio 
of  b/L  =  .05).  These  equivalent  values  are  then  used  in  the  forward  problem 
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with  k  =  1  sqs.  3.6,  9,  10  (10,  11,  14,  16)  to  determine  the  third  charac¬ 
teristic  value  due  to  the  equivalent  crack.  The  computation  thus  yields 

=  7.6176.  Then,  -  8^'  =  .0016.  If  we  assume  that  the 

/  3  \ 

accuracy  in  s'  '  is  not  better  than  ±.0016,  the  equivalent  single  crack  . 
solution  may  be  accepted  as  the  true  damage  diagnosis.  On  the  other  hand, 
the  given  R.  data  were  actually  generated  in  a  forward  problem  with  ten 

v 

cracks  of  equal  intensity  (0  was  taken  to  be  .01  fi.e.,  y  -  .1}  for  all  the 

cracks)  located  such  that  e^  =  .001,  ei  =  .094  for  i  =  2,3 . 10.  Thus,  in 

this  case,  even  though  the  actual  and  equivalent  S'  ’  values  match  closely, 
a  firm  diagnosis  cannot  be  reached. 

3.5.  Effect  of  Minor  Cracks  Associated  with  a  Major  Crack. 

For  the  purpose  of  illustrating  the  effect,  it  is  sufficient  to  use  the 
forward  problem  formulation  in  which  the  locations  and  intensities  of  the 
cracks  are  known.  The  character i sti  c  values  are  determined  from 

Equation  13.10).  The  relative  frequency  changes  are  then  determined  from 
the  characteristic  values,  namely. 


a.  First  the  analysis  will  start  with  one  crack  (Fig.  16a),  i.e., 
k  =  1,  at  e  =  .54  (in  the  peak  region  of  the  curve.  Figure  15)  with 
0  =  .174  (y  =  .6).  The  relative  frequency  changes  are  then  comDuted  as 
Rj  =  .029,  R2  =  .137,  R3  =  .014. 

b.  Now,  in  addition  to  the  crack  in  part  (a),  there  are  five  more 
cracks  of  intensity  0  =  .0026  (y  =  .1)  each,  distributed  on  the  cantilever 
such  the  =  s2  =  e3  =  e5  =  -1*  =  .24,  and  e^  =  .09  (Fig.  16b).  The 
actual  characteristic  values  are  computed  as  3^  =  1.840,  =  4.351,  and 
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3^  =  7.768  from  which  =  .037,  R^  =  -141,  R^  =  *022.  If  thesa  values  of 

R^  ware  obtained  from  actual  measured  frequencis  from  which  the  damage  were 

to  be  diagnosed  (the  inverse  problem),  3  - ^  and  B ‘ ^  ■  are  used  to  find  a 

solution  for  a  sinqle  crack,  which  is  (e  .91*  (.51  .133).  This  sinale 

'  eq  eq-  ' 

crack  would  produce  B^-  =  7.845  or  R0  =  .0024. 

eq  3eq 

c.  The  same  case  as  in  (b)  except  0.  are  changed  to  .01  (/  =  .2)  with 
the  major  crack  remaining  the  same,  i.e.,  9^  =  .174.  Hence,  the  deoths  of 
the  minor  cracks  are  doubled.  The  corresponding  frequency  changes  in  the 
forward  problem  are  R^  =  .058,  R^  =  .151,  R^  =  .043.  If  the  same  R^  and  R^ 
are  in  turn  used  in  the  inverse  problem  of  damage  diagnosis  with  k  =  1,  we 
obtain  (eeq,  0qq)  =  (.47,  .224),  B^  =  7.830,  R3eq  =  .0064. 

The  above  cases  indicate  that  the  value  of  R3  did  not  change  signifi¬ 
cantly  wnen  minor  cracks  were  added  to  the  beam.  If  this  were  an  actual 


damage  diagnosis  problem,  in  which  case  only  (a^)  and  (R.)  would  be  known 

.  J 

from  the  measured  data,  then  a  decision  would  have  to  be  made  about  the 

comouted  (e  .  8  1.  In  case  (b),  R0  =  .022  whereas  R-,  =  .0024.  In  cas 

ea  ea  3  ieq 

(c),  R-.  =  .043,  but  R0 _  =  .0064.  Hence,  based  onlv  on  comoarison  of  R-, 

3  3eq  •  3eq 

with  the  actual  values  of  R,,  the  solutions  for  (e  ,  0  )  would  be  re- 

3  '  eq  eq  ’ 

jected.  Nevertheless,  the  equivalent  damage  parameters  in  both  cases 
closely  identify  the  major  crack  located  at  the  distance  of  .54  from  the 
built-in  end  with  0  =  .174.  When  there  is  no  minor  crack  on  the  beam,  as  in 
case  (a),  (e  ,  ©eq)  would  be  computed  as  (.54,  .174).  The  effect  of  the 
minor  cracks  on  (egq,  0gq)  is  thus  seen  to  be  small.  Hence,  it  is  concluded 


that,  when  the  given  R.  values  exhibit  a  pattern  such  as  in  this  example, 

v 

the  discrepancy  in  R,  and  R,„„  or  in  B^  and  can  be  iqnc  ed 

3  3eq  ac  eq 

when  the  interest  is  in  diagnosing  the  major  crack. 
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3.6.  Peak  Modal  Response 

It  is  to  be  demonstrated  that,  when  multiple  cracks  exist  on  a  beam,  a 

crack  which  is  located  in  the  peak  region  of  an  R,  curve  affects  the  corre- 

J 

sponding  frequency,  the  most. 

J 

a.  Given  k  =  I ,  e^  =  0.001,  9-j  =  0.174.  The  frequency  changes  are 

computed  as  (the  forward  problem)  R-^  =  0.234,  =  0.162,  R^  =  0.125. 

b.  A  second  crack  is  added  to  case  (a)  such  that  e2  =  0.546  (i.e., 
near  the  peak  region  of  the  R2  curve  of  a  single  crack  as  in  Figure  15)  and 
92  =  0.174.  Then  R,  =  0.246,  R2  =  0.272,  R3  =  0.162. 

c.  A  third  crack  is  added  to  case  (b)  such  that  e^  =  0.185  (i.e.,  near 

a  peak  region  of  the  R3  curve  of  the  single  crack  in  Fig.  15)  and  9^  =  0.174 

[parameter  values  in  cases  ( a- b )  are  preserved].  Then  Rj  =  0.247, 

R2  =  0.273,  R^  =  0.167.  In  case  (b),  the  greatest  decrease  relative  to  case 
(a)  occurred  in  the  second  frequency ,  whereas  in  case  (c)  it  was  the  third 
frequency  that  suffered  the  largest  decrease  relative  to  case  (b)  although 
this  decrease  was  insignificant  for  9^  =  0.01.  The  example  illustrates  the 
effect  of  crack  location  and  intensity  on  the  frequencies  and  confirms  the 
usefulness  of  curves  such  as  Figure  15  to  diagnosis  of  damage  as  outlined  in 
the  guidelines  presented  earlier. 

It  should  be  noted  that  the  sensitivity  number  9  is  a  measure  of  how 
sensitive  the  natural  frequencies  are  to  given  crack  depth  and  location  as  a 
function  of  the  slenderness  ratio.  Oecrease  in  natural  frequencies  is 
greater  for  larger  values  of  9.  Between  two  cracked  beams  with  the  same 
crack  location  and  relative  crack  depth,  9  for  the  more  slender  one  will  be 
smaller  giving  rise  to  smaller  Rj  values.  For  example,  for  7  =  0.6,  (and 
v  =  0.3),  9  =  0.0697  when  b/l  =  0.02  and  9  =  0.0174  when  b/L  =  0.005.  Thus, 
in  practice  it  is  relatively  harder  to  diagnose  damage  in  slender  beams.  On 


the  other  hand,  two  beams  of  the  same  slenderness  ratio  and  Poisson's  ratio, 
but  or  difrerent  materials,  experience  the  same  relative  decrease  in  fre¬ 
quencies  for  -the  same  crack  location  and  relative  depth. 
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CHAPTER  4 

TRANSMISSIBILITY  THEORY 'OF  DAMAGE  DIAGNOSIS 


In  this  chapter,  a  new  method  is  proposed  for  diagnosing  damage  in 
large  structures.  The  method  utilizes  the  changes  in  transmissibi 1 ities  at 
several  locations,  called  the  response  stations,  on  the  structure.  The 
method  is  developed  by  computing  the  transmissibi lity  changes  for  known 
crack  configurations  and  comparing  the  magnitudes  of  the  changes  at  differ¬ 
ent  response  stations.  The  equations  derived  in  Chapter  2  through  the 
electrical  analogy  of  a  damped  frame  structure  are  used  in  the  computations. 


4.1.  Transmissibi lity  and  Relative  Transmissibi lity  Change 

In  this  study,  transmissibility  is  defined  as  the  magnitude  of  the 
ratio  of  acceleration  at  a  response  station  to  the  force  applied  at  an  arbi¬ 
trary  location  on  the  structure  (Fig.  17).  Namely, 


T  ^  ^  (4.1) 

where  yt  is  the  complex  deflection  at  the  response  station  of  interest. 
When  the  force  is  sinusoidal,  that  is,  F=Fme^wt,  the  deflection  becomes 


1  [u*t— 0] 
ym 


(4.2) 


where  y  is  the  complex  modal  shape,  ym  is  the  amplitude  of  transverse  de¬ 
flection,  and  0  is  the  phase  angle.  ym  and  0  are  functions  of  location  and 
excitation  frequency.  0  is  zero  for  an  undamped  structure.  Transmissi¬ 
bility  for  sinusoidal  excitation  is  then  given  by 
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Figure  1 7.  A  three-story  four-span  frame. 
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(4.3) 


The  amplitude  of  deflection,  ym,  at  the  response  station  can  be  computed  from 


equation  (2.63)  for  any  plane  frame  with  and  without  a  crack,  when  the  crack  location 
and  the  sensitivity  number,  0.  are  known.  In  practice,  transmissibilities  can  be 
measured  periodically  or  after  a  strong  motion  event  and  compared  to  baseline 
values  to  diagnose  the  damage.  To  facilitate  a  diagnosis  method,  the  relative 
transmissibility  change  is  introduced  as  follows: 


_  7  -  T,  Tf _ 

*r%  ”  rp  1 

<*0  «4  fl 


(4.4) 


where  the  subscripts  c  and  o  indicate  the  transmissibilities  with  and  without  the 
crack,  respectively,  TTnen  the  structure  is  excited  at  the  same  frequency  and  with 
the  same  force  magnitude  before  and  after  the  damage.  R?  becomes 


^1 


An  undamped,  structure  has  node  points,  that  is,  points  with  zero  transverse 
deflection  at  ail  times.  If,  for  such  a  structure,  the  response  station  is  at  a  node  point 
of  the  undamaged  modal  shape,  then  =  0  .  If,  on  the  other  band,  it  is  at  a  node 
point  of  the  post-damage  modal  shape,  then  y„_  =0.1 n  general,  then, 

-  1  <  Rt  <  “  (4.6) 

Damped  structures,  in  general,  do  not  have  true  node  points.  Every  point  on  the 
structure  deflects  somewhat  during  a  period  of  the  excitation.  R?,  however,  still  has 
minimum  2nd  maximum  values  with  respect  to  the  location  of  the  response  station. 
To  facilitate  referral  to  such  locations,  the  following  definition  is  introduced.  A 
oseudo-node  point.  ?N?,  of  a  damped  structure  is  a  point  with  locally  minimum 
amplitude  of  transverse  deflection.  Thus.  Rt  is  maximum  at  the  FN?s  of  the 
undamaged  structure  and  minimum  at  those  of  the  damaged  structure. 

Relative  trarsmissibility  change  is  a  function  of  the  excitation  location  and 
frequency,  the  location  of  the  response  station,  and  the  crack  location  and  severity, 
as  well  as  the  structural  parameters.  It  is  desirable  to  have  large  values  of  relative 
transncissibility  changes  for  as  small  cracks  as  possible.  This  may  be  accomplished  by 
optimizing  the  parameters  which  can  be  controlled. 

4.2.  Optimum  Parameter  Values 

It  is  found  that  exciting  a  frame  structure  cn  a  column  rather  than  cn  a  girder  is 
more  advantageous  for  diagnosis  purposes.  A  3x4  structure  (Le.,  n  =3  ,  m.=4)  is 
shown  in  Figure  17 .  The  frames  studied  throughout  this  chapter  are  uniform  unless 
stated  otherwise;  namely,  all  of  the  beam  elements  which  make  up  the  frames  are  of 
equal  length  and  have  the  same  material  and  geometric  properties.  Figure  IS 
illustrates  the  variation  of  the  relative  transmissibiiity  change  at  a  response  station 
as  a  function  of  the  excitation  location  which  is  varied  along  the  leftmost  wall  cf  the 
frame.  Each  third  of  the  abscissa  corresponds  to  one  of  the  columns  cn  that  wail,  with 
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N=3,  M=4.  CRACK  ON  17  AT  .1,  THETA=.G24,  EXC.  ON  WALL  1 
RESPONSE  STATION  ON  IS  AT  .51,  UNDAMFED  EETA=.S6 


Figure  13. Relative  transmissibility  change  vs  excitation  location. 


distance  measured  from  the  anchor  up.  The  response  station  is  or.  girder  IS  at  a 
normalized  distance  of  Q.S1  measured  from  the  left  end.  The  crack  is  on  girder  17  at 
a  distance  of  0.1  from  the  left  end  of  that  girder.  Sensitivity  number,  9,  for  the  crack 
is  0.024  (Eq.  2.40).  Damping  factors,  o,  (Sq.  2.2)  for  the  structure  is  taken  to  be  0.01. 
The  undamped  /?  value,  &,  which  is  a  function  of  the  excitation  frequency  (Eq.  2.S),  Is 
0.38,  nearly  the  same  as  the  first  undamped  characteristic  value  of  the  structure. 
The  figure  indicates  that  the  relative  transmissibiiity  change  for  a  given  crack  Is 
largest  when  the  excitation  Is  applied  near  a  frame  joint. 

Figure  19  illustrates  the  dependence  of  the  relative  transmissibiiity  change,  Rr, 
on  the  excitation  frequency  expressed  in  terms  of  the  undamped  value.  Other 
variables  have  the  values  given  above  with  the  excitation  applied  on  column  1  at  a 
distance  or  0.95  from  the  ground.  Rj-  is  seen  to  have  peaks  near  the  natural  medal 
frequencies  of  the  structure.  Thus,  it  Is  concluded  that  structures  should  be  tested 
near  resonance. 


N=3,  M=4;  TRANSMISSiBIUTY  STATION  ON  16  AT  .61 
EXCJT.  ON  1  AT  .95;  CRACK  ON  17  AT  .1  WITH  TH£TA=.024 


Figure  19.  Relative  transmissibility  change  vs  frequency. 


0 

The  amplitudes  of  dimensionless  transverse  deflections  ,\y  |,  along  girder  15 
•with  and  without  a  crack  on  girder  17  are  shown  in  Figure  20.  "where  dimensionless 
deflection  for  the  uniform  frame  is  defined  by 


SISP  _  Lo'J 

V  "  LzFm  V  ' 


(4.7) 


with  y  given  by  Equation  (2.63).  A  computer  program  is  written  to  compute  the 
dimensionless  deflections  on  a  frame  with  uniform  material  properties,  but  ciflerent 
beam  lengths  in  general.  The  listing  of  the  program  is  given  in  Appendix 


For  an  undamp  eu  (i.e.,  d=0)  structure,  the  curves  in  Figure  20  wouia  be  ceiled 
modal  shapes.  The  pseudo-node  point,  PNP,  of  the  girder  is  observed  to  have  shifted 
after  the  damage.  The  curves  in  Figure  20  also  represent  the  transmissibilities  (Ea. 
4.3)  with  and  without  the  crack  as  a  function  of  the  location  of  the  response  station  on 
girder  IS.  These  two  curves  yield  the  Rt  curve  in  Figure2i  which  indicates  that  in¬ 
takes  its  minimum  and  maximum  values  at  the  RXPs  of  the  post-damage  and  pre¬ 
damage  structure,  respectively.  Figure  22  shows  the  variation  of  R?  as  a  function  of 


N=3,  M=4;  EXC.  ON  It  CRACK  ON  17  AT  TKE7A=.C24 
AMFUTUOE  OF  TRANSVERSE  DEFLECTIONS  ON  16;  aETA=.S6 


\  LaS* 


LOCATION  ALONG  TKE  CIRCES  IS 


w/  o»*cx 


Figure20.  Transverse  deflections  on  girder  IS. 


N  =  3,  M  =  4,  EXCITATION  ON  1  AT  0.25,  UNDAMPED  BETA  =  .SS 
CRACK  ON  17  AT  0.1  WITH  THETA  =  0.024 


4  4.1  4J  4.J  Xi  4J  4.4  4.7  4J  4.4 

LOCATION  OF  TRANSWISSiEiUTY  STATION  ON  IS 


Figure  21.  Relative  transmissibUity  change  vs  response  location. 


the  location  of  the  resnonse  station  on  15  -rhea  there  is  a  crack  cn  column  4. 


ANSMISSIBIllIY  CHANGE  ^  RELATIVE TRANSUISSIDtUTY  CHANGE 


N=3,  W=4,  EXCITATION  ON  1  AT  Q.S5,UNDAMFED  EETA=0.S5 
CRACK  ON  A  AT  0.7  WITH  THETA=0.024 


gure  22.  Relative  transmissibility  change  vs  response  location. 


N=3,  M=4,  EXC.  ON  1  AT  0.95,  UNDAMPED  SE7A=0.25 
CRACK  ON  17  AT  0.1  WITH  TKE7A=0.C24 


0.01- 

0.00  {  -  f«>  |  ■  ■— .pf— —  |  1  1  1  1  '  T  —  — 1  *  .  > 

9  0.1  OJ  OJ  Ou  9-9  9.4  OJ  OU  OJ 

LOCATION  OF  RESPONSE  STATION  CM  COLUMN  4 

Figure  23.  Relative  transmissibility  change  vs  response  location. 
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At  lower  natural  frequencies,  the  columns  of  the  structure  do  not  have  pseudo- 
ncde  points.  Figure  23  shows  that,  at  the  fundamental  frequency  (Le.,  £*=0.36),  R? 
remains  nearly  constant  when  the  location  of  the  response  station  is  varied  along  the 
column  4. 

In  summary,  then,  Rf  can  be  maximized  when  the  excitation  is  on  a  column  near 
a  joint  and  at  a  resonant  frequency  and  when  the  response  station(s)  is(are)  near  the 
PNP(s)  of  a  girder(s). 

Table  1  lists  the  PN?  locations  at  the  fundamental  frequency,  on  some  of  the 
girders  of  the  frame  in  Figure  17  In  the  table,  denotes  the  excitation  location  on 
the  respective  column.  All  of  the  distances  are  normalized  and  measured  relative  to 
the  lower  (left)  end  of  the  respective  column  (girder).  As  the  table  indicates,  a  PJJP 
location  is  highly  insensitive  to  the  excitation  location  and  to  the  value  of  the  damping 
factor. 

4.3.  The  Number  of  Response  Stations 

It  is  desirable,  from  a  practical  point  of  view,  to  accomplish  damage  diagnosis 
with  a  minimum  number  of  transducers.  This  number  is  related  to  the  objectives  of  a 
diagnosis  process.  In  theory,  a  transducer  may  be  installed  on  every  girder  and  the 
crack  can  be  isolated  to  within  a  beam  element.  In  the  Light  of  the  high  sensitivity  of 
the  transmissibility  changes  to  the  response  location  in  the  vicinity  of  the  PNPs  and 
the  uncertainties  involved  in  practice,  however,  the  objective  in  this  study  is  to  isolate 
a  crack  to  within  one  or  two  cells  of  a  frame  structure.  A  local  method,  such  as  an 
ultrasonic  technique,  may  then  be  utilized  to  determine  the  exact  location  of  the 
crack.  With  this  objective  in  mind,  the  six  response  stations  on  the  frame  of  Figure  17 
were  chosen.  These  stations  are  on  girders  16,  19,  21,  22,  24,  and  27.  Numerical 
results  indicate  that,  in  general,  response  at  a  measuring  station  is  most  significantly 
ajfected  when  the  crack  is  at  a  frame  cell  neighboring  the  station.  In  Figure  17 
station  26,  for  example,  has  four  neighboring  ceils  while  station  21  has  six.  The 
criterion  used  here  for  selecting  the  stations  is  that  each  station  has  at  most  two 
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neighboring  cells  in  common  with  those  of  another  station.  This  implies  that  ever7 
other  girder  on  a  floor  level  has  a  station  on  it  and  that  every  other  girder  is  skipped 
as  one  moves  vertically  up  the  frame.  The  number  of  transducers  (response  stations) 

required  is,  then,  when  nm  is  even  and  when  nm  is  odd,  where  n  is  the 

■  number  of  stories  and  m  is  the  number  of  spans.  The  selection  of  stations  in  Figure  17 
violates  the  above  criterion,  even  though  the  number  of  stations  is  still  six  Most  of 
the  numerical  computations  were  carried  out  with  the  configuration  of  Figure  17 
therefore,  the  development  of  the  method  is  based  on  that  configuration. 


TABLE  1 

PNFs  of  Some  Girders  (No  crack) 

P*  = 

.96; 

=  .35  unless  noted  otherwise  j 

Girder  16 

5 

Column  with  the  excitation 

1 

7 

11  ( |/=.55) 

15 

.01 

.6062 

.6084 

1 

.6033 

.30 

.5994 

•  •• 

£H£| 

.6079 

.30* 

.5932 

•  •• 

■1 

.5016 

Girder  24 

5 

Column  with  the  excitation 

(1 

3 

6 

7 

.01  ! 

.6133 

.6140 

.5249 

.6137 

Other  GLrders(excitation  on  7) 

<5 

Girder  No. 

19 

21  ' 

22 

27 

.01 

.3916 

.4895 

.5105 

.3513  | 

*  (With  a  crack  of  9=.G24  on  girder  17  at  a  distance  of  .1) 
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4.4.  Selection  of  the  Column  for  Excitation 

Figures  24-27  show  the  relative  transaiissibility  change  curves  resulting  from, 
four  different  excitation  locations  when  the  other  parameter  values  are  kept  the 
same.  Trans missibility  change  at  station  27  for  a  crack  on  girder  17  is  small  in  all  the 
cases;  therefore,  the  corresponding  curve  is  omitted.  The  figures  indicate  that  the 
trends  of  Rr  values,  for  the  given  response  stations,  are  highly  insensitive  to  the 
excitation  location.  A  trend  of  relative  transmissibility  changes  is  defined  as  the 
order,  with  respect  to  the  magnitudes,  of  the  relative  transmissibility  changes  for  a 
given  crack  configuration.  The  word  "trend",  however,  will  be  used  in  a  more 
qualitative  manner  rather  than  with  an  implication  of  a  precise  ordering  of. numbers. 
For  a  crack  cn  a  certain  beam  element,  for  instance,  two  response  stations  may 
exhibit  the  largest  transmissibility  changes  without  one  response  clearly 
overwhelming  the  other,  in  which  case  the  trend  will  be  described  as  sucn.  From 
Figures  2S-31,  then,  the  trend  of  f?r  values  for  a  crack  on  girder  17  is  such  that 


N=3,  M=4,  CRACK  ON  17,  TH£TA=.Q24,  EXC.  ON  1  AT  .95,  BETA=.96 
X16  =  .618,  X21  =  .475,  X22  =  .525,  X24  =  .649,  X27  =  .351 


Legend 

IT  » l« 

IT  »  21 
IT  »  22 
jT-24  _ 
IT»  17 


Figure  24.  Relative  transmissibility  change  vs  crack  iocation(excitation  cn  l). 


RELATIVE  TRANSMlSSIOILflY  CHAIIGE 


Figure  26. Relative  transmissibility  change  vs  crack  location(escitation  on  7) 
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N=3,  M=4,  CRACK  ON  17,  THETA=.024,  EXC.  ON  13  AT  .95,  BETA=.9S  . 
X16  =  .618,  X21  =  .475,  X22  =  .525,  X24  =  .649,  X27  =  .351 


Legend 

IT  «  n 

it  » at  ■ 

ir»a 

IT 


Figure  27.  Relative  transmissibility  change  vs  crack  location( excitation  on  13). 


station  16  displays  a  much  more  significant  transmissibility  change, than  the  other 
stations  ever  a  large  portion  of  the  range  of  the  crack  location.  However,  when  the 
response  stations  are  very  close  to  the  PN?s,  unlike  in  the  figures,  the  trends  may  be 
more  sensitive  to  the  excitation  location.  It  should  be  noted  that  relative 
transmissibility  changes  and  transmissibilities  themselves  change  as  the  excitation 
location  is  varied  The  interest  of  this  study,  however,  is  in  the  trends  rather  than  the 
absolute  values.  In  what  follows,  the  general  trends  of  Rj  values  will  be  studied  by 
exciting  the  structure  on  column  7  at  a  distance  of  0.95  from  the  ground  and  by 
studying  only  the  left  half  of  the  frame  and  utilizing  the  symmetry. 

4.5.  Selection  of  Response  Locations  Relative  to  the  FNPs 

Figures  21  and  22  demonstrate  that  relative  transmissibility  changes  are  highly 
sensitive  to  the  response  location  in  the  vicinity  of  the  PNPs.  In  practice,  however,  it 
may  not  be  possible  to  determine  the  exact  location  of  a  PNP.  On  the  other  hand  it 
may  not  be  desirable  to  install  a  transducer  at  a  PITP  either  since  the  acceleration  or 
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deflection  there  may  be  too  small  to  be  measurable.  Transducers  then  must  be 
installed' at  some  distance  from  the  PNP s  and  yet  must  be  close  enough  to  the  PNPs 
to  yield  large  R?  values.  For  some  crack  configurations,  installing  the  transducers  to 
the  right  of  the  pre-damage  PNP  yields  better  results  (Fig.  21). .  whereas  for  others 
installing  them  to  the  left  is  better  (Fig.  22) .  This  is  due  to  the  fact  that  PNPs  may 
shift  to  the  right  or  left  after  the  damage.  The  situation  may  be  improved  by 
installing  two  transducers  in  the  vicinity  of  a  PNP  in  close  proximity  to  each  other. 
Such  a  measure  Is  likely  to  increase  the  probability  of  obtaining  a  large  Rf  value  for  a 
given  crack. 

4.5.1.  The  Amount  of  Shift  of  the  Response  Stations  Relative  to  the  PNPs  • 

Table  2  lists  the  amplitudes  of  dimensionless  deflections  of  selected  girders  on 

the  undamaged  frame  in  the  vicinity  of  PNPs  as  well  as  the  lateral  deflections  of  the 

floor  levels  for  the  frame  of  Figure  17  with  the  excitation  on  column  7.  The  undamped 

P  value,  Pi,  is  0.96  (Le.,  at  the  fundamental  frequency).  The  amplitude  of  complex 

# 

deflection  can  be  shewn  to  he 

IV*l  =  TTr-a  +  5*)‘'»!yi  (4.8) 

Pd *  m 

from  Equations  (4.7)  and  (2.5),  where  |y|  is  the  physical  deflection.  An  empirical 

20t7 

formula  for  the  fundamental  frequency  of  multistorv  buildings  is  y,  =  — —  in  rad/s 
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[33], which  is  the  Uniform  Building  Code  formula,  where  n  is  the  number  of  stories.  For 
a  three  story  building,  then,  ~  21  rad/s.  If  it  is- assumed  that  the  amplitude  of  the 
force,  Fm,  is  approximately  equal  to  the  weight,  Lpg ,  of  one  beam  element,  where  g  is 
the  gravitational  constant,  equation  (4.3)  can  be  written  as 

|y*|  »L2|y|  ,  |y|  in  inches.  (4.9) 

y  Eddy-current  displacement  transducers  are  capable  of  measuring  deflections  less 

than  10-3  inch  (0.025  mm).  If  this  value  is  accepted  as  the  base,  |y  #|  must  be  of  the 

order  10“3  in  order  for  |y  |  to  be  measurable.  From  Table  4,  then,  the  locations  of  the 

response  stations  are  determined  to  be  at  the  points  away  from  the  PNPs.  where  j  y  *| 
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TABLE  2 

Deflections  at  Selected  Locations 
Excitation  on  column  7  at  .95 
ft*  =  .96,  6  =  .01 


Side  sway  of  the  Floors 


Floor 

\vm\ 

1 

0.438 

2 

0.909 

3 

1.157 

at  the  excitation 

1  0.422 

In  the  vicinity  of  the  girder  P>J?s 

Girder 

Location 

ly'lxlO3 

.588 

3.04 

.5S8 

1.53 

16 

.60  S4(PN?) 

0.008 

.618 

1.40 

1 

.623 

2.33 

■  ! 

.475 

1.32 

21  I 

.4895(PN?) 

0.014 

i 

.505 

1.41 

.539 

1.14 

24 

.61S7(PNP) 

0.003 

.649 

1.09 

reaches  a  value  of  at  least  10-3.  This  translates  into  a  shift,  away  from  the  respective 
PNF,  of  0.01  for  girders  IS  and  19,  0.015  for  21  and  22,  and  0.03  for  24  and  27.  The 
response  locations  used  in  Figures  24-27  were  based  on  these  shift  values.  As  the 
response  stations  are  moved  further  away  from  the  PNPs,  Rf  values  get  smaller  and 
become  comparable  in  value  as  a  result  of  which  the  trends  are  gradually  lost.  For 
example,  when  the  locations  of  the  response  stations  are 
£ia  =  0.35,  £ig  =  0.15,  £21  =  0.20,  £22  =  0.50,  £24  =  0.35,  and  £27  =  0.15,  all  of  which 
are  measured  from  the  left  end  of  the  respective  girder,  all  of  the  relative 
transmissibiiity  changes  for  a  given  crack  are  either  small  (less  than  10%)  or  close  in 


N=3,  M=4,  CRACK  ON  17,  TH£TA=.024,  EXC.  ON  7  AT  .95,  8ETA=.96 
X1&=.653,  X21=.475,  X22=.525,  X24=.S49,  X27=.351 


Figure 28.  Relative  transmissibility  change  vs  crack  Location. 


value.  On  the  other  hand,  moving  the  response  stations  closer  to  the  PNPs  increases 
the  possibility  of  diagnosis  as  the  trends  become  more  definitive. 

It  is  desirable  to  know  how  far  away  from  the  ?NPs  the  trends  still  persist.  Figure 
28  illustrates  that  at  a  distance  of  0.05  from  the  ?N?  of  girder  16,  the  trend  for  a 
crack  on  girder  17  is  still  clear,  when  the  crack  is  on  the  left  half  of  the  frame  and 
when  the  other  response  stations  are  kept  at  the  same  locations  as  in  Figures  24-27. 
Further  away  from  the  PNP,  however,  the  trend  starts  becoming  Less  cLear-cut. 

4.5.2.  The  Direction  of  Shift  of  the  Response  Stations  Relative  to  the  FNPs 

♦ 

Another  question  related  to  the  optimum  Location  of  the  response  stations 

relative  to  the  PNPs  is  that  of  which  direction  to  shift  the  response  stations.  Tables  3 
v 

and  ‘4  list  the  relative  transmissibility  changes  at  the  sis  response  stations  as  a 
function  of  the  crack  location.  In  Table  3,  the  response  stations  on  girders  16  and  19 
are  shifted  from  the  corresponding  PNPs  by  a  distance  of  0.01  away  from  the  center, 
while  those  on  21  and  22  are  shifted  by  a- distance  of  0.015  toward  the  center  and  the 
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TABLE  3 


Relative  Transmissibility  Changes  vs  Crack  Location 


Excitation  on  column  7  at  a  distance  of  .95,  ^  =  .96,  6  =  .01,  0  =  .024. 


Response  Locations  : 

fa  =  .593  £21  =  .505  $24  =  .589 

$19  -.402  $22  =  .495  $r?  =  -411 


($e  -  crack  Location,  RS  :  response  station) 


1)  Crack  on  column  1 

RS  16 

RS  19 

RS  21 

RS  22 

RS  24 

RS  27 

0.001 

-0.440 

0.239 

0.257 

0.224 

0.179 

■-  0.235 

0.100 

0.163 

0.171 

0.153 

0.123 

0.160  - 

0.200 

-0.067 

0.103 

j.  104 

0.096 

0.037 

0.101  ■ 

0.300 

0.027 

0.059 

0.05S 

0.054 

0.054 

0.057 

0.400 

0.071 

0.023 

0.023 

0.025 

0.030 

0.02S 

0.500 

0.067 

0.008 

0.004 

0.007 

0.013 

0.003 

0.600 

0.013 

0.000 

-0.001 

0.000 

0.002 

0.000 

0.700 

-0.075 

0.003 

0.007 

0.003 

-0.003 

0.003 

O.SCO 

-0.214 

0.016 

0.023 

0.016 

-O.0C2 

0.017 

C.900 

-0.400 

0.040 

0.051 

0.040 

0.005 

0.041 

0.999 

-0.635 

0.074 

0.107 

0.073 

0.017 

0.076 

2)  Crack  on  column  2 

.  RS  16 

RS  19 

ES 

RS  22 

RS  24 

RS  27 

0.001 

-0.4O4 

0.043 

0.023 

C.049 

0.118 

0.038 

0.1C0 

-0.301 

0.024 

0.021 

0.029 

0.058 

0.022 

0.200 

-0.177 

0.011 

0.015 

0.015 

0.019 

0.010 

0.300 

-0.031 

0.003 

0.008 

0.005 

-0.001 

0.003 

0.400 

-0.013 

0.000 

0.002 

0.000 

-0.002 

0.000 

0.500 

0.029 

0.002 

-0.005 

0.000 

0.013 

0.001 

0.600 

0.046 

0.007 

-0.012 

0.004 

0.044 

0.005 

0.700 

0.037 

0.017 

-0.019 

0.091 

0.012- 

O.SOO 

0.004 

0.030 

-0.027 

0.021 

0.153 

0.022 

0.900 

-0.052 

0.045 

-0.036 

0.035 

0.223 

0.035 

0.999 

-0.132 

0.064 

-0.045 

0.050 

0.315 

0.049 
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TABLE  3 
(continued) 


3)  Crack  on  column  4 


0.001 

0.100 

0.200 

0.300 

0.400 

0.500 

0.600 

0.700 

0.300 

Q.SOO 

0.999 


RS  IS 


5)  Crack  on  girder  16 


RS  19 


RS  21 


RS  22 


.RS  24 


0.2S9 

0.170 

0.098 

0.048 

0.018 

0.002 

0.0C2 

0.013 

0.051 

0.102 

0.174 


0.299 

0.174 

0.037 

0.032 

0.002 

-0.004 

0.013 

0.054 


0.119 

0.214 

0.344 


0.233 

0.174 

0.096 

0.043 

0.012 

-0.000 

0.006 

0.030 

0.073 

0.139 

0.232 


0.  u46 
0.223 1  1 

0.133 


0.069 

0.027 

0.004 

-0.001 

0.011 

0.041 

0.092 

0.167 


RS  27 


0.253 
0.1S2 
0.092 
0.044 
0.015 
0.001 
0.003 
0.020 
0.C53 
0.1 05 
0.177 


t. 

RS  15  • 

RS  19 

RS  21 

'  RS  22 

RS  24 

RS  27 

0.001 

1.203 

0.349 

0.172 

0.269 

0.370 

0.323 

0.100 

0.680 

0.221 

a  123 

0.175 

0.231 

0.209 

0.200 

0.321 

0.123 

0.083 

0.105 

0.131 

0.121 

0.300 

0.09S 

0.063 

0.052 

0.056 

0.062 

0.060 

0.40C 

-0.013 

0.022 

0.023 

0.023 

0.020 

0.022 

0.500 

-0.035 

0.003 

0.009 

0.005 

0.001 

0.003 

0.600 

0.039 

0.002 

-0.005 

-0.000 

0.003 

0.0C2 

0.700 

0.207 

0.C20 

-0.014 

0.007 

0.027 

0.018 

0.300 

0.475 

0.058 

-0.019 

0.023 

0.072 

0.052 

0.900 

0.858 

0.117 

-0.021 

0.063 

0.141 

0.106 

0.999 

1.372 

0.200 

-0.020 

0.113 

0.238 

••  0.131 

4)  Crack  on  column  5 

_ _ _ . 

RS  IS 

RS  19 

RS  21 

RS  22 

RS  24 

RS  27 

0.001 

1.205 

0.171 

0.560 

0.205 

0.015 

0.133  1 

0.100 

0.773 

0.103 

0.237 

0.115 

0.025 

0.030 

0.200 

0.450 

0.054 

0.107 

0.053 

0.025 

0.043 

0.3C0 

0-219 

0.022 

0.008 

0.015 

0.022 

0.013 

0.400 

0.067 

0.004 

■0.021 

-0.001 

0.011 

0.004 

0.500 

- 

0.016 

0.000 

0.017  [ 

0.003 

-0.005 

-0.000 

0.500 

• 

0.034 

0.003 

0.118 

0.Q2S 

-0.023 

0.0C5 

0.700 

0.012 

0.029 

0.235 

0.067 

-0.058 

0.019 

0.300 

0.124 

0.061 

0.521 

0.123 

-0.097 

0.C42 

0.900 

0.306 

0.107 

0.335 

0.210 

-0.146 

0.074 

0.999 

0.563 

0.165 

1.231 

0.315 

-0.207 

0.115 
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!  TABLE  3  | 

(continued)  J 

S)  Crack  on  girder  17  J 

& 

RS  15 

RS  19 

RS  21 

.  RS  22 

RS  24 

RS  27 

0.001 

-0.3411 

0.130 

0.215 

0.098 

0.085 

0.135 

0.100 

-0.569 

• 

0.112 

0.146 

0.069 

0.050 

0.0S3 

0.200 

-0.354 

0.056 

0.091 

0.045 

0.025 

0.045 

0.300 

-0.190 

0.01S 

0.049 

0.025 

0.003 

0.013 

0.400 

-0.073 

0.001 

0.019 

0.010 

0.000 

0.004 

0.500 

0.003 

0.000 

-0.001 

-0.000 

O.CCO 

0.000 

Q.S00 

0.03S 

0.017 

-0.010 

-Q.007 

0.003 

C.007 

0.700 

0.032 

0.052 

-0.003 

-0.010 

0.025 

0.02S 

0.300 

-0.015 

0.107 

0.0C5 

-0.009 

0.C50 

0.056 

0.S00 

-0.108 

0.132 

0.029 

-0.004 

0.084 

0.093 

0.999 

-0.252 

0.230 

0.065 

0.C04 

0.129 

■.  0.155  j 

ones  on  24  and  27  by  a  distance  cf  0.03  away  from,  the  center.  In  Table  4,  the  shifts 
are  in  the  opposite  directions  by  the  same  amounts.  Table  3  shows  only  the  cases 
•when  the  crack  is  on  a  beam  element  in  the  neighborhood  or  response  station  25,  that 
is,  on  1,  2,  4,  5,  IS,  or  17.  A  negative  transmissibuity  change  which  is  large  in 
magnitude  is  considered  significant.  From  Table  3,  then,  the  changes  at  station  15 
are  the  most  significant  among  all  the  stations,  when  the  crack  is  on  an  element 
neighboring  station  16.  It  is  preferable,  however,  to  choose  a  response  station  at  a 
location  where  tr arsmis sib ility  (or  deflection)  increases  after  the  damage,  which 
increases  the  measurability  of  post-damage  deflections.  When  Table  3  is  compared 
with  the  corresponding  segments  of  Table  4,  it  is  observed  that  it  is  more 
advantageous  to  shift  response  station  18  to  the  right  for  a  crack  on  elements  1,  2,  15, 
and  17,  and  to  the  left  for  a  crack  cn  4  and  5.  Figure  20  clearly  illustrates  the  case  for 
a  crack  on  girder  17.  In  the  analysis  which  follows,  response  station  is  takan  to  be 
located  to  the  right  of  the  respective  FNIP  since  this  is  more  advantageous  for  more 
elements.  As  mentioned  earlier,  if  one  transducer  is  mounted  cn  each  side  of  a  ?NF, 
the  prospects  of  diagnosis  are  significantly  improved.  Similar  computafons  and 
comparisons  were  carried  cut  for  crack  locations  on  other  beam  elements.  The 
results  are  summarized  in  Table  5  which  shows  the  preferred  directions  of  shift 
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TABLE  4 

Relative  Transmissibility  Changes  vs  Crack  Location 


Excitation  on  column  7  at  a  distance  of  .95,  /?*  ~  .96,  d  =  .01,  0  =  .024. 


Response  Locations : 


|ia  =  .615  fet  =  .475  -  .649 

£19  =  .332  $22  =  -525  $zr  ~  .351 


(?«  =  crack  location,  RS  :  response  station) 


0.C01 

I  0.100 
j  0.200 
!  Q.300 
10.400 
1 0.500 
Q.SOO 
0.700 
0.300 
0.300 
0.S99 


RS  16 


RS  IS 


l)  Crack  on  column  1 


RS  19 

RS  21 

.  RS  22 

RS  24 

RS  27 

0.220 

0.196 

0.233 

0.277 

0.225 

0.150 

0.139 

0.160 

0.183 

0.154 

0.095 

0.093  • 

0.102 

0.111 

0.098 

0.C54 

0.057 

0.059 

■  0.059 

0.056 

C.025 

0.031 

0.023 

0.023 

0.027 

0.C08 

0.012 

0.009 

0.004 

0.008 

0.000 

0.001 

0.001 

-0.001 

o.cco 

0.003 

-0.002 

0.002 

0.009 

0.002 

0.015 

0.001 

0.014 

0.032 

0.014 

0.037 

0.011 

0.035 

0.069 

0.036 

0.063 

0.027 

0.066 

0.121 

0.067 

2)  Crack  on  column  2 


RS  19 


RS  21 


RS  22 


0.051 

0.030 

0.014 

0.005 

0.000 

0.001 

0.006 

0.016 

0.029 

0.04S 

0.065 


-0.001 

0.008 

0.027 

0.056 

0.093 

0.133 

0.191 


0.051 

0.029 

0.013 

0.003 

-0.000 

0.002 

0.010 

0.022 

0.038 

0.059 

0.082 


-0.015 

-0.000 

0.008 

0.008 

0.002 

-0.010 

-0.029 

-0.054 

-0.086 

-0.125 

-0.169 


RS  27 


0.055 

0.032 

0.015 

0.005 

0.000 

0.001 

0.007 

0.017 

0.031 

0.049 

0.071 


84 


0.001 

0.100 

0.200 

0.300 

0.400 

0.500 

0.600 

0.700 

0.300 

0.900 

0.999 


-0.653 
-0.321 
-0. 104 
0.020 
0.067 
0.046 
-0.041 
-0. 194 
-0.422 
-0.736 
-0.346 


TABLE  4 
(continued) 


3)  Crack  on  column  3 


tc 

RS  IS 

RS  19 

RS  21 

'  RS  22 

0.001 

-0.031 

0.007 

0.043 

0.012 

0.100 

-0.019’ 

0.003 

0.024 

0.008 

0.200 

-0.008 

0.001 

0.010 

0.002 

0.300 

0.000 

0.000 

-0.000 

-0.000 

0.400 

0.006 

0.001 

-0.006 

-0.000 

0.500 

0.011 

0.002 

-0.008 

0.001 

0.600 

0.013 

0.005 

-0.007 

0.003 

0.700 

0.013 

0.008 

-0.004 

0.005 

0.900 

0.011 

0.011 

0.001 

0.008 

0.900 

0.008 

0.014 

0.008 

0.012 

0.999 

0.003 

0.017  • 

0.016 

0.015 

4)  Crack  on  column  4 

fc 

RS  16 

RS  19 

RS  21 

RS  22 

0.2S5 

0.134 

0.1G9 

0.057 

0.022 

0.0C4 

0.000 

0.011 

0.038 

0.077 

0.135 


0477 

0.290 

0.157 

0.068 

0.017 

-0.003 

0.007 

0.C43 

0.120 

0.227 

0.373 


P367 

10.231 

0.132 

0.064 

0.022 

0.002 

0.003 

0.024 

0.066 

0.132 

0.223 


RS  24 


RS  24 


0.257 

0.176 

0.107 

0.058 

0.025 

0.005 

-0.001 

0.005 

0.023 

0.C54 

0.1C0 


RS  27 


0.006 
0.0C2 
0.000 
0.000 
0.001 
0.003 
0.005 
0.010 
0.013 
0.016 
0.0 


RS  27 


0.308 

0.1S9 

0.116 

0.053 

0.023 

0.004 

0.001 

0.014 

0.043 

0.Q30 

0.156 


0.001 

0.100 

0.200 

0.300 

0.400 

0.500 

0.600 

0.700 

0.800 

0.900 

0.999 


RS  16 


5)  Crack  on  column  5 


RS  19 


0.147 

0.086 

0.043 

0.016 

0.002 

0.001 

0.011 

0.032 

0.066 

0.111 

0.168 


RS  21 


RS  22 


0.129 

0.0S4 

0.050 

0.025 

0.008 

-0.002 

-0.0C6 

-0.005 

0.001 

0.012 

0.02S 


RS  24 


RS  27 


0.159 

0.110 

C.056 

0.021 

0.003 

0.001 

0.013 

0.040 

0.081 

0.138 

0.211 


TABLE  4 
(continued) 


6)  Crack  on  column  5 

£e 

RS  18 

RS  19 

RS  21 

'  RS  22 

RS  24 

RS  27 

0.001 

0.106 

0.036 

-0.424 

-0.017 

1 

0.021 

0.100 

0.064 

0.020 

-0.274 

-0.015 

S 

0.012 

0.200 

0.034 

0.009 

-0.155 

-0.011 

-0.036 

0.006 

0.300 

0.012 

■  0.002 

-0.064 

-0.006 

0.010 

0.002 

0.400 

0.000 

0.000 

-0.001 

-0.000 

0.000 

0.500 

-0.004 

0.001 

0.036 

0.007 

-0.060 

-0.001 

0.600 

-0.001 

0.005 

0.048 

0.014 

-0.163 

-0.000 

0.700 

0.008 

0.012 

0.038 

0.022 

-0.320 

0.001 

0.800 

0.023 

0.021 

0.029 

-0.510 

0.0C2 

0.900 

0.043 

0.031 

-0.051 

0.035 

-0.735 

0.004 

0.999 

0.067 

0.043 

-0.122 

0.041 

-0.981 

’•  0.005 

7)  Crack  on  column  7 


RS  16 

RS  19 

RS  21 

RS  22 

RS  24 

RS  27 

0.001 

0.414 

0.414 

0.192 

0.192 

0.100 

0.254 

0.254 

0.136 

0.135 

0.200 

0.139 

0.139 

0.092 

0.092 

0.300 

0.062 

0.062 

0.057 

0.057 

0.400 

0.015 

0.018 

0.030 

0.030 

0.500 

- 

0.002 

-0.C02 

0.010 

0.010 

0.600 

0.CC5 

•  -0.004 

-0.0C4 

0.700 

0.039 

0.039 

-0.013 

-0.013 

0.800 

0.099 

-0.017 

-0.017 

Q.SOQ 

0.188 

-0.016 

-0.016 

0.999 

0.301 

0.301 

-0.012 

-0.012 

0.317 

0.317 

0.2C3 

0.203 

0.1191 

0.119 

0.060 

0.022 

0.003 

0.001 

0.016 

0.048 

O.0S9 

0.164 


0.060 

0.022 

0.003 

0.C01 

0.016 

0.048 

0.099 

0.154 


A 

Sc 


S)  Crack  on  column  3 


RS  15 


RS  19 


RS  21 


RS  22  I  RS  24 


RS  27 


0.001 

0.279 

0.279 

0.525 

! 

0.525 

0.132 

0.132  ! 

0.100 

0.175 

0.175 

0.289 

0.269 

0.083 

0.083 

0.200 

0.098 

0.098 

0.099 

0.099 

0.046 

0.046 

0.300 

0.044 

0.044 

0.021 

0.021 

0.400 

0.012 

0.012 

-0.020 

-0.020 

0.0C6 

0.0C6 

0.500 

-0.002 

-Q.C02 

0.019 

0.019 

-0.001 

-0.001 

0.600 

0.001 

0.001 

0.119 

0.119 

0.000 

•  0.000 

0.700 

0.021 

0.021 

0.232 

0.282 

0.009 

0.009 

0.800 

0.059 

0.059 

0.512 

0.512 

0.025 

0.025 

0.900 

0.114 

0.114 

0.813 

0.313 

0.049 

0.049 

0.999 

0.183 

0.188 

1.139 

1.189 

0.079 

0.079 

86 


TABLE  4 
(continued) 


9)  Crack  on  column  9 


RS  21 


RS  22 


RS.  24 


RS  27 


0.423 

0.2SS 

0.146 


0.03S 


0.104 

0.1S6 


0.104 

0.196 


0.065 

0.029 

0.007 

-0.002 

0.000 

0.014 

0.037 

0.069 


0.151 

0.1S9 


0.065 

0.029 

0.007 

-0.002 

0.000 

0.014 

0.037 

0.068 


0.107 

0.151 

0.199 


RS  IS 


10)  Crack  on  girder  16 


RS  19 


RS  21 


RS  22 


RS  24 


RS  27 


0.001 

0.100 

0.200 

0.300 

0.400 

0.500 

0.600 

0.700 

0.3C0 

0.900 

0.S99 


0.279 

0.173 

0.097 

0.045 

0.014 

0.001 

0.005 

0.025 

0.066 

0.127 

0.212 


0.255 

0.174 

0.112 

0.C64 

0.030 

0.006 

-0.008 

-0.012 

-0.008 

0.007 

0.032 


0.273 

0.175 

0.102 

0.051 

0.019 

0.002 

0.001 

0.014 

0.044 

0.091 

0.159 


0.211 

0.126 

0.065 

0.025 

0.004 

-0.CQ2 

0.008 

0.033 

0.075 

0.137 

0.222 


0.299 

0.190 

0.101 

0.048 

0.015 

0.001 

0.004 

0.024 

0.063 

0.123 

0.207 


11)  Crack  on  girder  17 


RS  IS 


RS  19 


RS  21 


RS  22 


RS  24 

RS  27  ! 

0.190 

0.140 

0.121 

0.088 

0.069 

0.049 

0.032 

0.022 

0.009 

0.006 

-0.000 

-0.000 

0.003 

0.004 

0.020 

0.019 

0.051 

0.045 

0.097 

0.093 

0.160 

0.134 

0.001 

0.100 

0.200 

0.300 

0.400 

0.500 

0.500 

0.700 

0.200 

0.900 

0.939 


TABLE  4 
(continued) 


12)  Crack  on  girder  20 


$c 

RS  16 

RS  19 

0.001 

0.239 

0.123 

0.100 

0.168 

0.080 

0.200 

0.110 

0.047 

0.300 

0.064 

0.023 

0.400 

0.C29 

0.007 

0.500 

0.005 

0.000 

0.600 

-0.009 

0.002 

0.700 

-0.013 

0.011 

0.500 

-0.C07 

0.029 

0.900 

0.009 

0.055 

0.999 

0.036 

0.090 

RS  21 


RS  22 


0.170 

0.102 

0.053 

0.019 

0.001 

-0.002 

0.010 

0.035 

0.076 

0.133 

0.207 


£e 

RS  16 

RS  19  | 

0.001 

-0.037 

0.084 

0.100 

-0.030 

0.C52 

0.2C0 

-0.023 

•  0.027 

0.3G0 

-0.015 

0.011  1 

Q.4C0 

-0.008 

0.002  1 

0.500 

0.000 

o.ooo. 

0.600 

0.008 

0.006 

0.7C0 

0.017 

0.019 

0.300 

0.025 

0.040 

0.3CQ 

0.035 

0.069 

0.999 

0.C44 

0.1C6 

13)  Crack  on  girder  21 


RS  21 


14)  Crack  on  gird 


RS  18 


RS  19 


RS  21 


RS  22 


RS  24 


£348 

0.250 

0.168 

0.102 

0.049 

0.009 

-0.019 

-0.035 

-0.039 

-0.031 

-0.012 


RS  27 


0.123 

0.081 

0.048 

0.024 

0.008 

0.001 

0.001 

0.009 

0.025 

0.049 

0.082 


RS  24 


RS  27 


0.003 

0.001 

-0.001 

-0.001 

-0.001 

-0.001 

0.001 

0.003 

0.007 

0.011 

0.015 


0.003 

0.001 

0.000 

.0.000 

0.002 

0.004 

0.008 

0.013 


0.016 

0.015 

0.014 

0.011 

0.007 

0.002 

-0.004 

-0.011 

-0.019 

-0.023 

-0.038 


0.015 
0.011 
0.008 
0.005 
0.003 
0.001 
-0.001 
-0.002 
-0.0  02 
-0.003 
-0.002 


0.532 

0.223 


-0.031 

-0.166 

-0.131 

-0.075 

0.150 

0.325 

0.150 

-0.311 

-0.957 


0.019 

0.012 

0.006 

0.002 

-0.000 

-0.000 

0.001 

0.005 

0.011 

0.013 

0.027 


TABLE  4 
(continued) 


15)  Crack  on  girder  25 


RS  16 

RS  19 

RS  21 

•  RS  22 

RS  24 

— 

RS  27 

0.001 

0.017 

0.008 

-0.022 

0.020 

0.454 

-0.020 

0.100 

0.011 

0.005 

-o.oia 

0.010 

0.317 

-0.009 

0.200 

0.007 

0.003 

-0.015 

0.003 

Q.201 

-0.001 

0.300 

•  0.004 

0.002 

-0.010 

-0.001 

0.108 

0.003 

0.400 

0.001 

0.000 

-0.C05 

-0.002 

0.040 

0.C03 

0.500 

-0.000 

-0.000 

0.000 

0.000 

-0.002 

-0.000 

0.600 

-0.000 

0.000 

0.007 

0.006 

-0.008 

0.700 

0.001 

0.001 

0.013 

0.014 

-0.009 

-0.020 

0.800 

0.003 

0.002 

0.020 

0.025 

0.026 

-0.035 

0.900 

0.006 

0.004 

0.027 

0.039 

0.087 

-0.053 

0.999 

0.011 

0.006 

0.034  ' 

0.056 

0.171 

.  -0.075 

relative  to  the  PNFs  (to  the  left  or  to  the  right)  for  all  the  beam  elements  on  the  left 
half  of  the  frame.  The  preferred  direction  of  shift  is  indicated  only  when  the  cracked 
beam  element  is  In  the'  neighborhood  of  the  response  station  in  question.  From  Table 
5 ,  it  is  concluded  that  response  stations  15  and  24  should  be  shifted  to  the  right  and 
station  21  to  the  left  of  the  respective  ?NPs  of  the  undamaged  frame. 


4.6.  Development  of  the  General  Trends 

In  this  section  the  transmissibiiity  method  will  be  developed  with  the  3x4  frame 
of  Figure  17.  The  procedure  used  is  to  vary  the  location  of  a  single  crack  of  known 
sensitivity  number  throughout  the  frame,  on  each  column  and  girder  in  turn,  and  to 
compute  the  relative  transmissibiiity  change  at  each  one  of  the  six  response  stations 
corresponding  to  each  crack  location.  Trends  of  relative  transmissibiiity  changes  as  a 
function  of  the  crack  location  are  then  established,  which  can  serve  as  a  diagnosis 
procedure  in  a  practical  situation  wherein  the  crack  location  is  unknown  and  the 
transmissibiiity  changes  are  known  from  measured  data. 


4.6.1.  Exciting  the  Frame  at  its  Fundamental  Frequency 


Relative  transmissibiiity  changes  computed  for  the  frame  excited  at  the 
fundamental  freo.uency  are  shown  in  Table  4.  Table  4  is  arranged  in  15  groups  in  each 


Direction  of  Shift  of  the  Response  Stations 
(DJD  =  No.  of  the  beam,  element  with  the  crack, 

RS:  Response  Station, 

R:  Shift  the  response  station  to  the  right  of  the  PNP, 
L:  Shift  the  response  station  to  the  left  of  the  PNP.) 


of  which  numerical  data  for  a  crack  on  a  particular  beam  element  is  presented.  Crack 
location  in  each  group- of  data  is  varied  from  the  lower  (left)  end  of  the  element  to  the 
upper  (right)  end  for  a  column  (girder).  Significant  values  are  shown  framed  in 
the  table,  "significance"  being  determined  by  the  following  criteria. 

i)  An  Rr  value  whose  magnitude  is  greater  than  0. 1  is  considered  significant. 

ii)  If  all  or  most  of  the  response  stations  indicate  significant  transmissibility 
changes  for  a  certain  crack  location,  but  some  of  them  incioate  much  larger  changes 
than  others,  then  only  the  larger  R?  values  are  framed.  Examples  can  be  seen  in  dam 
groups  1  and  4. 

iii)  If  several  response  stations  indicate  significant  and  comparable 
transmissibility  changes  for  a  certain  crack,  the  overall  trend  of  transmissibility 
changes  for  that. crack  is  considered  insignificant.  An  example  in  Tabie  4  is  the  case  of 
a  crack  on  column  4  at  a  location  of  0.2  from  the  ground.  All  of  the  R?  values  for  that 
crack  have  magnitudes  between  0.10  and  0.16.  Such  a  crack  may  be  detectable,  but 
is  cot  diagncsabie. 
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A  crack  is  considered  diagnosabie  if  it  can  be  located  within  one  or  two  ceils  of 
the  frame;  that  is,  if  the  trend  of  R7  values  for  that  crack  is  as  described  in  (ii) 
above.  Thus,  the  ranges  of  crack  locations  Q.-0.29  and  0.71-1.  on  column  1.  and  the 
ranges  0.-0.2S  and  0.54-1.  on  column  6,  for  example,  ar?  diagnosabie.  The  diagnosabie 
regions  of  all  the  beam  elements  in  Table  4  and  the  response  stations  exhibiting 
significant  transmissibility  changes  in  each  one  of  those  regions  are  summarized  in 
Table  6.  Also  shown  are  the  percentages  of  the  lengths  of  the  beam  elements  which 
are  diagnosabie.  The  following  conclusions  are  drawn  from  Table  6 . 

i)  The  response  station(s)  closest  to  the  crack  exhibit  the  most  significant 
transmissibility  change(s)  with  the  exception  that  stations  24  and  27  experience 
larger  changes  than  stations  21  and  22  for  a  crack  on  column  7.  However,  since 
stations  16  and  19  exhibit  the  most  significant  changes,  the  diagncsabiiity  of  a  crack 
cn  7  is  unaffected.  For  a  crack  in  region  0. 9 1-0. 92  of  girder  IS,  the  most  significant 
changes,  which  are  about  0.13,  are  experienced  at  stations  19,  24,  and  27,  which  may 
cause  a  misdiagnosis. 

ii)  3eam  elements  with  diagnosabie  regions  larger  than  75%  are  5,  6,  16.  21,  and 
24,  the  last  three  of  which  have  response  stations  on  them.  Both  ends  of  the  columns 
5  and  5,  on  the  other  hand,  adjoin  the  end  of  a  girder  with  a  response  station.  That  is, 
column  5  is  neighbors  with  girders  16  and  21,  and  column  6  is  with  21  and  24. 

iii)  All  of  the  regions  near  the  joints  are  diagnosabie,  with  the  small  exception  of 
the  region  Q.-0.Q2  at  the  left  end  of  girder  21.  In  that  region,  response  station  22 
exhibits  the  most  significant  transmissibility  change  which  may  cause  a  misdiagnosis. 

iv)  The  mid-regions  of  all  the  elements,  with  the  exception  of  those  carrying  a 
response  station,  are  undiagnosable.  At  the  fundamental  frequency,  these  regions  are 
usually  regions  of  low  stress  or  low  resisting  moment.  Slope  discontinuity  in  these 
regions  is,  therefore,  small,  for  a  given  crack  depth  (Section  2.3),  which  reduces  the 
efiect  of  the  crack  on  the  structural  response. 


TABLE  6 
(continued) 


|  DR 

RS 

?D(%) 

.00  -  .14 

21.24 

20 

.14  -  .30 

24  . 

67 

.63  - 1.0 

21 

.00  -  .02 

22 

.02 -.03 

22,21 

.03  -.19 

21 

21 

.22  -  .47 
.53  -  .71 

21 

21 

91 

.71  -  .74 

21,22 

.74 -.31 

22 

.31  -  1.0 

21,22 

.00  -  .15 

24 

24 

.25  -  .47 
.53  -  .31 

24 

24 

75 

.35  -  1.0 

24 

25  | 

.00  -  .31 
.92  -  1.0 

24 

24 

39 

Diagnosafciiity  for  the  whole  frame  =  62. 7  %  .  j 


The  overall  diagncs  ability  of  the  frame  excited  at  its  fundamental  frequency  car 
be  computed  by  adding  all  the  percentages  in  Table  6,  multiplying  je  sum  by  two, 
subtracting  the  percentages  for  7  through  9  once  and  dividing  the  result  by  the  total 
number  of  beam  elements,  27.  This  yields  a  diagncsability  for  the  whole  frame  of 
62.7%.  In  other  words,  the  mid-regions  of  all  the  beam  elements  which  are 
undiagnosabie  make  up  the  37.3%  of  the  total  frame  area. 

4.S.2.  Exciting  the  Frame  at  the  Third  Modal  Frequency 

The  previous  section  indicated  that  mid-regions  of  the  beam  eiements,  which  are 
regions  of  low  resisting  moment  at  the  fundamental  frequency,  are  undiagnosabie.  If 
the  frame  is  excited  at  a  higher  modal  frequency  yielding  different  complex  medal 
shapes  and  moment  distributions,  more  segments  cf  the  frame  become  diagnesabie. 
The  PNPs  and  corresponding  dimensionless  defections  for  the  frame  excited  at  the 
second  and  third  modal  frequencies  are  summarized  in  Tabie  7. 


TABLE  7 


PNPs  at  the  Second  and  Third  Medal  Frequencies 
(undamaged  frame) 

Excitation  on  column  7  at  .35 

(Girders  19,  22,  27  are  symmetrical  with  IS,  21.  24.) 


At  =  1-72 

girder  no. 

IS  j|  21 

24 

PNP  location 

.S7S5 

.4798 

.5441 

\y  *|x!04  at  PNP 

.52 

Bi 

3.4 

u  -  2.25 

girder  no. 

16  | 

24 

PNP  location 

.6091  | 

.3255  ||  ■  .7174  | 

!  y  *  |  x!0a  at  PNP  j 

1.5 

1 

.72  j 

1.2 

It  will  be  assumed  that  the  response  stations  are  at  the  same  Locations  as  in  the 
previous  section.  Then,  at  the  second  modal  frequency,  the  dimensionless  deflections 
at  the  response  stations  on  girders  21  and  24  are  too  small  to  be  measurable  (Part  1 
of  Table 7)' .  This  is  a  result  of  the  fact  that  the  PNPs  at  the  first  two  modal  frecuencies 
are  close  to  each  ether  for  both  girders.  The  PNP  of  girder  15  at  the  second 
frequency,  however,  is  relatively  far  from  that  at  the  fundamental  frequency.  The 
situation  is  reversed  at  the  third  modal  frequency  (Part  2  of  Table  7) .  It  is  expected 
that  exciting  the  frame  at  the  third  frequency  will  alter  the  medal  shape  mors 
significantly  than  exciting  it  at  he  second  frequency,  thus  making  more  additional 
regions  of  the  frame  diagnesabie. 

Computations  similar  to  those  of  the  previous  section  were  carried  out  at  the 
third  modal  frequency.  Regions  which  become  diagnesabie  in  addition  to  those  in 
Table  6  are  presented  in  Table  g_  Also  shown  are  the  regions  which  may,  at  the 
fundamental  frequency,  be  misdiagnosed  and  are  now  diagnosable  correctly  at  the 
third  frequency.  Those  are  the  regions  0.-0. 14  on  column  7,  0.91*0.92  on  <drder  1c, 


94 


and  0.-0. 05  cn  girder  21.  Table  3  indicates  that  the  general  trend  established  at  the 
fundamental  frequency  Is  preserved  to  a  large  extent  at  the  third  frequency.  That  is. 
a  crack  afiects  the  transmissibility  at  the  closest  response  station  the  mesh 
Exceptions  are  for  the  cases  of  a  crack  on  elements  3  and  24,  in  -which  station  15 
exhibits  the  most  significant  transmissibility  changes  for  the  regions  sho-wn.  rather 

TABLE  3  j 

Additional  Diagnosable  Regions  at  the  Third  Modal  Frequency  | 


All  of  the  parameters  same  as  in  Table  4  except  5* =2.25. 
Nomenclature  same  as  in  Table  5  with  prefix  A  denoting  additional 


|  DJD 

ADR 

RS 

A?B(%) 

j  1 

.47 -.71 

!  is 

i 

.29  -  .47 

15 

I 

1  2 

.55  - .73 

15 

43 

I 

j 

.78  -  .35 

IS 

1 

j  3 

.04-  .51 

16 

47  | 

4  ■ 

i  ’ 

.13  -  .29 
.47  -  .54 

IS 

15 

•  33 

1 

|  5  i 

.38  -  .47 
.51  -  .59 

IS 

IS 

17  ! 

6  ; 

.25  -  .41 

21 

15  j 

i  7  i 

.00  -  .14* 

i 

1 7  j 

1  .  1 

.53  -  .20 

! 

| 

q 

.20  -  .29 

15,19 

l  • 

.56  -  .58 

21,22 

ii 

9 

.25  -  .47 

21,22 

22  J 

.23  -  .25 

15 

1 

18 

.51  -  .52 
.79 -.21 

15 

16 

s  1 

I 

.91  -  .92* 

18 

21 

.00  -  .03* 

21 

...  | 

mm 

.15  -  .25 

IS 

10  | 

17,20,25  | 

APD  less  than  1%  { 

•  (not  included  in  APD) 


MA 


i 


than  station  24.  This  is  not  surprising  since  station  16  is  much  closer  to  the  ?N?  of 
girder  IS  at  the  third  frequency  than  station  24  is  to  that  of  girder  24.  On  the  other 
hand,  it  is  expected  that  more  additional  regions  on  girder  20,  for  example,  "would  be 
diagnosabie  if  station  21  -were  closer  to  the  PN?  of  girder  21  at  the  third  frequency. 

With  the  addition  of  the  new  diagnosabie  regions,  the  diagnosability  for  the  whole 
frame  now  becomes  79%,  with  the  diagnosability  of  the  beam  elements  7,  9,  9.  17,  IS, 
20,  23,  25,  and  2S  remaining  below  the  average.  It  is  also  expected  that  if  the  criterion 
stated  in  Section  4  is  fully  obeyed  in  the  selection  of  the  beam  elements  for  the 
response  stations,  the  diagnosability  of  the  structure  will  increase.  According  to  that 
criterion,  girders  19,  23,  and  25  would  carry  the  response  stations  rather  than  IS,  22, 
and  27,  respectively. 
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The  general  trends  established  above  are  presented  in  the  form  of  "influence 
regions"  in  Figure  29.  Influence  region  of  a  response  station  is  defined  as  the  region  in 
•whose  diagnosability  the  station  in  question  piays  the  most  influential  role.  In  other 
words,  a  crack  in  that  region  aflects  the  transraissibility  at  that  particular  station  the 
most. 

4.7.  The  Effect  of  Crack  Severity  on  the  Trends  • 

The  previous  computations  were  carried  out  for  a  crack  of  sensitivity  number 
9=0.024  (Eq.  2.4C).  For  beams  with  slenderness  ratios  (length  to  beam-depth  ratio)  of 
10  and  20  and  for  a  Poisson's  ratio  of  0.3,  this  9  corresponds  to  the  relative  crack 
depths  of  0.3  and  0.4,  respectively.  Numerical  results  show  that  changing  9  does  not 
have  significant  effects  on  the  trends.  Only  the  transmissibility  changes  are  reduced 
or  amplified  depending  on  whether  9  is  decreased  or  increased,  respectively.  As  9 
gets  smaller  (or  as  the  crack  gets  shallower),  however,  the  diagnosable  regions  of  the 
beam'  elements  are  squeezed  toward  the  frame  joints,  thereby  reducing  the 
diagnosability  of  the  frame.  To  illustrate  the  effect,  diagnosable  regions  on  beam 
elements  3,  5,  20,  and  24  at  the  fundamental  frequency  are  computed  for  a  crack  of 
0=0.005  and  presented  la  Table  9.  This  0  corresponds  to  the  relative  crack  depths  cf 
0.14  an  0.20  for  beams  with  slenderness  ratios  of  10  and  20,  respectively. 

Two  different  sets  of  response  locations  are  used  in  generating  Table  9.  One 
set(group  1  cf  the  table)  is  the  same  as  that  used  previously.  In  the  second  set(group 
2),  the  response  stations  are  moved  closer  to  the  PNPs.  In  both  cases,  diagnosable 
regions  of  the  beam  elements  studied  are  smaller  than  these  for  a  crack  with 
0=0.024.  However,  when  the  response  stations  move  closer  to  the  PNPs,  the 
diagnosability  increases  significantly.  For  a  crack  with  0=0.001,  corresponding  to  the 
relative  crack  depths  of  0.062  and  0.09S  for  slenderness  ratios  of  10  and  20,  no 
segments  of  the  elements  3,  5,  20,  or  24  are  diagnosable  when  the  response  stations 
are  at  the  locations  given  in  Table  9.  In  conclusion,  if  smaller  deflections  can  be 
placed  closer  to  the  PNPs,  larger  regions  become  diagnosable. 


TABLE  9 

Diagnosabie  Regions  for  a  Crack  with  0=.GO5 
(Nomenclature  same  as  in  Table  S) 


l)  Parameter  values  same  as  in  Table  6  except  for  0 


IND 

DR 

RS 

?D(J5) 

3 

.SO  -  1.0 

24 

10 

S 

.73  -  1.0 

24 

22 

20 

.94-1.0 

21 

18 

H 

.92  -  1.0 

24 

8 

2)  Parameter  values  same  as  in  (l)  except  for  the  response  locations 


*ia  -  -813 
=  -287 


*21  =  -4S3 
*22  -  .517 


*2+  =  -634 
*27  =  .356 


DR 

RS  j 

.70  -  1.0 

24 

.00  -  .12 

21 

.57  -  1.0 

24 

.00  -  .09 

21,24 

.72  -  1.0 

21 

.00  -  .09 

24 

.68  -  .74 

24 

.88  -  1.0 

24 

?D(%) 


3 


4.B.  The  Effect  of  Damping  on  the  Trends 


As  the  damping  factor,  o,  is  increased,  the  amplitudes  of  deflections,  jy*j,  and, 
therefore,  the  transmissibilities  are  reduced.  Table  10  lists  four  groups  of  relative- 
transmissibility-change  data  for  a  o’  value  of  0.10.  Other  parameters  are  kept  the 
same  as  in  Table  4  which  contains  data  for  <5=0.01.  The  comparison  of  Table  10  with 
the  corresponding  data  groups  in  Table  4  yields  the  following  conclusions: 


i)  Positive  Rf  values  are  invariably  reduced  as  <5  increases.  An  example  is  the 
case  of  a  crack  on  girder  20.  Although  all  cf  the  responses  are  influenced  in  this  case, 
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TABLE  10 


Relative  Transmissibility  Changes  vs  Crack  Location 


Excitation  on  column  7  at  a  distance  of  .35,  fa  =  .96,  5  =  .10,  0  =  .024. 


Response  Locations : 

=  .SIS  =  -475  «  .649 

^is  ”  .382  £22  ~  .525  £27  =  .351 


(£c  —  crack  location,  RS  :  response  station) 


1)  Crack  on  column  3 


RS  16 

RS  19 

RS  21 

RS  22 

RS  24 

RS  27 

0.001 

-0.037 

0.001 

0.036 

0.008 

BBS 

-0.001 

0.100 

-0.021 

0.000 

0.021 

0.003 

0.049 

-0.001 

0.200 

-0.009 

0.000 

0.009 

0.001 

0.008 

-0.000 

0.300 

0.000 

0.000 

-0.000 

-0.000 

0.000 

0.400 

0.006 

0.000 

-0.006 

-0.C01 

. 

0.024 

0.001 

0.500 

0.009 

0.000 

-0.010 

-0.001 

0.001 

0.600 

0.009 

0.001 

-0.011 

-0.002 

•See] 

0.C02 

0.700 

0.006 

0.001 

-0.011 

-0.001 

0.240 

0.003 

0.300 

0.002 

0.001 

-0.003 

-0.001 

0.345 

0.CC4 

0.900 

-0.004 

0.002 

-0.004 

-0.000 

0.457 

0.004 

0.599 

-0.011 

0.002 

0.001 

0.000 

0.571 

0.005 

2)  Crack  on  column  6 

£c 

RS  16 

RS  19 

RS  21 

RS  22 

RS  24 

RS  27 

0.001 

0.074 

0.007 

K3.429 

-0.044 

-0.325 

-0.C08 

0.100 

0.046 

0.004 

-0.2S0 

-0.030 

-0.153 

-0.C03 

0.200 

0.027 

0.002 

-0.156 

-0.013 

■0.042 

-o.coo 

0.300 

0.011 

0.001 

hO.  065 

-0.008 

0.008 

0.001 

0.400 

0.000 

0.000 

-0.001 

-0.000 

0.001 

0.000 

0.500 

-0.0C5 

-0.000 

0.034 

- 

0.005 

-0.061 

-0.002 

0.600 

-0.006 

0.000 

0.042 

0.009 

-0.171 

,-0.006 

0.700 

-0.003 

0.001 

0.024 

0.010 

-0.325 

-0.010 

o.aoo 

0.005 

0.002 

-0.015 

0.010 

-0.513 

-0.018 

0.900 

0.015 

0.004 

-0.075 

0.008 

-0.725 

-0.023 

0.999 

0.029 

0.008 

-0.150 

0.004 

-0.883 

-0.C30 

TABLE  10 
(continued; 
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1 

3)  Crack  on  girder  20  j 

& 

RS  IS 

RS  19 

RS  21 

■  RS  22 

RS  24 

RS  27 

0.001 

0.115 

0.010 

0.306 

0.052 

ItTSl? 

0.011 

0.100 

0.0S9 

0.C07 

0.155 

0.C2S 

0.166 

0.003 

0.200 

0.QS5 

0.C05 

0.045 

0.010 

0.121 

0.CQ6 

0.300 

0.043 

0.003 

-0.021 

-0.001 

0.060 

0.004 

0.400 

0.C23 

0.001 

-0.041 

-0.005 

0.042 

0.002 

0.500 

0.005 

0 .000 

-0.016 

-0.002 

0.0C9 

0.000 

0.S00 

-0.011 

-0.000 

0.055 

0.007 

-0.020 

-0.001 

0.700 

-0.023 

-0.000 

0.171 

0.024 

- 

0.045 

-0.002 

0.800 

-0.C34 

0.001 

0.330 

0.047 

-O.0S4 

•  -0.003 

0.900 

-0.041 

.  0.002 

0.531 

0.076 

-0.079 

'-0.003 

0.999 

-0.046 

0.004 

0.788 

0.110 

-0.089 

-0.003 

4)  -Crack  on  girder  24 


£c 

RS  16 

RS  19 
« 

RS  21 

RS  22 

RS  27 

0.001 

-0.011 

0.002 

0.C01 

0.000 

Kg 

0.005  ] 

0.100 

-0.009 

0.002 

0.C05 

0.001 

0.C02  j 

0.200 

-0.005 

0.001 

0.006 

0.002 

-0.000  1 

0.300 

-C.C04 

0.000 

0.008 

0.002 

frMEEl 

-0.001 

0.400 

-0.002 

0.000 

0.006 

0.0C2 

C'Btil 

-0.001 

0.500 

-0.001 

-0.000 

0.002 

0.001  ■ 

-0.075 

-0.001 

0.600 

•  0.001 

Q.Q00 

-0.004 

-0.001 

0.149 

0.001 

a  700 

a  002 

0.000 

-0.003 

0.323 

0.004 

0.800 

0.003 

0.001 

-0.022 

-0.0C6 

0.146 

0.C07 

0.900 

0.004 

0.002 

-0.034 

-0.009 

-0.313 

0.011 

0.999 

0.005 

0.003 

-0.047 

-0.012 

-0.972 

0.016 

tie  most  adversely  affected  transmissihility  changes  are  those  at  the  stations  IS  and 
27,  the  farthest  two  stations  from  the  cracked  girder,  which  were  reduced  to  near 
zero  values.  The  stations  near  the  crack  are  not  affected  so  significantly.  Example: 
the  case  of  a  crack  on  girder  24. 

ii)  Negative  P.7  values  are  affected  less  by  an  increase  in  <5.  The  magnitude  of  a 
negative  transmissihility  change  can  be  amplified  or  reduced,  both  not  very 
significantly.  Example:  P7  values  at  station  24  for  a  crack  cn  beam  elements  8  and 
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20. 

The  diagnosable  regions  with  5=0.10  are  computed  and  presented  in  Table  13. 
When  these  regions  are  compared  with  the  corresponding  data  in  Table  6 ,  it  is 
observed  that  an  increase  in  the  damping  factor  from  0.01  to  0.10  leaves  the  trends 

A 

essentially  unaffected  and  reduces  the  diagnosable  regions  slightly  (between  0  to  6 
percentage  points  for  the  cases  presented  in  Table  n  j . 

4.9.  A  Response  Station  on  a  Column 

At  the  fundamental  frequency,  the  transmissibility  change  on  a  column  for  a 
given  crack  is  found  to  be  almost  independent  of  the  response  location  on  that 
column.  Relative  transmissibility  changes  at  a  response  station  on  column  8  for  a 
crack  on  various  beam  elements  are  shown  in  Table  12  .In  this  examole,  the  excitation 
is  on  column  1.  Only  one  crack  location  on  each  beam  is  considered.  Beams  are 
grouped  based  on  the  closeness  of  the  corresponding  Rr  values.  For  example.  R ? 
values  for  a  crack  at  a  location  of  0.15  cn  each  of  the  beams  1.  7,  13  and  16* range 
from  0. 125  to  0. 150.  The  R?  value  for  a  particular  crack  is  nearly  the  same  regardless 
of  where  the  response  station  is  located  on  column  3.  That  nearly  constant  value  for 
the  crack  on  column  7,  for  example,  is  0.15.  Furthermore,  the  transmissibility 
changes  at  a  station  on  column  9  are  very  close  to  those  on  column  8,  as  well  as  being 
nearly  constant  along  the  length  of  column  9. 

Table  12  Indicates  that  as  the  crack  is  moved  to  upper  levels  on  the  frame,  the 
Rt  value  on  column  8  decreases.  There  is  especially  a  significant  decrease  from  the 
first  level  to  the  second(Fig.  17) .  The  overall  Rf  values,  however,  are  not  very  large. 
A  transducer  on  a  column  may,  nevertheless,  helD  in  diagnosis  by  differentiating 
between  a  crack  on  the  first  level  and  on  the  upper  levels. 

4.10.  Effect  of  the  Crack  Model  on  the  *TVwr)grm‘«wri'hiHty  Method 

In  this  study,  a  Mode-I  crack  has  been  modeled  as  a  torsional  spring.  The  crack 
surface  nas  been  assumed  to  be  perpendicular  to  the  neutral  of  the  beam 
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TABLE  li 

Diagnosable  Regions  for  a  Frame  with  Heavier  Damping 
Parameter  values  same  as  in  Table  S  except  <5  =  .  10 
Nomenclature  same  as  in  Table  8. 


3 

.00 -.03 
.54-1.0  ■ 

8 

.00 -.28 
.54-  1.0 

20 

.00 -.15 
.15 -.25 
.64-1.0 

24 

.00 -.14 
.25  -  .48 
.58  -  .81 
.85-1.0 

PD(JS) 
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TABLE  12 

Relative  Trans  missibility  Changes  at  a  Response  Station  on  Column  9 
<5=.01,  excitation  on  column  1  at  .95,  .9*  =  .96 
0  =  .024,  crack  location  is  .15  on  all  the  beam  elements 
IND  =  no.  of  the  beam  element  ■with  the  crack 


1 

IND  ! 

t 

1.7,13,16 

5,11 

20,21.23 

6,9 

24.25, 2S 

15 

RT  \ 

\ 

.125 -.150 

.069 

.034-  .060 

.010-  .011 

.004 -.008 

.002 

element  on  which  the  crack  is  located.  It  can  be  shown,  however,  that  the  trends  o; 
relative  trans missibility  changes  are  independent  of  the  crack  model  as  long  as  the 
crack  introduces  a  reduction  in  the  flexural  rigidity  of  the  beam.  To  this  end,  a  one- 
story  four-span  frame  excited  at  the  mid-column  is  analyzed  (Figure  30') .  The 
fundamental  characteristic  value  of  such  a  uniform  frame  is  LS87.  In  the 


computations,  the  frame  is  excited  at  a  ft*  value  of  1.58.  Damping  factor  is  taken  to 
be  0.01.  Two  response  stations  are  placed  on  girders  S  and  3  at  the  locations  0.63  and 
0.56,  respectively.  (?NFs  of  these  girders  for  a  uniform  frame  are  at  0.681  and  0.551, 
respectively,  when  the  excitation  is  on  column  3  at  a  distance  of  0.95  from  the 
ground.) 

First,  the  relative  transmissibility  changes  at  the  two  stations  of  the  uniform 
frame  were  computed  for  a  crack  of  ©=0.01  on  girder  9,  as  the  crack  location  was 
varied  along  the  girder.  It  was  observed  that  the  station  on  girder  3  exhibited  much 
more  significant  transmissibility  changes  except  when  the  crack  was  in  the  mid¬ 
region  of  girder  9.  In  the  latter  case,  both  transmissibility  changes  were  small.  Next, 
the  length  of  girder  9  was  assumed  to  be  20%  longer  than  the  lengths  of  the  other 
columns  and  girders,  which  made  the  former  less  rigid  than  the  others  and 
qualitatively  simulated  a  crack.  With  the  other  parameters  kept  the  sane,  the 
transmissibility  changes  at  the  two  stations  relative  to  the  transmissibilities  across 
the  undamaged  uniform  frame  were  computed.  The  relative  transmissibility  changes 
due  to  the  lengthening  of  girder  9  were  found  to  be  -0.41  and  0.39  at  stations  3  and  8, 
respectively.  As  the  result  indicates,  softening  in  a  region  of  a  frame  structure 
afiects  the  response  at  the  nearest  station  the  most. 
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Figure  30 .  A  one-story  four-span  frame. 
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CHAPTER  5 

Discussion  and  Conclusions 


Nondestructive  testing  of  large-scale  structures  is  still  in  its  infancy.  Local 
methods  currently  used  require  a  sweeping  of  the  wfaoLs  structure  and  can  become 
costly  and  time-consuming.  A  global  method  whereby  a  structure  can  be  continually 
monitored  and  damage  can  be  isolated  at  least  within  a  section  of  the  structure  is 
needed-  Such  a  method,  the  transmissibiiity  method,  is  proposed  La  this  study.  The 
transmissibility  method,  which  is  based  on  the  relative  transmissibiiity  changes  at 
various  response  stations,  isolates  a  crack  to  within  one  or  two  ceils  of  a  frame 
structure  after  which  a  local  method  may  be  used  to  pinpoint  the  crack.  In  the  study, 
the  effects  of  excitation  Location  and  frequency,  location  of  response  stations,  crack 
Location  and  severity,  and  damping  on  the  relative  transmissibiiity  changes  are 
investigated.  General  trends  of  the  values  of  transmissibiiity  changes  as  a  function  of 
these  factors  are  established.  In  practice,  these  trends  can  be  used  fcr  locating  a 
crack,  given  the  transmissibiiity  changes.  The  following  conclusions  are  reached  from 
the  transmissibiiity  analysis: 

i)  The  most  significant  transmissibiiity  change  is  usually  observed  at  a  response 
station  in  the  neighborhood  of  the  crack. 

ii)  The  trends  as  stated  in  (i)  are  highiy  insensitive  to  the  excitation  location. 

ill)  The  trends  are  insensitive  to  the  crack  depth  and  material  properties  such  as 
damping. 

iv)  The  trends  are  highiy  insensitive  to  the  location  of  the  response  stations  when 
the  stations  are  hi  the  vicinity  of  the  PNPs.  Trends,  however,  are  lost  away  from  the 
P NPs,  in  which  case  damage  may  be  detected  but  cannot  be  located.  On  the  other 
hand,  the  closer  the  response  stations  to  the  PNPs,  the  larger  are  the  transmissibiiity 
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changes  and  the  chances  of  diagnosing  a  crack. 

v)  Amplitude  of  deflections  and  accelerations  are  small  near  the  ?N?s.  Therefore, 
how  close  the  response  stations  can  be  located  to  the  PNPs  depends  on  the  signal 
levels  which  can  be  measured. 

vi)  Cracks  at  or  near  the  joints  of  a  frame  are  most  easily  diagncsable.  A  crack 
located  away  from  the  joints,  in  the  mid-regions  of  the  beam  elements,  may  not  be 
diagnosable. 

vii)  Diagnos  ability  of  a  crack  increases  if  transmissibility  data  is  available  at 
more  than  one  resonant  frequency. 

Conclusion  (vi)  is  especially  encouraging  in  the  Light  of  the  fact  that  structures 
such  as  offshore  platforms  usually  fail  at  the  foundation  or  develop  cracks  at  the 
welded  joints. 

The  electrical  analogy  method  as  developed  for  frame  structures  in  this  study 
provides  an  economical  test  for  the  analysis  of  large  structures.  Tne  responses 
computed  with  the  analogy  are  exact  Insofar  as  3emoulli-Euler  theory  is  exact.  Tne 
structures  analysed  in  the  study  are  assumed  rigidly  fixed  at  the  ground  level  In  the 
case  of  non-rigid  foundations,  foundation  flexibility  can  be  accounted  for  by  inserting 
resistances  in  the  electrical  circuit  at  points  corresponding  to  the  anchors[34,  pp. 
82-95]. 

The  other  nondestructive  method  of  damage  diagnosis  investigated  in  this  study 
is  the  frequency  method.  In  complex  structures,  the  changes  in  modal  frequencies 
due  to  a  crack  are  small.  Therefore,  the  frequency  method  is  not  suitable  for  such 
structures.  For  simpLe  structures  such  as  beams,  however,  it  is  a  feasible  method 

It  is  shown  in  the  study  that  closely  spaced  multiple  cracks  can  be  effectively 
represented  by  a  single  crack  for  which  the  sensitivity  number,  9,  is  approximately 
equal  to  the  sum  of  the  individual  sensitivity  numbers.  Uncertainty  exists,  however,  in 
diagnosing  such  a  damage  as  to  whether  there  is  only  one  major  crack  or  several 
closely  spaced  minor  cracks.  Nevertheless,  iocadon  of  the  damage  can  be  identified 
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quite  accurately. 

In  the  case  of  a  structure  with  several  cracks,  only  cne  of  which  is  severe,  it  is 
possible  to  diagnose  the  major  crack.  The  contributions  of  minor  cracks  to  frequency 
decreases  is  small  compared  to  that  of  the  severe  crack.  The  relative  frequency 
changes,  f?3-,  in  such  a  situation,  exhibit  a  trend  similar  to  that  which  would  be 
observed  were  the  major  crack  the  only  crack  present  on  the  structure.  If  there  is 
more  than  one  major  crack,  not  closely  spaced,  on  the  structure,  it  is,  in  general,  not 
possible  to  diagnose  deterministically  the  damage  with  only  a  few  values  known. 

In  the  frequency  method,  curves  of  relative  frequency  changes,  versus  crack 
location  yield  valuable  information  regarding  damage  diagnosis.  It  may  be  possible  to 
estimate  damage  location,  provided  there  is  only  one  major  crack,  by  only  looking  at 
the  relative  magnitudes  of  frequency  decreases  and  using  the  relative-change  curves. 

In  diagnosing  damage  in  a  structure  using  the  frequency  method,  given  three 
post-damage  characteristic  values,  £,  the  solution  is  sought  for  a  single  crack  since 
there  are  not  enough  equations  to  solve  for  a  large  number  of  cracks.  In  solving  for  a 
single  crack  (L.e.  for  ef<7  and  two  measured  j3  values  are  actually  used  hrst  The 
third  measured  £  value  is  then  compared  with  the  computed  counterpart  which  would 
be  produced  by  the  crack  with  the  characteristic  pair  ,  S#7).  It  is  shewn  that  the 
two J3  values  may  not  match  closely  ever,  though  the  diagnosis  is  accurate  (that  is,  ea? 
and  9t?  identify  the  major  crack  correctly),  the  discrepancy  possibly  being  due  to  the 
presence  of  ether  minor  cracks.  If  the  discrepancy  is  large,  then  the  solution  is 
rejected  with  the  conclusion  that  :here  is  more  than  cne  major  crack.  Cn  the  ether 
hand,  the  choice  of  the  pair  of  measured  0  values  to  be  used  in  computing  ett}  and 
may  affect  the  diagnosis.  In  some  cases,  choosing  the  two  characteristic  values  which 
correspond  to  the  largest  two  J?j  values  leads  to  the  correct  diagnosis.  It  is  not  clear, 
however,  whether  there  is  a  Tight  choice  in  each  case. 
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APPENDIX 


The  Program  to  Compute  Dimensionless  Deflections  in  a  Frame 
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SUBROUTINE  T?JL\B(N.M.K.ES,EC,THE.IND,IE,XE.IT,XT,NT,EUN,DAM. 
UXXmZ.RH.IUR.IXR.WA.'ffK.IWK) 

C 

COMPLEX  "IS  U(IURIUR).  XTR(IURIXR),  Z(JXHJXR).  X(K3.IUR). 

•  RH(IXR).  WA(l).  DCMPLX.  CD  SIN,  CDCOS,  CDEX?, 

•  B,  BI.  W.  WP.  Q,  QP,  P.  PP.  EXPl.  EXP2,  ALF 

REAL-9  XT(NT).  ES(i).  EC(l).  THE(l).  WK(l),  XE,  BUN,  DAM.  Dl.  D2,  D3,  D4  . 
INTEGER  IT(NT).  IWK(l),  IND(l),  N.  M.  X.  IE.  NT.  IUR  IXR 
C 

C  THIS  SUBROUTINE  COMPUTES  THE  DIMENSIONLESS  COMPLEX  DEFLECTIONS 
C  AT  SELECTED  LOCATIONS  ON  AN  NxM  PLANAR  FRAME  STRUCTURE  "WITH  OR 
C  WITHOUT  CRACKS  EXCITED  AT  AN  ARBITRARY  LOCATION  BY  A  SINUSOIDAL 
C  FORCE.  THE  FRAME  IS  ASSUMED  TO  HAVE  UNIFORM  PROPERTIES  EXCEPT  ' 

C  FOR  THE  FLOOR  HEIGHTS  AND  SPAN  WIDTHS.  THE  BEAM  ELEMENTS 
C  COMPOSING  THE  FRAME  ARE  NUMBERED  AS  EXPLAINED  IN  SECTION  2.4. 

C  THE  FRAME  IS  'SIMULATED  BY  II  CIRCUITS. 

C 

C  WRITTEN  BY  MEKMET  A.  ARGUN 
C  DATE  LAST  REVISED  :  APRIL  19B5 
C 

C  INPUT  :  (ALL  OF  THE  INPUT  VARLASLES  ARE  PRESERVED) 

C 

C  N  NUMBER  OF  STORIES  OF  THE  FRAME 
C  M  NUMBER  OF  SPANS  OF  THE  FRAME  (NO.  OF  ANCHORS  -  1). 

C  ■  X  NUMBER  OF  CRACKS  ON  THE  FRAME  (ZERO  FOR  NO-CRACK  CASE). 

C  ES  A  VECTOR  OF  LENGTH  N+M.  THE  FIRST  N  LOCATIONS  SHOULD  CONTAIN 
C  THE  RELATIVE  FLOOR  HEIGHTS,  Q/Lc.  IN  THE  ORDER  FROM  GROUND 

C  UP,  WHERE  Lc  IS  THE  REFERENCE  LENGTH  CHOSEN  FOR  THE  FRAME  AND 

C  Li  IS  THE  HEIGHT  OF  THE  ith  FLOOR  THE  LAST  M  LOCATIONS 

C  SHOULD  CONTAIN  THE  RELATIVE  SPAN  WIDTHS,  Q /  Lq,  IN  THE  ORDER 

C  FROM  LEFT  TO  RIGHT.  WHERE  Q  IS  THE  WIDTH  OF  THE  j  th  SPAN. 

C  EC  A  VECTOR  OF  LENGTH  MAX(  l.K)  CONTAINING  THE  RELATIVE  CRACK 
C  LOCATIONS,  LQ/ Q,  WHERE  Lj  IS  THE  LENGTH  OF  THE  ;  th 

C  BEAM  ELEMENT  WHICH  CARRIES  THE  ith  CRACK.  WHEN  THE  ith  CRACK 

C  •  IS  THE  ONLY  CRACK  OR  THE  FIRST  ONE  OF  A  GROUP  OF  CRACKS  CN 
C  THE  j  th  ELEMENT,  THE  DISTANCE  LQ  OFTKS  CRACK  IS  MEASURED 

C  FROM  THE  LEFT(LOWER)  END  OF  THE  GIRDER( COLUMN)  ON  WHICH  THE 
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C  CRACK  IS  LOCATED.  OTHERWISE.  LQ,  IS  TEE  DISTANCE  BETWEEN  THE 

C  (i -l)th  AND  ith  CRACKS.  THE  ORDERING  OF  THE  CRACK  LOCATIONS 

C  WITHIN  EC  SHOULD  FOLLOW  THE  NUMBERING  SCHEME  EXPLAINED  IN 

C  SECTION  2.4.  NAMELY,  FIRST  THE  LOCATIONS  OF  THE  COLUMN  CRACKS. 

C  THEN  THOSE  OF  THE  GIRDER  CRACKS  SHOULD  BE  WRITTEN. 

C  0  <  £C(/)  <  1.  EC  IS  IGNORED  WHEN  K=Q. 

C  .  TEE  A  VECTOR  OF  LENGTH  MAX(l.K)  CONTAINING  THE  SENSITIVITY 
C  NUMBERS.  0,  OF  TEE  CRACKS  IN  THE  SAME  ORDER  AS  EC. 

C  TEE(I)  *  0  for  1  £  I  is  K  IGNORED  WHEN  K=0. 

C  IND  A  VECTOR  OF  LENGTH  MAX(  1. K)  CONTAINING  TEE  NUMBERS  OF  TEE 
C  BEAM  ELEMENTS  ON  WHICH  TEE  CRACKS  ARE  LOCATED.  ORDERED  IN 

C  CORRESPONDENCE  WITH  EC.  IND(l)  £  IND(2)  £  ...  £  IND(K). 

C  IE  TEE  NUMBER  OF  THE  BEAM  ELEMENT  ON  WHICH  THE  EXCITATION  IS ' 

C  APPLIED. 

C  XE  TEE  RELATIVE  LOCATION  OF  TEE  EXCITATION.  LE/L(3E).  WHERE  LE 
C  IS  TEE  DISTANCE  TO  THE  EXCITATION  FROM  THE  LEFT(LOWER)  END 

C  OF  THE  GIRDER(COLUMN)  ON  WHICH  TEE  EXCITATION  IS  APPLIED. 

C  L(IE)  IS  THE  LENGTH  OF  TEAT  GIRD ER( COLUMN) .  0  <  XE  <  i  . 

C  NT  TOTAL  NUMBER  OF  RESPONSE  STATIONS  ON  THE  FRAME. 

C  IT  A  VECTOR  OF  LENGTH  NT  CONTAINING  THE  NUMBERS  OF  THE  BEAM 
C  ELEMENTS  WHICH  HAVE  A  RESPONSE  STATION  ON  THEM.  ONLY  ONE 

C  RESPONSE  STATION  FER  BEAM  ELEMENT  IS  ALLOWED. 

C  IT(1)  <  IT(2)  <  ...  <  IT(NT). 

C  XT  A  VECTOR  OF  LENGTH  NT  CONTAINING  THE  RELATIVE  LOCATIONS. 

C  '  LT(J)/L(IT(J)),  OF  THE  RESPONSE  STATIONS  ON  THE 

C  CORRESPONDING  3EAM  ELEMENTS,  WHERE  L(IT(J))  IS  THE  LENGTH 

C  OF  THE  PARTICULAR  3EAM  ELEMENT. 

C  BUN  UNDAMPED  BETA  VALUE. 

DAMPING  FACTOR 

N“(M+1)  +  2*K  +  NTrL  THE  ORDER  OF  MATRIX  U. 

N  -r  K  +  NT  +  1.  THE  ORDER  OF  MATRIX  Z. 

C 

C  OUTPUT: 

C  7 

C  RH  A  COMPLEX  VECTOR  OF  LENGTH  JXR  CONTAINING  THE  NONDIMENSIONAL 
C  '  TRANSVERSE  DEFLECTIONS,  y  *.  THE  FIRST  N  LOCATIONS  CONTAIN  THE 
C  LATERAL  DISPLACEMENTS  OF  THE  N  FLOORS.  ORDERED  FROM  THE  FIRST 
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C  FLOOR  UP.  THE  NEXT  X  LOCATIONS  CONTAIN  THE  DEFLECTIONS  AT  THE 
C  CRACX3,  ORDERED  IN  ACCORDANCE  WITH  EC.  THE  DEFLECTIONS  AT  THE 

C  RESPONSE  STATIONS  ARE  CONTAINED  IN  THE  NEXT  NT  LOCATIONS.  THE 

C  LAST  LOCATION  CONTAINS  THE  DEFLECTION  AT  THE  EXCITATION  POINT. 

C  THE  CHANGES  IN  TRANSMISSEIUTIES  CAN  BE  OBTAINED  BY  COMPUTING 
C  '  THE  AMPLITUDES  OF  DEFLECTIONS  AT  THE  RESPONSE  STATIONS  WITH 
C  AND  WITHOUT  THE  CRACK(S). 

C 

C  WORK  ARRAYS: 

n 

C  U  A  COMPLEX  ARRAY  OF  ORDER  IUR  X  IUR. 

C  X  A  COMPLEX  ARRAY  OF  ORDER  IXRXIUR. 

C  XTR  A  COMPLEX  ARRAY  OF  ORDER  IUR  X  KP- 
C  Z  A  COMPLEX  ARRAY  OF  ORDER  DCR  X IXR. 

C  WA  A  COMPLEX  VECTOR  OF  LENGTH  IUR»(IUR  +  2). 

C  WX  A  REAL  VECTOR  OF  LENGTH  IUR. 

C  '  IWX  AN  INTEGER  VECTOR  OF  LENGTH  X  +  2»NT  +■  LIHS  FIRST  X  LOCATIONS 
C  ARE  RESERVED  FOR  THE  MULTIPLICITIES  OF  THE  CRACKS  ON  A  BEAM 
C  ELEMENT.  THE  NEXT  NT  +  1  LOCATIONS  SOMETIMES  CONTAIN  CRACK 

C  INDICES  WHEN  A  CRACK  AND  A  RESPONSE  STATION  OR  THE  EXCITATION 

C  ARE  LOCATED  ON  THE  SAME  BEAM  ELEMENT.  THE  NUMBERS  OF  THE 
C  ROWS  IN  Z  MATRIX  WHICH  CORRESPOND  TO  THE  DEFLECTIONS  AT 
C  THE  RESPONSE  STATIONS  .ARE  SAVED  IN  THE  LAST  NT  LOCATIONS. 

C 

C  REQUIRED  IM5L  ROUTINES  : 

C 

C  LSQ2C  (COMPLEX  LINEAR  EQUATION  SOLVER)  AND  OTHER  ROUTINES  CALLED  BY  IT. 
C  LEQ2C  FACTORS  THE  COEFFICIENT  MATRIX  INTO  THE  L-U  DECOMPOSITION 
C  OF  A  ROWWISE  PERMUTATION  OF  A  AND  SOLVES  THE  SYSTEM.  LEQ2C  COMPUTES 
C  THE  RESIDUALS  AND  USES  ITERATIVE  IMPROVEMENT  UNTIL  THE  SOLUTION 
C  IS  ACCURATE  TO  MACHINE  PRECISION.  (REF:  IMSL  LIBRARY  REFERENCE 
C  MANUAL  9TH  ED..  V.2,  CKPT.  L  IMSL  INC..  1932.) 

C 

C  REQUIRED  FORTRAN  FUNCTION  SUBPROGRAMS  : 

C 

C  CD  SIN,  CDCOS.  CD  EXP,  DCMPLX  DSQRT,  DMINl.  DMAXL 
C 


C  ERROR  CONDITIONS  : 

C 

C  IF  IE  =  rr(J)  FOR  ANY  J  AND  A3S(XE  -  X7(J))  IS  SMALLER  THAN 
C  5.D-4  OR  IF  DJD(I)  =  IS  OR  IT(J)  FOR  ANY  I  AND  J  AND 
C  EC(I)  =  XS  OR  XT(J) .  THE  SUBROUTINE  RETURNS  TO  TEE  CALLING  PROGRAM. 

C  •*••••******•***•***•**•■*■***•*»•■■*******•***••»**■»*»•»**•*••»••*••••••*•«»« 

c 

C  COMPUTE  THE  DAMPED  3  VALUE 
C 

Dl  =  DSQR7(l.D0+DAM«2) 

D2  =  .5D0/DSQRT(Dl) 

D3  =  .5D0*DSQRT(.5D0-t-.5D0*Dl)/Dl 
B  =  SUN-DCM?LX(DSQKr(D2+D3).-DSQRT(D2-D3)) 

C 

C  FIND  THE  MULTIPLICITY  OF  CRACKS  ON  EACH  CRACKED  ELEMENT 
C 

IWX(l)  =  1 

IF  (K  .LE.  1)  GO  TO  5 

J3  =  1 

KM1  =  K  -  1 

DO  4  Jr  =1.  KM1 

IF  (IND(Jl-rl)  .NS.  IND(Jl))  GO  TO  3 
rWH(J3)  =  IWK(J3)  +  I 
GO  TO  4 

3  J3  =  J3  +  1 
IWK(J3)  =  1 

4  CONTINUE 
C 

C  CLEAR  THE  MATRICES 
C 

5  DO  8  J1  =  1,  IUR 
DO  6  II  =  1.  IXR 

S.X(IUl)  =  0.D0 
DO  7  II  a  1.  IUR 

7  U(I1.J1)  =  0.D0 

8  CONTINUE 

DO  9  Jl  =  1,  KR 


DO  9  II  =  1.  KR 
9  Z(U.J1)  =  0.D0 
NMP1  =  N-(M-rl) 

C 

C  SET  VARIABLES 
C 

JET  —  1 
KR  =  1 
J  =  1 
J3  =  1 

J4  =  NMP 1  +  1 
JS  =  IWK(l) 

JADD  =  0 
NT?1  =  NT  +  1 
DO  12  L*  1.  NTPl 
12  IWX(K+L)  =  0 
ICOMP  =  4 
NOW  s  1 
LUaM+1 
IU  =  N 
C 

C  ESTA3LISH  THE  MATRICES 
C 

15  DO  230  L  =  1,  LU 

11  =  JADD 
J1  =  L 

DO  220  I  =  1.  FJ 

12  =  1 
J2  =  J 

IF  (ICOMP  .SQ.  0)  J2  =  I*N  +  L 
IGIRD  =  1  ■ 

D3  =  O.DO 

IF  (IND(XR)  .NE.  J  .OR.  XR  .GT.  K)  GO  TO  120 
Dl  =  EC(KR) 

IF  (IE  .NE.  J)  GO  TO  20 
D2  =  XS 

IF  (NOW  .GT.  NT  .OR  IE  .NE.  IT(NOW))  GO  TO  40 


IF  (DABS(XE  -  XT(NOW))  .LT.  5.D-4)  GO  TO  240 
D2  =  D3flNlCXKXT(NOW)) 

D3  =  D  1LAX  1  (XS.XT  ( N  QW) ) 

GO  TO  40 

20  IF  (NOW  .GT.  NT. OR.  IT{NOW)  .NE.  J)  GO  TO  SO 
D2  =  XT(NOW) 

40  IF  (D1  -  D2)  50.  240,  70 

50  BI  =  EC(KR)  *3 
12  S  N  +  KR 
J2  =  J4 
J4  =  J4  4  2 

U(J2,J2)  =  -l.D0/(3*THE(KR)*DCMPLX(l.D0,DAM)) 
U(J2,J24l)  s  -U(J2,J2) 

U(J2+1.J2)  =  U(J2.J24i) 

U(J2+1,J2+1)  a  U(J2.J2) 

KR  =  XR  +  1 
ASSIGN  ISO  TO  JUMP 
ASSIGN  40  TO  JCR 
IF  (KR  .LE.  J0)  GO  TO  60 
J3  =  J3  4  1 
ja  =  iwx(J3) 

ASSIGN  80  TO  JCR 
GO  TO  130 

60  D1  a  D1  4  EC(KR) 

GO  TO  130 

70  31  =  (EC(KR)  4  D2  -  Dl)’S 
WK(JET)  =  EC(KR) 

IWK(X4JET)  =  KR 
'  EC(XR)  =  Dl  -D2 
J2  =  NMPl  4  K  4  K  4  JET 
12  a  N  4  K  4  JET 
JET  =  JET  4  l 
ASSIGN  200  TO  JUMP 
04  =  D2 

IF  (D3  .EQ.  O.EO)  GO  TO  71 
ASSIGN  40  TO  JCR 
02  =  D3 
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D3  =  O.DO 
GO  TO  72 

71  ASSIGN  SO  TO  ICR 

72  IF  (IT(NCW)  .NS.  IND(XR)  .OR.  D4  .NS.  XT(NOW))  GO  TO  73 
IWK(X+NTPl+NOW)  =  12 

NOW  =  NOW  +  1 
GO  TO  130 

73  17  =  12 

GO  TO  130 
SO  31  a  (02  -  Dl)*3 
12  =  N  +  K  +  1ST 
J2  *  NMPl  +  K  +  K  +  IET 
1ST  =  JS7  -r  1 
Dl  =  D2 

ASSIGN  200  TO  JUMP 
ASSIGN  110  TO  JCR 
IF  (03  .Sq.  O.DO)  GO  TO  84 
02  a  C3 
03  a  0.00 
ASSIGN  90  TO  JCR 

94  IF  (IT(NOW)  .NS.  I  .OR.  Dl  .NS.  XT(NCW))  GO  TO  85 
IWK(X+  NTP 1 + N  OW)  =  12 

NOW  =  NOV f  +  1 
GO  TO  130 

95  17  a  12 

GO  TO  130 
90  31  a  SC(KR)  "S 
12  =  N  *  KR 
J2  a  J4 
J4  a  J4  +  2 

U(J2J2)  a  -  1.D 0/(3 *THS(KR)  *DCMPLX(  l.DO.D AM) ) 
U(J2,J2+1)  a  -U(J2,J2) 

17(12+1.12)  =  U (12,12+1) 

U(J2+1.J2+1)  =  U(J2.J2) 

KR  a  KR  +  1 
ASSIGN  90  TO  JCR 
ASSIGN  190  TO  JUNG3 
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IF  (XR  .LE.  JB)  GO  TO  100 
J3  =  J3  +  1 
J8  =  IWX(J3) 

ASSIGN  110  TO  ICR 
GO  TO  130 

ICO  Dl  =  Dl  +  EC(KR) 

GO  TO  130 

110  BI  =  (1.D0  -  Dl)*3 
IGIFJ3  =  2 
12  =  1 
*2  = ; 

IF  (ICOM?  .EQ.  0)  J2  =  I-N  +  L 
ASSIGN  210  TO  JUMP 
GO  TO  130 

120  Dl  =  0.D0 

IF  (IE  .NE.  J)  GO  TO  121 
D2  =  XE 

IF  (NCVT  .GT.  NT  .OR.  IE  .NE.  IT(NOW))  GO  TO  SO 
IF  (DA3S(XZ  -  XT(NO’.V))  .17.  5.D-4)  GO  TO  240 
D2  =  D'GN  1(XE.XT(NC’V)) 

03  =  0  MAX  1  (XE.  XT  ( N  C  W) ) 

GO  TO  SO 

121  IF  (NC’V  .GT.  NT  .CP.  IT(NOW)  .NE.  J)  GO  TO  122 
02  a  XT(NCW) 

GO  TO  SO 

122  31  =  3 
igird  =  3 

ASSIGN  210  TO  JUMP 
120  31  =  ES(WAOD)*5I 
C 

C  COMPUTE  THE  H-CIRCUIT  PAPJLVETSRS 
C 

EXP1  =  C0EXP(3I) 

EXP2  =  CDEX?(-3I) 

W  =  (SXP1  -  EX?2)/2.C0 
T?  =  (EXPl  +  EX?2)/2.DQ 
EXP1  =  CDSIN(BI) 
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EXP2  =  CECGS(BI) 

ALF  =  WP*EX?2  *  1.C0 
P  =  VP-SXPl 
Qa(P- WSX?2)/ALF 
p  =  (P  +  W*EXP2)/ALF 
QP  =  WEXPl/ALF 
pp  =  (-W  +  EX?l)/ALF 
W  =  (W-EXPl)/ALF 
W?  =  (WP  -  EX?  2)/ ALF 
C 

U(J2.J2)  =  U(J2,J2)  +  Q 
IF  (II  .EQ.  0)  GO  TO  160 
U(JUl)  =  U(J1.J1)  +  Q 
U(J1.J2)  =  W 
U(J2.J1)  =  W 
13  =  ICCMP  +  IGIRD 
GO  TO  (150, 170, 130), 13 
Z(IUl)  =  2(11.11)  +  ? 

Z(U.I2)  =  2(11.12)  -  ?? 

Z(I2,U)  =  Z (12,11)  -  ?? 

X(1UD  =  X(U.Jl)  -  Q? 

X(U.J2)  =  -WP 
150  X(I2, J 1)  =  W? 

ISO  X(I2,J2)  =  Q? 

Z(I2.I2)  =  Z(I2J2)  +  ? 

GO  TO  150 

170  X(U.J1)  -  X(U.Jl)  -  QP 
X(1U2)  =  -VP 
Z(I1.I1)  =  2(11.11)  +  ? 

ISO  J1  =  J2 

11  =12 

GO  TO  JUMP. (30,  110.  ISO.  200.  210) 
7  *190  J1  *  J1  +  1 

200  IGIRD  =  4 

GO  TO  JCR,(40,  SO.  90,  110) 

210  J  =  J  +  1 
220  CONTINUE 


A. 
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220  CONTINUE 

IF  (JADD  .EQ.  N)  GO  TO  2S0 
JADD  =  N 
ICOM?  a  0 
LU  =  N 

iu=  m 

GO  TO  15 
240  WRITE(6,25Q) 

250  FORMAT( IX,' LOCATIONS  SPECIFIED  FOR  DIFFERENT  THINGS  TOO  CLOSE’) 
RETURN 
2S0  Dls  0.D0 

DO  2S5  I  =  1.  M 
2S5  Dl  =  Dl  +  ES(N-rl) 

El  a  B-Dl 
DO  2S6  I  a  l.  N 
2SS  2(1.1)  =  Z(I.I)  -  BI 
C 

C  MATRICES  ESTABLISHED 
C 

DO  270  I  a  l,  NTPl 

IF  (IWK(K+N7?1+1-I)  .EQ.  0)  GO  TO  270 
EC(IWK(X+NT?l-rl-I))  =  WK(NT?1+1-I) 

270  CONTINUE 

DO  300  J  a  i.  IXR 
DO  290  I  a  l,  IUR 
290  XTR(I.J)  =  X(J,I) 

300  CONTINUE 
C 

C  COMPUTE  ITlX 
C 

IJOB  =  0 

CALL  LEQ2C(U,  IUR,  IUR,  XTR,  E<R,  IUR,  UCB,  WA,  T(VK,  ER) 

IF  (ER  .NE.  0)  RETURN 
C 

C  COMPUTE  Z-Xl Tl2f 
C 


DO  380  J  =  1,  IXR 


[ATE 
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